Search

Abstract Orobanche coerulescens is a parasitic plant that cannot complete its life cycle without a host and is incapable of photosynthesis. The habitats of O. coerulescens span the coasts of Korea and its volcanic islands, Ulleungdo and Dokdo. Those on the volcanic islands exhibit morphological differences and have distinct hosts compared to those on the peninsula. The family of Orobanchaceae, encompassing both autotrophic and parasitic species, serves as a model for evolutionary studies of parasitic states. However, there are limited genome assemblies for the Orobanche genus. In our study, we produced approximately 100x ONT long reads to construct a chromosome-level genome of O. coerulescens. The resulting assembly has a total size of 3,648 Mb with an N50 value of 195 Mb, and 82.0% of BUSCO genes were identified as complete. Results of the repeat annotation revealed that 86.3% of the genome consisted of repeat elements, and 29,395 protein-coding genes were annotated. This chromosome-level genome will be an important biological resource for conserving biodiversity and further understanding parasitic plants.

ABSTRACT Precise identification of species is fundamental in microbial genomics and is crucial for understanding the microbial communities. While the 16S rRNA gene, particularly its V3-V4 regions, has been extensively employed for microbial identification, however has limitations in achieving species-level resolution. Advancements in long-read sequencing technologies have highlighted the rRNA operon as a more accurate marker for microbial classification and analysis than the 16S rRNA gene. This study aims to compare the accuracy of species classification and microbial community analysis using the rRNA operon versus the 16S rRNA gene. We evaluated the species classification accuracy of the rRNA operon,16S rRNA gene, and 16S rRNA V3-V4 regions using a BLAST-based method and a k-mer matching-based method with public data available from NCBI. We further performed simulations to model microbial community analysis. We accessed the performance using each marker in community composition estimation and differential abundance analysis. Our findings demonstrate that the rRNA operon offers an advantage over the 16S rRNA gene and its V3-V4 regions for species-level classification within the genus. When applied to microbial community analysis, the rRNA operon enables a more accurate determination of composition. Using the rRNA operon yielded more reliable results in differential abundance analysis as well.IMPORTANCEWe quantitatively demonstrated that the rRNA operon outperformed the 16S rRNA and its V3-V4 regions in accuracy for both individual species identification and species-level microbial community analysis. Our findings can provide guidelines for selecting appropriate markers in the field of microbial research.

IntroductionProbiotics, such as Lacticaseibacillus paracasei EG005, are gaining attention for their health benefits, particularly in reducing oxidative stress. The goal of this study was to reinforce the antioxidant capacity of EG005, along with comprehensive genomic analysis, with a focus on assessing superoxide dismutase (SOD) activity, acid resistance and bile tolerance, and safety.MethodsEG005 was screened for SOD activity and change of SOD activity was tested under various pH conditions. Its survival rates were assessed in acidic (pH 2.5) and bile salt (0.3%) conditions and the antibiotic MIC test and hemolysis test were performed to evaluate safety. Genetic analyses including functional identification and phylogenetic tree construction were performed. The SOD overexpression system was constructed using Ptuf, Pldh1, Plhd2, and Pldh3 strong promoters.ResultsEG005 demonstrated higher SOD activity compared to Lacticaseibacillus rhamnosus GG, with optimal activity at pH 7.0. It showed significant acid and bile tolerance, with survival rates recovering to 100% after 3 h in acidic conditions. Phylogenetic analysis confirmed that EG005 is closely related to other L. paracasei strains with ANI values above 98%. Overexpression of SOD using the Ptuf promoter resulted in a two-fold increase in activity compared to the controls. Additionally, EG005 exhibited no hemolytic activity and showed antibiotic susceptibility within safe limits.DiscussionOur findings highlight EG005’s potential as a probiotic with robust antioxidant activity and high tolerance to gastrointestinal conditions. Its unique genetic profile and enhanced SOD activity through strong promoter support its application in probiotic therapies and functional foods. Further research should be investigated to find the in vivo effects of EG005 on gut health and oxidative stress reduction. In addition, attB and attP-based recombination, combined with CRISPR-Cas9 technologies, could offer a more stable alternative for long-term sodA gene expression in commercial and medical applications.

Abstract Background False duplications in genome assemblies lead to false biological conclusions. We quantified false duplications in popularly used previous genome assemblies for platypus, zebra finch, and Anna’s Hummingbird, and their new counterparts of the same species generated by the Vertebrate Genomes Project, of which the Vertebrate Genomes Project pipeline attempted to eliminate false duplications through haplotype phasing and purging. These assemblies are among the first generated by the Vertebrate Genomes Project where there was a prior chromosomal level reference assembly to compare with. Results Whole genome alignments revealed that 4 to 16% of the sequences are falsely duplicated in the previous assemblies, impacting hundreds to thousands of genes. These lead to overestimated gene family expansions. The main source of the false duplications is heterotype duplications, where the haplotype sequences were relatively more divergent than other parts of the genome leading the assembly algorithms to classify them as separate genes or genomic regions. A minor source is sequencing errors. Ancient ATP nucleotide binding gene families have a higher prevalence of false duplications compared to other gene families. Although present in a smaller proportion, we observe false duplications remaining in the Vertebrate Genomes Project assemblies that can be identified and purged. Conclusions This study highlights the need for more advanced assembly methods that better separate haplotypes and sequence errors, and the need for cautious analyses on gene gains.

