Comparison of characteristics of long noncoding RNA in Hanwoo according to sex.

OBJECTIVE: Cattle were some of the first animals domesticated by humans for the production of milk, meat, etc. Long noncoding RNA (lncRNA) is defined as longer than 200 bp in non-protein coding transcripts. lncRNA is known to function in regulating gene expression and is currently being studied in a variety of livestock including cattle. The purpose of this study is to analyze the characteristics of lncRNA according to sex in Hanwoo cattle. METHODS: This study was conducted using the skeletal muscles of 9 Hanwoo cattle include bulls, steers and cows. RNA was extracted from skeletal muscle of Hanwoo. Sequencing was conducted using Illumina HiSeq2000 and mapped to the Bovine Taurus genome. The expression levels of lncRNAs were measured by DEGseq and quantitative trait loci (QTL) data base was used to identify QTLs associated with lncRNA. The python script was used to match the nearby genes. RESULTS: In this study, the expression patterns of transcripts of bulls, steers and cows were identified. And we identified significantly differentially expressed lncRNAs in bulls, steers and cows. In addition, characteristics of lncRNA which express differentially in muscles according to the sex of Hanwoo were identified. As a result, we found differentially expressed lncRNAs according to sex were related to shear force and body weight. CONCLUSION: This study was classified and characterized lncRNA which differentially expressed by sex in Hanwoo cattle. We believe that the characterization of lncRNA by sex of Hanwoo will be helpful for future studies of the physiological mechanisms of Hanwoo cattle.

Analysis of genetic characteristics of pig breeds using information on single nucleotide polymorphisms.

OBJECTIVE: This study was undertaken to investigate the genetic characteristics of Berkshire (BS), Landrace (LR), and Yorkshire (YS) pig breeds raised in the Great Grandparents pig farms using the single nucleotide polymorphisms (SNP) information. METHODS: A total of 25,921 common SNP genotype markers in three pig breeds were used to estimate the expected heterozygosity (HE), polymorphism information content, F-statistics (FST), linkage disequilibrium (LD) and effective population size (Ne). RESULTS: The chromosome-wise distribution of FST in BS, LR, and YS populations were within the range of 0-0.36, and the average FST value was estimated to be 0.07±0.06. This result indicated some level of genetic segregation. An average LD (r2) for the BS, LR, and YS breeds was estimated to be approximately 0.41. This study also found an average Ne of 19.9 (BS), 31.4 (LR), and 34.1 (YS) over the last 5th generations. The effective population size for the BS, LR, and YS breeds decreased at a consistent rate from 50th to 10th generations ago. With a relatively faster Ne decline rate in the past 10th generations, there exists possible evidence for intensive selection practices in pigs in the recent past. CONCLUSION: To develop customized chips for the genomic selection of various breeds, it is important to select and utilize SNP based on the genetic characteristics of each breed. Since the improvement efficiency of breed pigs increases sharply by the population size, it is important to increase test units for the improvement and it is desirable to establish the pig improvement network system to expand the unit of breed pig improvement through the genetic connection among breed pig farms.

In silico approaches to identify the functional and structural effects of non-synonymous SNPs in selective sweeps of the Berkshire pig genome.

OBJECTIVE: Non-synonymous single nucleotide polymorphisms (nsSNPs) were identified in Berkshire selective sweep regions and then were investigated to discover genetic nsSNP mechanisms that were potentially associated with Berkshire domestication and meat quality. We further used bioinformatics tools to predict damaging amino-acid substitutions in Berkshire-related nsSNPs. METHODS: nsSNPs were examined in whole genome resequencing data of 110 pigs, including 14 Berkshire pigs, generated using the Illumina Hiseq2000 platform to identify variations that might affect meat quality in Berkshire pigs. RESULTS: Total 65,550 nsSNPs were identified in the mapped regions; among these, 319 were found in Berkshire selective-sweep regions reported in a previous study. Genes encompassing these nsSNPs were involved in lipid metabolism, intramuscular fatty-acid deposition, and muscle development. The effects of amino acid change by nsSNPs on protein functions were predicted using sorting intolerant from tolerant and polymorphism phenotyping V2 to reveal their potential roles in biological processes that may correlate with the unique Berkshire meat-quality traits. CONCLUSION: Our nsSNP findings confirmed the history of Berkshire pigs and illustrated the effects of domestication on generic-variation patterns. Our novel findings, which are generally consistent with those of previous studies, facilitated a better understanding of Berkshire domestication. In summary, we extensively investigated the relationship between genomic composition and phenotypic traits by scanning for nsSNPs in large-scale whole-genome sequencing data.

Characterization of gene expression and genetic variation of horse ERBB receptor feedback inhibitor 1 in Thoroughbreds.

OBJECTIVE: This study aimed to test the expression patterns of ERBB receptor feedback inhibitor 1 (ERRFI1) before and after exercise and the association of non-synonymous single-nucleotide polymorphisms (nsSNPs) of horse ERRFI1 with racing traits in Thoroughbreds. METHODS: We performed bioinformatics and gene expression analyses for horse ERRFI1. Transcription factor (TF) binding sites in the 5'-regulatory region of this gene were identified through a tool for prediction of TF-binding site (PROMO). A general linear model was used to detect the association between the nsSNP (LOC42830758 A to G) and race performance. RESULTS: Quantitative polymerase chain reaction analysis showed that expression level of ERRFI1 after exercise was 1.6 times higher than that before exercise. Ten transcription factors were predicted from the ERRFI1 regulatory region. A novel nsSNP (LOC42830758 A to G) was found in ERRFI1, which was associated with three racing traits including average prize money, average racing index, and 3-year-old starts percentile ranking. CONCLUSION: Our analysis will be helpful as a basis for studying genes and SNPs that affect race performance in racehorses.

Validation of selection accuracy for the total number of piglets born in Landrace pigs using genomic selection.

OBJECTIVE: This study was to determine the relationship between estimated breeding value and phenotype information after farrowing when juvenile selection was made in candidate pigs without phenotype information. METHODS: After collecting phenotypic and genomic information for the total number of piglets born by Landrace pigs, selection accuracy between genomic breeding value estimates using genomic information and breeding value estimates of best linear unbiased prediction (BLUP) using conventional pedigree information were compared. RESULTS: Genetic standard deviation (σa ) for the total number of piglets born was 0.91. Since the total number of piglets born for candidate pigs was unknown, the accuracy of the breeding value estimated from pedigree information was 0.080. When genomic information was used, the accuracy of the breeding value was 0.216. Assuming that the replacement rate of sows per year is 100% and generation interval is 1 year, genetic gain per year is 0.346 head when genomic information is used. It is 0.128 when BLUP is used. CONCLUSION: Genetic gain estimated from single step best linear unbiased prediction (ssBLUP) method is by 2.7 times higher than that the one estimated from BLUP method, i.e., 270% more improvement in efficiency.

Exploring evidence of positive selection reveals genetic basis of meat quality traits in Berkshire pigs through whole genome sequencing.

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.

Single nucleotide polymorphism analysis of Korean native chickens using next generation sequencing data.

There are five native chicken lines in Korea, which are mainly classified by plumage colors (black, white, red, yellow, gray). These five lines are very important genetic resources in the Korean poultry industry. Based on a next generation sequencing technology, whole genome sequence and reference assemblies were performed using Gallus_gallus_4.0 (NCBI) with whole genome sequences from these lines to identify common and novel single nucleotide polymorphisms (SNPs). We obtained 36,660,731,136 ± 1,257,159,120 bp of raw sequence and average 26.6-fold of 25-29 billion reference assembly sequences representing 97.288 % coverage. Also, 4,006,068 ± 97,534 SNPs were observed from 29 autosomes and the Z chromosome and, of these, 752,309 SNPs are the common SNPs across lines. Among the identified SNPs, the number of novel- and known-location assigned SNPs was 1,047,951 ± 14,956 and 2,948,648 ± 81,414, respectively. The number of unassigned known SNPs was 1,181 ± 150 and unassigned novel SNPs was 8,238 ± 1,019. Synonymous SNPs, non-synonymous SNPs, and SNPs having character changes were 26,266 ± 1,456, 11,467 ± 604, 8,180 ± 458, respectively. Overall, 443,048 ± 26,389 SNPs in each bird were identified by comparing with dbSNP in NCBI. The presently obtained genome sequence and SNP information in Korean native chickens have wide applications for further genome studies such as genetic diversity studies to detect causative mutations for economic and disease related traits.

Genome Wide Association Studies Using Multiple-lactation Breeding Value in Holsteins.