Abstract Background Many short-read genome assemblies have been found to be incomplete and contain mis-assemblies. The Vertebrate Genomes Project has been producing new reference genome assemblies with an emphasis on being as complete and error-free as possible, which requires utilizing long reads, long-range scaffolding data, new assembly algorithms, and manual curation. A more thorough evaluation of the recent references relative to prior assemblies can provide a detailed overview of the types and magnitude of improvements. Results Here we evaluate new vertebrate genome references relative to the previous assemblies for the same species and, in two cases, the same individuals, including a mammal (platypus), two birds (zebra finch, Anna’s hummingbird), and a fish (climbing perch). We find that up to 11% of genomic sequence is entirely missing in the previous assemblies. In the Vertebrate Genomes Project zebra finch assembly, we identify eight new GC- and repeat-rich micro-chromosomes with high gene density. The impact of missing sequences is biased towards GC-rich 5′-proximal promoters and 5′ exon regions of protein-coding genes and long non-coding RNAs. Between 26 and 60% of genes include structural or sequence errors that could lead to misunderstanding of their function when using the previous genome assemblies. Conclusions Our findings reveal novel regulatory landscapes and protein coding sequences that have been greatly underestimated in previous assemblies and are now present in the Vertebrate Genomes Project reference genomes.

Abstract Over the past decades, accumulating evidences have highlighted the gut microbiota as a key player in the brain functioning via microbiota–gut–brain axis, and accordingly, the beneficial role of several probiotic strains in cognitive ability also have been actively investigated. However, the majority of the research have demonstrated the effects against age-related cognitive decline or neurological disease. To this end, we aimed to investigate lactic acid bacteria strains having beneficial effects on the cognitive function of healthy young mice and elucidate underlying characteristics by carrying out nanopore sequencing-based genomics and metagenomics analysis. 8-week consumption of Streptococcus thermophilus EG007 demonstrated marked enhancements in behavior tests assessing short-term spatial and non-spatial learning and memory. It was revealed that EG007 possessed genes encoding various metabolites beneficial for a health condition in many aspects, including gamma-aminobutyric acid producing system, a neurotransmitter associated with mood and stress response. Also, by utilizing 16S–23S rRNA operon as a taxonomic marker, we identified more accurate species-level compositional changes in gut microbiota, which was increase of certain species, previously reported to have associations with mental health or down-regulation of inflammation or infection-related species. Moreover, correlation analysis revealed that the EG007-mediated altered microbiota had a significant correlation with the memory traits.

Abstract Background Africa is an important watershed in the genetic history of domestic cattle, as two lineages of modern cattle, Bos taurus and B. indicus, form distinct admixed cattle populations. Despite the predominant B. indicus nuclear ancestry of African admixed cattle, B. indicus mitochondria have not been found on the continent. This discrepancy between the mitochondrial and nuclear genomes has been previously hypothesized to be driven by male-biased introgression of Asian B. indicus into ancestral African B. taurus. Given that this hypothesis mandates extreme demographic assumptions relying on random genetic drift, we propose a novel hypothesis of selection induced by mitonuclear incompatibility and assess these hypotheses with regard to the current genomic status of African admixed cattle. Results By analyzing 494 mitochondrial and 235 nuclear genome sequences, we first confirmed the genotype discrepancy between mitochondrial and nuclear genome in African admixed cattle: the absence of B. indicus mitochondria and the predominant B. indicus autosomal ancestry. We applied approximate Bayesian computation (ABC) to assess the posterior probabilities of two selection hypotheses given this observation. The results of ABC indicated that the model assuming both male-biased B. indicus introgression and selection induced by mitonuclear incompatibility explains the current genomic discrepancy most accurately. Subsequently, we identified selection signatures at autosomal loci interacting with mitochondria that are responsible for integrity of the cellular respiration system. By contrast with B. indicus-enriched genome ancestry of African admixed cattle, local ancestries at these selection signatures were enriched with B. taurus alleles, concurring with the key expectation of selection induced by mitonuclear incompatibility. Conclusions Our findings support the current genome status of African admixed cattle as a potential outcome of male-biased B. indicus introgression, where mitonuclear incompatibility exerted selection pressure against B. indicus mitochondria. This study provides a novel perspective on African cattle demography and supports the role of mitonuclear incompatibility in the hybridization of mammalian species.

Unlike other Cirsium in Korea, Cirsium nipponicum (Island thistle) is distributed only on Ulleung Island, a volcanic island off the east coast of the Korean Peninsula, and a unique thistle with none or very small thorns. Although many researchers have questioned the origin and evolution of C. nipponicum, there is not much genomic information to estimate it. We thus assembled the complete chloroplast of C. nipponicum and reconstructed the phylogenetic relationships within the genus Cirsium. The chloroplast genome was 152,586 bp, encoding 133 genes consisting of 8 rRNA genes, 37 tRNA genes, and 88 protein-coding genes. We found 833 polymorphic sites and eight highly variable regions in chloroplast genomes of six Cirsium species by calculating nucleotide diversity, as well as 18 specific variable regions distinguished C. nipponicum from other Cirsium. As a result of phylogenetic analysis, C. nipponicum was closer to C. arvense and C. vulgare than native Cirsium in Korea: C. rhinoceros and C. japonicum. These results indicate that C. nipponicum is likely introduced through the north Eurasian root, not the mainland, and evolved independently in Ulleung Island. This study contributes to further understanding the evolutionary process and the biodiversity conservation of C. nipponicum on Ulleung Island.

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.

ABSTRACT Evidence for the concept of the “gut-brain axis” (GBA) has risen. Many types of research demonstrated the mechanism of the GBA and the effect of probiotic intake. Although many studies have been reported, most were focused on neurodegenerative disease and, it is still not clear what type of bacterial strains have positive effects. We designed an experiment to discover a strain that positively affects brain function, which can be recognized through changes in cognitive processes using healthy mice. The experimental group consisted of a control group and three probiotic consumption groups, namely, Lactobacillus acidophilus, Lacticaseibacillus paracasei, and Lacticaseibacillus rhamnosus. Three experimental groups fed probiotics showed an improved cognitive ability by cognitive-behavioral tests, and the group fed on L. acidophilus showed the highest score. To provide an understanding of the altered microbial composition effect on the brain, we performed full 16S-23S rRNA sequencing using Nanopore, and operational taxonomic units (OTUs) were identified at species level. In the group fed on L. acidophilus, the intestinal bacterial ratio of Firmicutes and Proteobacteria phyla increased, and the bacterial proportions of 16 species were significantly different from those of the control group. We estimated that the positive results on the cognitive behavioral tests were due to the increased proportion of the L. acidophilus EG004 strain in the subjects’ intestines since the strain can produce butyrate and therefore modulate neurotransmitters and neurotrophic factors. We expect that this strain expands the industrial field of L. acidophilus and helps understand the mechanism of the gut-brain axis. IMPORTANCE Recently, the concept of the “gut-brain axis” has risen and suggested that microbes in the GI tract affect the brain by modulating signal molecules. Although many pieces of research were reported in a short period, a signaling mechanism and the effects of a specific bacterial strain are still unclear. Besides, since most of the research was focused on neurodegenerative disease, the study with a healthy animal model is still insufficient. In this study, we show using a healthy animal model that a bacterial strain (Lactobacillus acidophilus EG004) has a positive effect on mouse cognitive ability. We experimentally verified an improved cognitive ability by cognitive behavioral tests. We performed full 16S-23S rRNA sequencing using a Nanopore MinION instrument and provided the gut microbiome composition at the species level. This microbiome composition consisted of candidate microbial groups as a biomarker that shows positive effects on cognitive ability. Therefore, our study suggests a new perspective for probiotic strain use applicable for various industrialization processes.

Lactobacillus acidophilus (L. acidophilus) is a representative probiotic and is widely used in many industrial products for its beneficial effects on human and animal health. This bacterium is exposed to harsh environments such as high temperatures for manufacturing industrial products, but cell yield under high temperatures is relatively low. To resolve this issue, we developed a new L. acidophilus strain with improved heat resistance while retaining the existing beneficial properties through the adaptive laboratory evolution (ALE) method. The newly developed strain, L. acidophilus EG008, has improved the existing limit of thermal resistance from 65°C to 75°C. Furthermore, we performed whole-genome sequencing and comparative genome analysis of wild-type and EG008 strains to unravel the molecular mechanism of improved heat resistance. Interestingly, only two single-nucleotide polymorphisms (SNPs) were different compared to the L. acidophilus wild-type. We identified that one of these SNPs is a non-synonymous SNP capable of altering the structure of MurD protein through the 435th amino acid change from serine to threonine. We believe that these results will directly contribute to any industrial field where L. acidophilus is applied. In addition, these results make a step forward in understanding the molecular mechanisms of lactic acid bacteria evolution under extreme conditions.