A genome wide association study was conducted using estimated breeding value (EBV) for milk production traits from 1st to 4th lactation. Significant single nucleotide polymorphism (SNP) markers were selected for each trait and the differences were compared by lactation. DNA samples were taken from 456 animals with EBV which are Holstein proven bulls whose semen is being sold or the daughters of old proven bulls whose semen is no longer being sold in Korea. High density genome wide SNP genotype was investigated and the significance of markers associated with traits was tested using the breeding value estimated by a multiple lactation model as a dependent variant. As the result of significance comparisons by lactations, several differences were found between the first lactation and subsequent lactations (from second to 4th lactation). A similar trend was noted in mean deviation and correlation of the estimated effects by lactation. Since there was a difference in the genes associated with EBV for each trait between first and subsequent lactations, a multi-lactation model in which lactation is considered as a different trait is genetically useful. Also, significant markers in all lactations and common markers for different traits were detected, which can be used as markers for quantitative trait loci exploration and marker assisted selection in milk production traits.

Differential regulation of microRNA transcriptome in chicken lines resistant and susceptible to necrotic enteritis disease.

Necrotic enteritis (NE) is a re-emerging disease as a result of increased restriction on the use of antibiotics in poultry. However, the molecular mechanisms underlying the pathogenesis of NE are unclear. Small RNA transcriptome analysis was performed using spleen and intestinal intraepithelial lymphocytes (IEL) from 2 inbred chicken lines selected for resistance or susceptibility to Marek's disease (MD) in an experimentally induced model of avian NE to investigate whether microRNA (miRNA) control the expression of genes associated with host response to pathogen challenge. Unique miRNA represented only 0.02 to 0.04% of the total number of sequences obtained, of which 544 were unambiguously identified. Hierarchical clustering revealed that most of miRNA in IEL were highly expressed in the MD-susceptible line 7.2 compared with MD-resistant line 6.3. Reduced CXCL14 gene expression was correlated with differential expression of several unique miRNA in MD-resistant chickens, whereas TGFßR2 gene expression was correlated with altered gga-miR-216 miRNA levels in MD-susceptible animals. In conclusion, miRNA profiling and deep sequencing of small RNA in experimental models of infectious diseases may be useful for further understanding of host-pathogen interactions, and for providing insights into genetic markers of disease resistance.

Copy number deletion has little impact on gene expression levels in racehorses.

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.

Genetic traceability of black pig meats using microsatellite markers.

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.

A genomic estimated breeding value-assisted reduction method of single nucleotide polymorphism sets: a novel approach for determining the cutoff thresholds in genome-wide association studies and best linear unbiased prediction.

Traditionally, the p-value is the criterion for the cutoff threshold to determine significant markers in genome-wide association studies (GWASs). Choosing the best subset of markers for the best linear unbiased prediction (BLUP) for improved prediction ability (PA) has become an interesting issue. However, when dealing with many traits having the same marker information, the p-values' themselves cannot be used as an obvious solution for having a confidence in GWAS and BLUP. We thus suggest a genomic estimated breeding value-assisted reduction method of the single nucleotide polymorphism (SNP) set (GARS) to address these difficulties. GARS is a BLUP-based SNP set decision presentation. The samples were Landrace pigs and the traits used were back fat thickness (BF) and daily weight gain (DWG). The prediction abilities (PAs) for BF and DWG for the entire SNP set were 0.8 and 0.8, respectively. By using the correlation between genomic estimated breeding values (GEBVs) and phenotypic values, selecting the cutoff threshold in GWAS and the best SNP subsets in BLUP was plausible as defined by GARS method. 6,000 SNPs in BF and 4,000 SNPs in DWG were considered as adequate thresholds. Gene Ontology (GO) analysis using the GARS results of the BF indicated neuron projection development as the notable GO term, whereas for the DWG, the main GO terms were nervous system development and cell adhesion.

Genome-wide association studies on collagen contents trait for meat quality in Hanwoo.