Abstract Background Indigenous cattle in Africa have adapted to various local environments to acquire superior phenotypes that enhance their survival under harsh conditions. While many studies investigated the adaptation of overall African cattle, genetic characteristics of each breed have been poorly studied. Results We performed the comparative genome-wide analysis to assess evidence for subspeciation within species at the genetic level in trypanotolerant N’Dama cattle. We analysed genetic variation patterns in N’Dama from the genomes of 101 cattle breeds including 48 samples of five indigenous African cattle breeds and 53 samples of various commercial breeds. Analysis of SNP variances between cattle breeds using wMI, XP-CLR, and XP-EHH detected genes containing N’Dama-specific genetic variants and their potential associations. Functional annotation analysis revealed that these genes are associated with ossification, neurological and immune system. Particularly, the genes involved in bone formation indicate that local adaptation of N’Dama may engage in skeletal growth as well as immune systems. Conclusions Our results imply that N’Dama might have acquired distinct genotypes associated with growth and regulation of regional diseases including trypanosomiasis. Moreover, this study offers significant insights into identifying genetic signatures for natural and artificial selection of diverse African cattle breeds.

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.

Although the chemical, physical, and nutritional properties of bovine milk have been extensively studied, only a few studies have attempted to characterize milk-synthesizing genes using RNA-seq data. RNA-seq data was collected from 21 Holstein samples, along with group information about milk production ability; milk yield; and protein, fat, and solid contents. Meta-analysis was employed in order to generally characterize genes related to milk production. In addition, we attempted to investigate the relationship between milk related traits, parity, and lactation period. We observed that milk fat is highly correlated with lactation period; this result indicates that this effect should be considered in the model in order to accurately detect milk production related genes. By employing our developed model, 271 genes were significantly (false discovery rate [FDR] adjusted p-value

Selective sweep can cause genetic differentiation across populations, which allows for the identification of possible causative regions/genes underlying important traits. The pig has experienced a long history of allele frequency changes through artificial selection in the domestication process. We obtained an average of 329,482,871 sequence reads for 24 pigs from three pig breeds: Yorkshire (n = 5), Landrace (n = 13), and Duroc (n = 6). An average read depth of 11.7 was obtained using whole-genome resequencing on an Illumina HiSeq2000 platform. In this study, cross-population extended haplotype homozygosity and cross-population composite likelihood ratio tests were implemented to detect genes experiencing positive selection for the genome-wide resequencing data generated from three commercial pig breeds. In our results, 26, 7, and 14 genes from Yorkshire, Landrace, and Duroc, respectively were detected by two kinds of statistical tests. Significant evidence for positive selection was identified on genes ST6GALNAC2 and EPHX1 in Yorkshire, PARK2 in Landrace, and BMP6, SLA-DQA1, and PRKG1 in Duroc.These genes are reportedly relevant to lactation, reproduction, meat quality, and growth traits. To understand how these single nucleotide polymorphisms (SNPs) related positive selection affect protein function, we analyzed the effect of non-synonymous SNPs. Three SNPs (rs324509622, rs80931851, and rs80937718) in the SLA-DQA1 gene were significant in the enrichment tests, indicating strong evidence for positive selection in Duroc. Our analyses identified genes under positive selection for lactation, reproduction, and meat-quality and growth traits in Yorkshire, Landrace, and Duroc, respectively.

African indigenous cattle breeds have been reared within the continent for millennia. Due to the harsh tropical environmental conditions (e.g. sever disease and parasite prevalence, high temperature, feed and water scarcity) where they evolved, they have developed various levels of tropical environment adaptation attributes. In order to explore the genomic signatures of tropical environment adaptation in African cattle, we compared the whole genomes of East African Indicus cattle breeds with European and Asian Taurine cattle breeds using XP-EHH and XP-CLR population statistical methods. Several genes involved in various biological processes and pathways related to domestication and behavior (dopaminergic and glutamatergic synapse), feeding and metabolism (gastric acid secretion, metabolic pathways), thermotolerance (sphingolipid and Wnt signaling), immune system response (T cell receptor signaling), and growth and reproduction (osteoblast differentiation, fibroblast migration) were identified from our genome analysis. Genes associated with tick and parasite resistance traits such as keratin genes, collagen genes, calcium signaling, and tumor necrosis factor proteins were also identified. The genes and pathways identified in this study improve our understanding of the biological mechanisms of tropical environment adaptation of African Indicus cattle breeds, which may allow us to use them for genomic selection programs. This result presents a basis for further study and may help to develop vaccines for tick and gastrointestinal parasite challenge.

The missing heritability has been a major problem in the analysis of best linear unbiased prediction (BLUP). We introduced the traditional genome-wide association study (GWAS) into the BLUP to improve the heritability estimation. We analyzed eight pork quality traits of the Berkshire breeds using GWAS and BLUP. GWAS detects the putative quantitative trait loci regions given traits. The single nucleotide polymorphisms (SNPs) were obtained using GWAS results with p value

At least 150 indigenous African cattle breeds have been named, but the majority of African cattle populations remain largely uncharacterized. As cattle breeds and populations in Africa adapted to various local environmental conditions, they acquired unique features. We know now that the history of African cattle was particularly complex and while several of its episodes remain debated, there is no doubt that African cattle population evolved dramatically over time. Today, we find a mosaic of genetically diverse population from the purest Bos taurus to the nearly pure Bos indicus. African cattle are now found all across the continent, with the exception of the Sahara and the river Congo basin. They are found on the rift valley highlands as well as below sea level in the Afar depression. These unique livestock genetic resources are in danger to disappear rapidly following uncontrolled crossbreeding and breed replacements with exotic breeds. Breeding improvement programs of African indigenous livestock remain too few while paradoxically the demand of livestock products is continually increasing. Many African indigenous breeds are endangered now, and their unique adaptive traits may be lost forever. This paper reviews the unique known characteristics of indigenous African cattle populations while describing the opportunities, the necessity and urgency to understand and utilize these resources to respond to the needs of the people of the continent and to the benefit of African farmers.

Copy number variation (CNV) or single nucleotide phlyorphism (SNP) is useful genetic resource to aid in understanding complex phenotypes or deseases susceptibility. Although thousands of CNVs and SNPs are currently avaliable in the public databases, they are somewhat difficult to use for analyses without visualization tools. We developed a web-based tool called the VCS (visualization of CNV or SNP) to visualize the CNV or SNP detected. The VCS tool can assist to easily interpret a biological meaning from the numerical value of CNV and SNP. The VCS provides six visualization tools: i) the enrichment of genome contents in CNV; ii) the physical distribution of CNV or SNP on chromosomes; iii) the distribution of log2 ratio of CNVs with criteria of interested; iv) the number of CNV or SNP per binning unit; v) the distribution of homozygosity of SNP genotype; and vi) cytomap of genes within CNV or SNP region.