Beef consumers valued meat quality traits such as texture, tenderness, juiciness, flavor, and meat color that determining consumers' purchasing decision. Most research on meat quality has focused on marbling, a key characteristic related to meat eating quality. However, other important traits such as meat texture, tenderness, and color have not much studied in cattle. Among these traits, meat tenderness and texture of cattle are among the most important factors affecting quality evaluation of consumers. Collagen is the main component of connective tissues.It greatly affects meat tenderness. The objective of this study was to determine significant variants and candidate genes associated with collagen contents trait (total collagen) through genome-wide association studies (GWAS). Phenotypic and genomic data from 135 Hanwoo were used. The BLUPF90 family program and GRAMMAR method for GWAS were applied in this study. A total of 73 potential single nucleotide polymorphisms (SNPs) showed significant associations with collagen content. They were located in or near 108 candidate genes. TMEM135 and ME3 genes were identified to have the most significant SNPs associated with collagen contents trait. Data indicated that these genes were related to collagen. Biological processes and pathways for the prediction of biological functions of candidate genes were confirmed. We found that candidate genes were involved in positive regulation of CREB transcription factor activity and actin cytoskeleton related to tenderness and texture of beef. Three genes (CRTC3, MYO1C and MYLK4) belonging to these biological functions were related to tenderness. These results provide a basis for improving genomic characteristics of Hanwoo for the production of tender beef. Furthermore, they could be used they could be used as an index to select desired traits for consumers.

Molecular modeling and dynamic simulation of chicken Mx protein with the S631N polymorphism.

Myxovirus resistance (Mx) proteins are antiviral GTPases induced by type I interferons (IFNs). In chickens, a single Mx protein variant, S631N, has been suggested to possess antiviral activity. However, the impact of this variant on chicken Mx (chMx) protein structure and conformation has not been investigated. Hence, in this study, we applied computational methods such as molecular modeling, molecular dynamic simulation, inter domain motion and residue networks to examine the structure and dynamic behavior of wild-type and mutant chMx. At first, we built 3-dimensional structural models for both wild-type and mutant chMx proteins, which revealed that the structural organization of chMx was similar to that of human Mx proteins. Subsequently, molecular dynamics simulations revealed that angle variation around the hinge1 region led to the different stalk domain conformations between the wild-type and mutant chMx proteins. Domain motion analysis further suggested that the conformational differences in the loop region surrounded by the mutant residue may lead to an inclined stalk domain conformation in the mutant compared to the wild-type protein. In addition, we performed betweenness centrality analysis from residue interaction networks, to identify the crucial residues for intramolecular signal flow in chMx. The results of this study provided information on the differences in structure and dynamics between wild-type and mutant chMx, which may aid in understanding the structural features of the S631N mutant, that may be associated with chMx protein antiviral activity.Communicated by Ramaswamy H. Sarma.

High-throughput sequencing-based metagenomic and transcriptomic analysis of intestine in piglets infected with salmonella.

Salmonella enterica serovar Typhimurium isolate HJL777 is a virulent bacterial strain in pigs. The high rate of salmonella infection are at high risk of non-typhoidal salmonella gastroenteritis development. Salmonellosis is most common in young pigs. We investigated changes in gut microbiota and biological function in piglets infected with salmonella via analysis of rectal fecal metagenome and intestinal transcriptome using 16S rRNA and RNA sequencing. We identified a decrease in Bacteroides and increase in harmful bacteria such as Spirochaetes and Proteobacteria by microbial community analysis. We predicted that reduction of Bacteroides by salmonella infection causes proliferation of salmonella and harmful bacteria that can cause an intestinal inflammatory response. Functional profiling of microbial communities in piglets with salmonella infection showed increasing lipid metabolism associated with proliferation of harmful bacteria and inflammatory responses. Transcriptome analysis identified 31 differentially expressed genes. Using gene ontology and Innate Immune Database analysis, we identified that BGN, DCN, ZFPM2 and BPI genes were involved in extracellular and immune mechanisms, specifically salmonella adhesion to host cells and inflammatory responses during infection. We confirmed alterations in gut microbiota and biological function during salmonella infection in piglets. Our findings will help prevent disease and improve productivity in the swine industry.

Antitumorigenic effect of insect-derived peptide poecilocorisin-1 in human skin cancer cells through regulation of Sp1 transcription factor.

Malignant melanoma is highly resistant to conventional treatments and is one of the most aggressive types of skin cancers. Conventional cancer treatments are limited due to drug resistance, tumor selectivity, and solubility. Therefore, new treatments with fewer side effects and excellent effects should be developed. In previous studies, we have analyzed antimicrobial peptides (AMPs), which showed antibacterial and anti-inflammatory effects in insects, and some AMPs also exhibited anticancer efficacy. Anticancer peptides (ACPs) are known to have fewer side effects and high anticancer efficacy. In this study, the insect-derived peptide poecilocorisin-1 (PCC-1) did not induce toxicity in the human epithelial cell line HaCaT, but its potential as an anticancer agent was confirmed through specific effects of antiproliferation, apoptosis, and cell cycle arrest in two melanoma cell lines, SK-MEL-28 and G361. Additionally, we discovered a novel anticancer mechanism of insect-derived peptides in melanoma through the regulation of transcription factor Sp1 protein, which is overexpressed in cancer, apoptosis, and cell cycle-related proteins. Taken together, this study aims to clarify the anticancer efficacy and safety of insect-derived peptides and to present their potential as future therapeutic agents.