Best linear unbiased prediction (BLUP) has been used to estimate the fixed effects and random effects of complex traits. Traditionally, genomic relationship matrix-based (GRM) and random marker-based BLUP analyses are prevalent to estimate the genetic values of complex traits. We used three methods: GRM-based prediction (G-BLUP), random marker-based prediction using an identity matrix (so-called single-nucleotide polymorphism [SNP]-BLUP), and SNP-SNP variance-covariance matrix (so-called SNP-GBLUP). We used 35,675 SNPs and R package "rrBLUP" for the BLUP analysis. The SNP-SNP relationship matrix was calculated using the GRM and Sherman-Morrison-Woodbury lemma. The SNP-GBLUP result was very similar to G-BLUP in the prediction of genetic values. However, there were many discrepancies between SNP-BLUP and the other two BLUPs. SNP-GBLUP has the merit to be able to predict genetic values through SNP effects.

The aim of this study was to investigate which of the gut microbes respond to probiotic intervention, as well as study whether they are associated with gastrointestinal symptoms in a healthy adult human. For the experimental purpose, twenty-one healthy adults were recruited and received probiotic mixture, which is composed of five Lactobacilli strains and two Bifidobacteria strains, once a day for 60 days. Defecation survey and Bioelectrical Impedance Analysis were conducted pre- and post-administration to measure phenotypic differences. Stool samples of the subjects were collected twice.The statistical analysis was performed for pair designed metagenome data with 11 phenotypic records of the bioelectrical impedance body composition analyzer and 6 responses of the questionnaires about gastrointestinal symptom. Furthemore, correlation-based network analysis was conducted for exploring complex relationships among microbiome communities.The abundances of Citrobacter, Klebsiella, and Methanobrevibacter were significantly reduced, which are strong candidates to be highly affected by the probiotic administration. In addition, interaction effects were observed between flatulence symptom attenuation and decreasing patterns of the Methanobrevibacter abundance.These results reveal that probiotic intervention modulated the composition of gut microbiota and reduced the abundance of potential pathogens (i.e. Citrobacter and Klebsiella). In addition, methanogens (i.e. Methanobrevibacter) associated with the gastrointestinal symptom in an adult human.

The whole-genome sequence of watermelon (Citrullus lanatus (Thunb.) Matsum. & Nakai), a valuable horticultural crop worldwide, was released in 2013. Here, we compared a de novo-based approach (DBA) to a reference-based approach (RBA) using RNA-seq data, to aid in efforts to improve the annotation of the watermelon reference genome and to obtain biological insight into male-sterility in watermelon. We applied these techniques to available data from two watermelon lines: the male-sterile line DAH3615-MS and the male-fertile line DAH3615. Using DBA, we newly annotated 855 watermelon transcripts, and found gene functional clusters predicted to be related to stimulus responses, nucleic acid binding, transmembrane transport, homeostasis, and Golgi/vesicles. Among the DBA-annotated transcripts, 138 de novo-exclusive differentially-expressed genes (DEDEGs) related to male sterility were detected. Out of 33 randomly selected newly annotated transcripts and DEDEGs, 32 were validated by RT-qPCR. This study demonstrates the usefulness and reliability of the de novo transcriptome assembly in watermelon, and provides new insights for researchers exploring transcriptional blueprints with regard to the male sterility.

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.