Risk prediction and marker selection in nonsynonymous single nucleotide polymorphisms using whole genome sequencing data.

Despite the various existing studies about nonsynonymous single nucleotide polymorphisms (nsSNPs), genome-wide studies based on nsSNPs are rare. NsSNPs alter amino acid sequences, affect protein structure and function, and have deleterious effects. By predicting the deleterious effect of nsSNPs, we determined the total risk score per individual. Additionally, the machine learning technique was utilized to find an optimal nsSNP subset that best explains the complete nsSNP effect. A total of 16,100 nsSNPs were selected as the best representatives among 89,519 regressed nsSNPs. In the gene ontology analysis encompassing the 16,100 nsSNPs, DNA metabolic process, chemokine- and immune-related, and reproduction were the most enriched terms. We expect that our risk score prediction and nsSNP marker selection will contribute to future development of extant genome-wide association studies and breeding science more broadly.

Conformational Changes Induced by S34Y and R98C Variants in the Death Domain of Myd88.

Myeloid differentiating factor 88 (Myd88) is a universal adaptor protein that plays a critical role in innate immunity by mediating TLR downstream signaling. Myd88 death domain (DD) forms an oligomeric complex by association with other DD-containing proteins such as IRAK4. Despite its universal role, polymorphisms in Myd88 can result in several diseases. Previous studies have suggested that, out of several non-synonymous single-nucleotide polymorphisms (nsSNPs), the variants S34Y and R98C in the DD of Myd88 disrupt the formation of the Myddosome complex. Therefore, we performed molecular dynamics (MD) simulations on wild-type (Myd88WT) and mutant (Myd88S34Y, Myd88R98C) DDs to evaluate the subtle conformational changes induced by these mutations. Our results suggest that the S34Y variant induces large structural transitions compared to the R98C variant as evidenced by residual flexibility at the variable loop regions, particularly in the H1-H2 loop, and variations in the collective modes of motion observed for wild-type and mutant Myd88 DDs. The residue interaction network strongly suggests a distortion in the interaction pattern at the location of the mutated residue between the wild type and mutants. Moreover, betweenness centrality values indicate that variations in the distribution of functionally important residues may be reflected by distinct residue signal transductions in both wild-type and mutant Myd88 DDs, which may influence the interaction with other DDs in TLR downstream signaling.

Comparative genomics of canine Lactobacillus reuteri reveals adaptation to a shared environment with humans.

BACKGROUND: Lactobacillus reuteri is a gram-positive, non-motile bacterial species that has been used as a representative microorganism model to describe the ecology and evolution of vertebrate gut symbionts. OBJECTIVE: Because the genetic features and evolutionary strategies of L. reuteri from the gastrointestinal tract of canines remain unknown, we tried to construct draft genome canine L. reuteri and investigate modified, acquired, or lost genetic features that have facilitated the evolution and adaptation of strains to specific environmental niches by this study. METHODS: To examine canine L. reuteri, we sequenced an L. reuteri strain isolated from a dog in Korea. A comparative genomic approach was used to assess genetic diversity and gain insight into the distinguishing features related to different hosts based on 27 published genomic sequences. RESULTS: The pan-genome of 28 L. reuteri strains contained 7,369 gene families, and the core genome contained 1070 gene families. The ANI tree based on the core genes in the canine L. reuteri strain (C1) was very close to those for three strains (IRT, DSM20016, JCM1112) from humans. Evolutionarily, these four strains formed one clade, which we regarded as C1-clade in this study. We could investigate a total of 32,050 amino acid substitutions among the 28 L. reuteri strain genomes. In this comparison, 283 amino acid substitutions were specific to strain C1 and four strains in C1-clade shared most of these 283 C1-strain specific amino acid substitutions, suggesting strongly similar selective pressure. In accessory genes, we could identify 127 C1-clade host-specific genes and found that several genes were closely related to replication, recombination, and repair. CONCLUSION: This study provides new insights into the adaptation of L. reuteri to the canine intestinal habitat, and suggests that the genome of L. reuteri from canines is closely associated with their living and shared environment with humans.