Indigenous (native) breeds of livestock have higher disease resistance and adaptation to the environment due to high genetic diversity. Even though their extinction rate is accelerated due to the increase of commercial breeds, natural disaster, and civil war, there is a lack of well-established databases for the native breeds. Thus, we constructed the native pig and chicken breed database (NPCDB) which integrates available information on the breeds from around the world. It is a nonprofit public database aimed to provide information on the genetic resources of indigenous pig and chicken breeds for their conservation. The NPCDB (http://npcdb.snu.ac.kr/) provides the phenotypic information and population size of each breed as well as its specific habitat. In addition, it provides information on the distribution of genetic resources across the country. The database will contribute to understanding of the breed’s characteristics such as disease resistance and adaptation to environmental changes as well as the conservation of indigenous genetic resources.

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.

Copy number variations (CNVs), important genetic factors for study of human diseases, may have as large of an effect on phenotype as do single nucleotide polymorphisms. Indeed, it is widely accepted that CNVs are associated with differential disease susceptibility. However, the relationships between CNVs and gene expression have not been characterized in the horse. In this study, we investigated the effects of copy number deletion in the blood and muscle transcriptomes of Thoroughbred racing horses. We identified a total of 1,246 CNVs of deletion polymorphisms using DNA re-sequencing data from 18 Thoroughbred racing horses. To discover the tendencies between CNV status and gene expression levels, we extracted CNVs of four Thoroughbred racing horses of which RNA sequencing was available. We found that 252 pairs of CNVs and genes were associated in the four horse samples. We did not observe a clear and consistent relationship between the deletion status of CNVs and gene expression levels before and after exercise in blood and muscle. However, we found some pairs of CNVs and associated genes that indicated relationships with gene expression levels: a positive relationship with genes responsible for membrane structure or cytoskeleton and a negative relationship with genes involved in disease. This study will lead to conceptual advances in understanding the relationship between CNVs and global gene expression in the horse.

Pork from Jeju black pig (population J) and Berkshire (population B) has a unique market share in Korea because of their high meat quality. Due to the high demand of this pork, traceability of the pork to its origin is becoming an important part of the consumer demand. To examine the feasibility of such a system, we aim to provide basic genetic information of the two black pig populations and assess the possibility of genetically distinguishing between the two breeds. Muscle samples were collected from slaughter houses in Jeju Island and Namwon, Chonbuk province, Korea, for populations J and B, respectively. In total 800 Jeju black pigs and 351 Berkshires were genotyped at thirteen microsatellite (MS) markers. Analyses on the genetic diversity of the two populations were carried out in the programs MS toolkit and FSTAT. The population structure of the two breeds was determined by a Bayesian clustering method implemented in structure and by a phylogenetic analysis in Phylip. Population J exhibited higher mean number of alleles, expected heterozygosity and observed heterozygosity value, and polymorphism information content, compared to population B. The FIS values of population J and population B were 0.03 and −0.005, respectively, indicating that little or no inbreeding has occurred. In addition, genetic structure analysis revealed the possibility of gene flow from population B to population J. The expected probability of identify value of the 13 MS markers was 9.87×10−14 in population J, 3.17×10−9 in population B, and 1.03×10−12 in the two populations. The results of this study are useful in distinguishing between the two black pig breeds and can be used as a foundation for further development of DNA markers.

One of the most important traits for both animal science and livestock production is the number of offspring for a species. This study was performed to identify differentially evolved genes and their distinct functions that influence the number of offspring at birth by comparative analysis of eight monotocous mammals and seven polytocous mammals in a number of scopes: specific amino acid substitution with site-wise adaptive evolution, gene expansion and specific orthologous group. The mutually exclusive amino acid substitution among the 16 mammalian species identified five candidate genes. These genes were both directly and indirectly related to ovulation. Furthermore, in monotocous mammals, the EPH gene family was found to have undergone expansion. Previously, the EPHA4 gene was found to positively affect litter size in pigs and supports the possibility of the EPH gene playing a role in determining the number of offspring per birth. The identified genes in this study offer a basis from which the differences between monotocous and polytocous species can be studied. Furthermore, these genes may harbor some clues to the underlying mechanism, which determines litter size and may prove useful for livestock breeding strategies.

The increasing importance of meat quality has implications for animal breeding programs. Research has revealed much about the genetic background of pigs, and many studies have revealed the importance of various genetic factors. Since meat quality is a complex trait which is affected by many factors, consideration of the overall phenotype is very useful to study meat quality. For integrating the phenotypes, we used principle component analysis (PCA). The significant SNPs refer to results of the GRAMMAR method against PC1, PC2 and PC3 of 14 meat quality traits of 181 Duroc pigs. The Genome-wide association study (GWAS) found 26 potential SNPs affecting various meat quality traits. The loci identified are located in or near 23 genes. The SNPs associated with meat quality are in or near five genes (ANK1, BMP6, SHH, PIP4K2A, and FOXN2) and have been reported previously. Twenty-five of the significant SNPs also located in meat quality-related QTL regions, these result supported the QTL effect indirectly. Each single gene typically affects multiple traits. Therefore, it is a useful approach to use integrated traits for the various traits at the same time. This innovative approach using integrated traits could be applied on other GWAS of complex-traits including meat-quality, and the results will contribute to improving meat-quality of pork.

To increase plumage color uniformity and understand the genetic background of Korean chickens, we performed a genome-wide association study of different plumage color in Korean native chickens. We analyzed 60K SNP chips on 279 chickens with GEMMA methods for GWAS and estimated the genetic heritability for plumage color. The estimated heritability suggests that plumage coloration is a polygenic trait. We found new loci associated with feather pigmentation at the genome-wide level and from the results infer that there are additional genetic effect for plumage color. The results will be used for selecting and breeding chicken for plumage color uniformity.

In this study, protein domains with cellulase activity in goat rumen microbes were investigated using metagenomic and bioinformatic analyses. After the complete genome of goat rumen microbes was obtained using a shotgun sequencing method, 217,892,109 pair reads were filtered, including only those with 70% identity, 100-bp matches, and thresholds below E−10 using METAIDBA. These filtered contigs were assembled and annotated using blastN against the NCBI nucleotide database. As a result, a microbial community structure with 1431 species was analyzed, among which Prevotella ruminicola 23 bacteria and Butyrivibrio proteoclasticus B316 were the dominant groups. In parallel, 201 sequences related with cellulase activities (EC.3.2.1.4) were obtained through blast searches using the enzyme.dat file provided by the NCBI database. After translating the nucleotide sequence into a protein sequence using Interproscan, 28 protein domains with cellulase activity were identified using the HMMER package with threshold E values below 10−5. Cellulase activity protein domain profiling showed that the major protein domains such as lipase GDSL, cellulase, and Glyco hydro 10 were present in bacterial species with strong cellulase activities. Furthermore, correlation plots clearly displayed the strong positive correlation between some protein domain groups, which was indicative of microbial adaption in the goat rumen based on feeding habits. This is the first metagenomic analysis of cellulase activity protein domains using bioinformatics from the goat rumen.

Although there have been many studies of native Koreancattle, Hanwoo, there have been no selective sweep studies inthese animals. This study was performed to characterizegenetic variation and identify selective signatures. We sequencedthe genomes of 12 cattle, and identified 15125420 SNPs,1768114 INDELs, and 3445 CNVs. The SNPs, INDELs, andCNVs were similarly distributed throughout the genome, andhighly variable regions were shown to contain the BoLA familyand GPR180, which are related to adaptive immunity. We alsoidentified the domestication footprints of the Hanwoopopulation by searching for selective sweep signatures, whichrevealed the RCN2 gene related to BPV resistance. The resultsof this study may contribute to genetic improvement of theHanwoo population in Korea. [BMB Reports 2013; 46(7):346-351]

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 phylogeny of classical swine fever virus (CSFV), the causative agent of classical swine fever (CSF), has been investigated extensively. However, no evolutionary research has been performed using the whole CSFV genome. In this study, we used 37 published genome sequences to investigate the time-calibrated phylogenomics of CSFV. In phylogenomic trees based on Bayesian inference (BI) and Maximum likelihood (ML), the 37 isolates were categorized into five genetic types (1.1, 1.2, 2.1, 2.3, and 3.4). Subgenotype 1.1 is divided into 3 groups and 1 unclassified isolate, 2.1 into 4 groups, 2.3 into 2 groups and 1 unclassified isolate, and subgenotype 1.2 and 3.4 consisted of one isolate each. We did not observe an apparent temporal or geographical relationship between isolates. Of the 14 genomic regions, NS4B showed the most powerful phylogenetic signal. Results of this evolutionary study using Bayesian coalescent approach indicate that CSFV has evolved at a rate of 13×.010-4 substitutions per site per year. The most recent common ancestor of CSFV appeared 2770.2 years ago, which was about 8000 years after pig domestication. The effective population size of CSFV underwent a slow increase until the 1950s, after which it has remained constant.

Japanese quail (Coturnix coturnix japonica) reach sexual maturity earlier, breed rapidly and successfully, and cost less and require less space than other birds raised commercially. Given the value of this species for food production and experimental use, more studies are necessary to determine chromosomal regions and genes associated with gender and breed-differentiation. This study employed Trinity and edgeR for transcriptome analysis of next-generation RNA-seq data, which included 4 tissues obtained from 3 different breeding lines of Japanese quail (random bred control, heavy weight, low weight). Differentially expressed genes shared between female and male tissue contrast groups were analyzed to identify genes related to sexual dimorphism as well as potential novel candidate genes for molecular sexing. Several of the genes identified in the present study as significant sex-related genes have been previously found in avian gene expression analyses (NIPBL, UBAP2), and other genes found differentially expressed in this study and not previously associated with sex-related differences may be considered potential candidates for molecular sexing (TERA, MYP0, PPR17, CASQ2). Additionally, other genes likely associated with neuronal and brain development (CHKA, NYAP), as well as body development and size differentiation (ANKRD26, GRP87) in quail were identified. Expression of homeobox protein regulating genes (HXC4, ISL1) shared between our two sex-related contrast groups (Female Brain vs. Male Brain and Ovary vs. Testis) indicates that these genes may regulate sex-specific anatomical development. Results reveal genetic features of the quail breed and could allow for more effective molecular sexing as well as selective breeding for traits important in commercial production.

Consumer attention to food safety has increased rapidly due to animal-related diseases; therefore, it is important to identify their places of origin (POO) for safety purposes. However, only a few studies have addressed this issue and focused on machine learning-based approaches. In the present study, classification analyses were performed using a customized SNP chip for POO prediction. To accomplish this, 4,122 pigs originating from 104 farms were genotyped using the SNP chip. Several factors were considered to establish the best prediction model based on these data. We also assessed the applicability of the suggested model using a kinship coefficient-filtering approach. Our results showed that the LogitBoost-based prediction model outperformed other classifiers in terms of classification performance under most conditions. Specifically, a greater level of accuracy was observed when a higher kinship-based cutoff was employed. These results demonstrated the applicability of a machine learning-based approach using SNP chip data for practical traceability.

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.

Previous studies of horse RNA-seq were performed by mapping sequence reads to the reference genome during transcriptome analysis. However in this study, we focused on two main ideas. First, differentially expressed genes (DEGs) were identified by de novo-based analysis (DBA) in RNA-seq data from six Thoroughbreds before and after exercise, here-after referred to as "de novo unique differentially expressed genes" (DUDEG). Second, by integrating both conventional DEGs and genes identified as being selected for during domestication of Thoroughbred and Jeju pony from whole genome re-sequencing (WGS) data, we give a new concept to the definition of DEG. We identified 1,034 and 567 DUDEGs in skeletal muscle and blood, respectively. DUDEGs in skeletal muscle were significantly related to exercise-induced stress biological process gene ontology (BP-GO) terms: 'immune system process'; 'response to stimulus'; and, 'death' and a KEGG pathways: 'JAK-STAT signaling pathway'; 'MAPK signaling pathway'; 'regulation of actin cytoskeleton'; and, 'p53 signaling pathway'. In addition, we found TIMELESS, EIF4A3 and ZNF592 in blood and CHMP4C and FOXO3 in skeletal muscle, to be in common between DUDEGs and selected genes identified by evolutionary statistics such as FST and Cross Population Extended Haplotype Homozygosity (XP-EHH). Moreover, in Thoroughbreds, three out of five genes (CHMP4C, EIF4A3 and FOXO3) related to exercise response showed relatively low nucleotide diversity compared to the Jeju pony. DUDEGs are not only conceptually new DEGs that cannot be attained from reference-based analysis (RBA) but also supports previous RBA results related to exercise in Thoroughbred. In summary, three exercise related genes which were selected for during domestication in the evolutionary history of Thoroughbred were identified as conceptually new DEGs in this study.

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

Recent studies in population of European ancestry have shown that 30% ~ 50% of heritability for human complex traits such as height and body mass index, and common diseases such as schizophrenia and rheumatoid arthritis, can be captured by common SNPs and that genetic variation attributed to chromosomes are in proportion to their length. Using genome-wide estimation and partitioning approaches, we analysed 49 human quantitative traits, many of which are relevant to human diseases, in 7,170 unrelated Korean individuals genotyped on 326,262 SNPs. For 43 of the 49 traits, we estimated a nominally significant (P

Obesity represents a major global public health problem that increases the risk for cardiovascular or metabolic disease. The pigs represent an exceptional biomedical model related to energy metabolism and obesity in humans. To pinpoint causal genetic factors for a common form of obesity, we conducted local genomic de novo sequencing, 18.2 Mb, of a porcine QTL region affecting fatness traits, and carried out SNP association studies for backfat thickness and intramuscular fat content in pigs. In order to relate the association studies in pigs to human obesity, we performed a targeted genome wide association study for subcutaneous fat thickness in a cohort population of 8,842 Korean individuals. These combined association studies in human and pig revealed a significant SNP located in a gene family with sequence similarity 73, member A (FAM73A) associated with subscapular skin-fold thickness in humans (rs4121165, GC-corrected p-value = 0.0000175) and with backfat thickness in pigs (ASGA0029495, p-value = 0.000031). Our combined association studies also suggest that eight neuronal genes are responsible for subcutaneous fat thickness: NEGR1, SLC44A5, PDE4B, LPHN2, ELTD1, ST6GALNAC3, ST6GALNAC5, and TTLL7. These results provide strong support for a major involvement of the CNS in the genetic predisposition to a common form of obesity.