Genome-wide association study using a single-step approach for teat number in Duroc, Landrace and Yorkshire pigs in Korea.

We investigated the genetic basis of teat number in sows, which is an important factor in their reproductive performance. We collected genotyping data from 20 353 pigs of three breeds (Duroc, Landrace and Yorkshire) using the Porcine SNP60K Bead Chip, and analyzed phenotypic data from 240 603 pigs. The heritability values of total teat number were 0.33 ± 0.02, 0.51 ± 0.01 and 0.50 ± 0.01 in Duroc, Landrace and Yorkshire pigs, respectively. A genome-wide association study was used to identify significant chromosomal regions associated with teat number in SSC7 and SSC9 in Duroc pig, SSC3, SSC7 and SSC18 in Landrace pig, and SSC7, SSC8 and SSC10 in Yorkshire pig. Among the markers, MARC0038565, located between the vertnin (VRTN) and synapse differentiation-inducing 1-like (SYNDIG1L) genes, showed the strongest association in the Duroc pig and was significant in all breeds. In Landrace and Yorkshire pigs, the most significant markers were located within the apoptosis resistant E3 ubiquitin protein ligase 1 (AREL1) and latent transforming growth factor beta-binding protein 2 (LTBP2) genes in SSC7, respectively. VRTN is a candidate gene regulating the teat number. Most markers were located in SSC7, indicating their significance in determining teat number and their potential as valuable genomic selection targets for improving this trait. Extensive linkage disequilibrium blocks were identified in SSC7, supporting their use in genomic selection strategies. Our study provides valuable insights into the genetic architecture of teat numbers in pigs, and helps identify candidate genes and genomic regions that may contribute to this economically important trait.

Pusillimonas thiosulfatoxidans sp. nov., a thiosulfate oxidizer isolated from activated sludge.

A Gram-stain-negative, motile bacterium, designated strain YE3T, was isolated from activated sludge obtained from a municipal wastewater treatment plant in Daejeon Metropolitan City, Republic of Korea. The cells were oxidase- and catalase-positive, and grew under aerobic conditions at 10-40 °C (optimum, 30 °C), with 1.0-8.0 % (w/v) NaCl (1.0 %) and at pH 5.5-9.0 (pH 7.0). Phylogenetic analysis based on the 16S rRNA gene sequence indicated that strain YE3T was most closely related to Pusillimonasharenae KACC 14927T (98.2 % sequence similarity) and Pusillimonasginsengisoli KCTC 22046T (98.0 %). DNA-DNA relatedness values for strain YE3T and P. harenae KACC 14927T, P. ginsengisoli KCTC 22046T and P. soli KCTC 22455T were 28.7±2.27 %, 21.3±1.16 %, and 14.0±0.67 %, respectively. The genomic G+C content of the type strain YE3T was 59.3 mol%, as determined by whole-genome sequencing. The dominant fatty acids were C16 : 0 (39.2 %) and C17 : 0cyclo (37.5 %). The major polar lipids of strain YE3T were diphosphatidylglycerol, phosphatidylglycerol and phosphatidylethanolamine. Two aminophospholipids and four unidentified lipids were also detected. Furthermore, strain YE3T was able to oxidize thiosulfate under heterotrophic conditions. Based on the phenotypic, genotypic, chemotaxonomic and phylogenetic analyses, strain YE3T represents a novel species of the genus Pusillimonas, for which the name Pusillimonas thiosulfatoxidans sp. nov. is proposed. The type strain is YE3T (=KCTC 62737T=NBRC 113113T).

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.

Effect of single nucleotide polymorphism on the total number of piglets born per parity of three different pig breeds.

OBJECTIVE: To determine the effects of genomic breeding values (GBV) and single nucleotide polymorphisms (SNP) on the total number of piglets born (TNB) in 3 pig breeds (Berkshire, Landrace, and Yorkshire). METHODS: After collecting genomic information (Porcine SNP BeadChip) and phenotypic TNB records for each breed, the effects of GBV and SNP were estimated by using single step best linear unbiased prediction (ssBLUP) method. RESULTS: The heritability estimates for TNB in Berkshire, Landrace, and Yorkshire breeds were 0.078, 0.107, and 0.121, respectively. The breeding value estimates for TNB in Berkshire, Landrace, and Yorkshire breeds were in the range of -1.34 to 1.47 heads, -1.79 to 1.87 heads, and -2.60 to 2.94 heads, respectively. Of sows having records for TNB, the reliability of breeding value for individuals with SNP information was higher than that for individuals without SNP information. Distributions of the SNP effects on TNB did not follow gamma distribution. Most SNP effects were near zero. Only a few SNPs had large effects. The numbers of SNPs with absolute value of more than 4 standard deviations in Berkshire, Landrace, and Yorkshire breeds were 11, 8, and 19, respectively. There was no SNP with absolute value of more than 5 standard deviations in Berkshire or Landrace. However, in Yorkshire, four SNPs (ASGA 0089457, ASGA0103374, ALGA0111816, and ALGA0098882) had absolute values of more than 5 standard deviations. CONCLUSION: There was no common SNP with large effect among breeds. This might be due to the large genetic composition differences and the small size of reference population. For the precise evaluation of genetic performance of individuals using a genomic selection method, it may be necessary to establish the appropriate size of reference population.

Exercise induced upregulation of glutamate-cysteine ligase catalytic subunit and glutamate-cysteine ligase modifier subunit gene expression in Thoroughbred horses.

OBJECTIVE: This study was performed to reveal the molecular structure and expression patterns of horse glutamate-cysteine ligase catalytic subunit (GCLC) and glutamate-cysteine ligase modifier subunit (GCLM) genes whose products form glutamate cysteine ligase, which were identified as differentially expressed genes in the previous study. METHODS: We performed bioinformatics analyses, and gene expression assay with quantitative polymerase chain reaction (qPCR) for horse GCLC and GCLM genes in muscle and blood leukocytes of Thoroughbred horses. RESULTS: Expression of GCLC showed the same pattern in both blood and muscle tissues after exercise. Expression of GCLC increased in the muscle and blood of Thoroughbreds, suggesting a tissue-specific regulatory mechanism for the expression of GCLC. In addition, expression of the GCLM gene increased after exercise in both the blood and muscle of Thoroughbreds. CONCLUSION: We established the expression patterns of GCLC and GCLM in the skeletal muscle and blood of Thoroughbred horses in response to exercise. Further study is now warranted to uncover the functional importance of these genes in exercise and recovery in racehorses.

Genome-wide analysis of DNA methylation patterns in horse.

BACKGROUND: DNA methylation is an epigenetic regulatory mechanism that plays an essential role in mediating biological processes and determining phenotypic plasticity in organisms. Although the horse reference genome and whole transcriptome data are publically available the global DNA methylation data are yet to be known. RESULTS: We report the first genome-wide DNA methylation characteristics data from skeletal muscle, heart, lung, and cerebrum tissues of thoroughbred (TH) and Jeju (JH) horses, an indigenous Korea breed, respectively by methyl-DNA immunoprecipitation sequencing. The analysis of the DNA methylation patterns indicated that the average methylation density was the lowest in the promoter region, while the density in the coding DNA sequence region was the highest. Among repeat elements, a relatively high density of methylation was observed in long interspersed nuclear elements compared to short interspersed nuclear elements or long terminal repeat elements. We also successfully identified differential methylated regions through a comparative analysis of corresponding tissues from TH and JH, indicating that the gene body regions showed a high methylation density. CONCLUSIONS: We provide report the first DNA methylation landscape and differentially methylated genomic regions (DMRs) of thoroughbred and Jeju horses, providing comprehensive DMRs maps of the DNA methylome. These data are invaluable resource to better understanding of epigenetics in the horse providing information for the further biological function analyses.

Thoroughbred Horse Single Nucleotide Polymorphism and Expression Database: HSDB.

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

Characteristics of Purified Horse Oil by Supercritical Fluid Extraction with Different Deodorants Agents.

This study investigated the impact of activated carbon, palm activated carbon, and zeolite on horse oil (HO) extracted from horse neck fat using supercritical fluid extraction with deodorant-untreated HO (CON) as a comparison. The yield and lipid oxidation of deodorant untreated HO (CON) were not significantly affected by the three deodorants. However, deodorant-treated HOs exhibited significantly elevated levels of α-linolenic acid (C18:3n3) and eicosenoic acid (C20:1n9) compared to CON (p<0.05), while other fatty acids remained consistent. Zeolite-purified HO demonstrated significantly lower levels of volatile organic compounds (VOCs) than other treatments (p<0.05). Remarkably, zeolite decreased the concentration of pentane, 2,3-dimethyl (gasoline odor), by over 90%, from 177.17 A.U. ×106 in CON to 15.91 A.U. ×106. Zeolite also effectively eliminates sec-butylamine (ammonia and fishy odor) as compared to other deodorant-treated HOs (p<0.05). Additionally, zeolite reduced VOCs associated with the fruity citrus flavor, such as nonanal, octanal, and D-limonene in HO (p<0.05). This study suggests that integrating zeolite in supercritical fluid extraction enhances HO purification by effectively eliminating undesirable VOCs, presenting a valuable approach for producing high-quality HO production in the cosmetic and functional food industries.

Whole transcriptome analyses of six thoroughbred horses before and after exercise using RNA-Seq.

BACKGROUND: Thoroughbred horses are the most expensive domestic animals, and their running ability and knowledge about their muscle-related diseases are important in animal genetics. While the horse reference genome is available, there has been no large-scale functional annotation of the genome using expressed genes derived from transcriptomes. RESULTS: We present a large-scale analysis of whole transcriptome data. We sequenced the whole mRNA from the blood and muscle tissues of six thoroughbred horses before and after exercise. By comparing current genome annotations, we identified 32,361 unigene clusters spanning 51.83 Mb that contained 11,933 (36.87%) annotated genes. More than 60% (20,428) of the unigene clusters did not match any current equine gene model. We also identified 189,973 single nucleotide variations (SNVs) from the sequences aligned against the horse reference genome. Most SNVs (171,558 SNVs; 90.31%) were novel when compared with over 1.1 million equine SNPs from two SNP databases. Using differential expression analysis, we further identified a number of exercise-regulated genes: 62 up-regulated and 80 down-regulated genes in the blood, and 878 up-regulated and 285 down-regulated genes in the muscle. Six of 28 previously-known exercise-related genes were over-expressed in the muscle after exercise. Among the differentially expressed genes, there were 91 transcription factor-encoding genes, which included 56 functionally unknown transcription factor candidates that are probably associated with an early regulatory exercise mechanism. In addition, we found interesting RNA expression patterns where different alternative splicing forms of the same gene showed reversed expressions before and after exercising. CONCLUSION: The first sequencing-based horse transcriptome data, extensive analyses results, deferentially expressed genes before and after exercise, and candidate genes that are related to the exercise are provided in this study.

Analyses of porcine public SNPs in coding-gene regions by re-sequencing and phenotypic association studies.

The Porcine SNP database has a huge number of SNPs, but these SNPs are mostly found by computer data-mining procedures and have not been well characterized. We re-sequenced 1,439 porcine public SNPs from four commercial pig breeds and one Korean domestic breed (Korean Native pig, KNP) by using two DNA pools from eight unrelated animals in each breed. These SNPs were from 419 protein-coding genes covering the 18 autosomes, and the re-sequencing in breeds confirmed 690 public SNPs (47.9%) and 226 novel mutations (173 SNPs and 53 insertions/deletions). Thus, totally, 916 variations were found from our study. Of the 916 variations, 148 SNPs (16.2%) were found across all the five breeds, and 199 SNPs (21.7%) were breed specific polymorphisms. According to the SNP locations in the gene sequences, these 916 variations were categorized into 802 non-coding SNPs (785 in intron, 17 in 3'-UTR) and 114 coding SNPs (86 synonymous SNPs, 28 non-synonymous SNPs). The nucleotide substitution analyses for these SNPs revealed that 70.2% were from transitions, 20.0% from transversions, and the remaining 5.79% were deletions or insertions. Subsequently, we genotyped 261 SNPs from 180 genes in an experimental KNP × Landrace F2 cross by the Sequenom MassARRAY system. A total of 33 traits including growth, carcass composition and meat quality were analyzed for the phenotypic association tests using the 132 SNPs in 108 genes with minor allele frequency (MAF)>0.2. The association results showed that five marker-trait combinations were significant at the 5% experiment-wise level (ADCK4 for rear leg, MYH3 for rear leg, Hunter B, Loin weight and Shearforce) and four at the 10% experiment-wise level (DHX38 for average daily gain at live weight, LGALS9 for crude lipid, NGEF for front leg and LIFR for pH at 24 h). In addition, 49 SNPs in 44 genes showing significant association with the traits were detected at the 1% comparison-wise level. A large number of genes that function as enzymes, transcription factors or signalling molecules were considered as genetic markers for pig growth (RNF103, TSPAN31, DHX38, ABCF1, ABCC10, SCD5, KIAA0999 and FKBP10), muscling (HSPA5, PTPRM, NUP88, ADCK4, PLOD1, DLX1 and GRM8), fatness (PTGIS, IDH3B, RYR2 and NOL4) and meat quality traits (DUSP4, LIFR, NGEF, EWSR1, ACTN2, PLXND1, DLX3, LGALS9, ENO3, EPRS, TRIM29, EHMT2, RBM42, SESN2 and RAB4B). The SNPs or genes reported here may be beneficial to future marker assisted selection breeding in pigs.

Complete genome of Halomonas aestuarii Hb3, isolated from tidal flat.

Halomonas aestuarii Hb3, a moderately halophilic bacterium belonging to the class Gammaproteobacteria, was isolated from a tidal flat. Herein, we report the complete genome sequence of its strain Hb3. Its size is estimated at 3.54Mbp with a mean G+C content of 67.9%. The genome includes 3238 open reading frames, 65 transfer RNAs, and four ribosomal RNA gene operons. Genes related to the degradation of monoaromatic compounds, detoxification of arsenic, and production of polymers were identified. These features indicate that this strain may be important for ecological and industrial application.

Genetic Diversity and Population Structure of Ethiopian Sheep Populations Revealed by High-Density SNP Markers.

Sheep in Ethiopia are adapted to a wide range of environments, including extreme habitats. Elucidating their genetic diversity is critical for improving breeding strategies and mapping quantitative trait loci associated with productivity. To this end, the present study investigated the genetic diversity and population structure of five Ethiopian sheep populations exhibiting distinct phenotypes and sampled from distinct production environments, including arid lowlands and highlands. To investigate the genetic relationships in greater detail and infer population structure of Ethiopian sheep breeds at the continental and global levels, we analyzed genotypic data of selected sheep breeds from the Ovine SNP50K HapMap dataset. All Ethiopian sheep samples were genotyped with Ovine Infinium HD SNP BeadChip (600K). Mean genetic diversity ranged from 0.29 in Arsi-Bale to 0.32 in Menz sheep, while estimates of genetic differentiation among populations ranged from 0.02 to 0.07, indicating low to moderate differentiation. An analysis of molecular variance revealed that 94.62 and 5.38% of the genetic variation was attributable to differences within and among populations, respectively. Our population structure analysis revealed clustering of five Ethiopian sheep populations according to tail phenotype and geographic origin-i.e., short fat-tailed (very cool high-altitude), long fat-tailed (mid to high-altitude), and fat-rumped (arid low-altitude), with clear evidence of admixture between long fat-tailed populations. North African sheep breeds showed higher levels of within-breed diversity, but were less differentiated than breeds from Eastern and Southern Africa. When African breeds were grouped according to geographic origin (North, South, and East), statistically significant differences were detected among groups (regions). A comparison of population structure between Ethiopian and global sheep breeds showed that fat-tailed breeds from Eastern and Southern Africa clustered together, suggesting that these breeds were introduced to the African continent via the Horn and migrated further south.

Comparative Transcriptomic Analyses by RNA-seq to Elucidate Differentially Expressed Genes in the Muscle of Korean Thoroughbred Horses.

The athletic abilities of the horse serve as a valuable model to understand the physiology and molecular mechanisms of adaptive responses to exercise. We analyzed differentially expressed genes in triceps brachii muscle tissues collected from Eonjena Taeyang and Jigusang Seryeok Thoroughbred horses and their co-expression networks in a large-scale RNA-sequence dataset comparing expression before and after exercise. High-quality horse transcriptome data were generated, with over 22 million 90-bp pair-end reads. By comparing the annotations, we found that MYH3, MPZ, and PDE8B genes in Eonjena Taeyang and PDE8B and KIF18A genes in Jigusang Seryeok were upregulated before exercise. Notably further, we observed that PPP1R27, NDUFA3, TNC, and ANK1 in Eonjena Taeyang and HIF1A, BDNF, ADRB2, OBSCN, and PER3 in Jigusang Seryeok have shown upregulation at the postexercise period. This investigation suggested that genes responsible for metabolism and oxidative phosphorylations associated with endurance and resistance exercise were highly expressed, whereas genes encoding structural proteins were generally suppressed. The expression profile of racehorses at pre- and postexercise will provide credible reference for further studies on biological effects such as responses to stress and adaption of other Thoroughbred horse, which might be useful for selective breeding for improvement of traits in commercial production.

Genome-wide analysis of DNA methylation before-and after exercise in the thoroughbred horse with MeDIP-Seq.

Athletic performance is an important criteria used for the selection of superior horses. However, little is known about exercise-related epigenetic processes in the horse. DNA methylation is a key mechanism for regulating gene expression in response to environmental changes. We carried out comparative genomic analysis of genome-wide DNA methylation profiles in the blood samples of two different thoroughbred horses before and after exercise by methylated-DNA immunoprecipitation sequencing (MeDIP-Seq). Differentially methylated regions (DMRs) in the pre-and post-exercise blood samples of superior and inferior horses were identified. Exercise altered the methylation patterns. After 30 min of exercise, 596 genes were hypomethylated and 715 genes were hypermethylated in the superior horse, whereas in the inferior horse, 868 genes were hypomethylated and 794 genes were hypermethylated. These genes were analyzed based on gene ontology (GO) annotations and the exercise-related pathway patterns in the two horses were compared. After exercise, gene regions related to cell division and adhesion were hypermethylated in the superior horse, whereas regions related to cell signaling and transport were hypermethylated in the inferior horse. Analysis of the distribution of methylated CpG islands confirmed the hypomethylation in the gene-body methylation regions after exercise. The methylation patterns of transposable elements also changed after exercise. Long interspersed nuclear elements (LINEs) showed abundance of DMRs. Collectively, our results serve as a basis to study exercise-based reprogramming of epigenetic traits.

Estimation of effective population size using single-nucleotide polymorphism (SNP) data in Jeju horse.

This study was conducted to estimate the effective population size using SNPs data of 240 Jeju horses that had raced at the Jeju racing park. Of the total 61,746 genotyped autosomal SNPs, 17,320 (28.1%) SNPs (missing genotype rate of >10%, minor allele frequency of <0.05 and Hardy-Weinberg equilibrium test P-value of <10(-6)) were excluded after quality control processes. SNPs on the X and Y chromosomes and genotyped individuals with missing genotype rate over 10% were also excluded, and finally, 44,426 (71.9%) SNPs were selected and used for the analysis. The measures of the LD, square of correlation coefficient (r(2)) between SNP pairs, were calculated for each allele and the effective population size was determined based on r(2) measures. The polymorphism information contents (PIC) and expected heterozygosity (HE) were 0.27 and 0.34, respectively. In LD, the most rapid decline was observed over the first 1 Mb. But r(2) decreased more slowly with increasing distance and was constant after 2 Mb of distance and the decline was almost linear with log-transformed distance. The average r(2) between adjacent SNP pairs ranged from 0.20 to 0.31 in each chromosome and whole average was 0.26, while the whole average r(2) between all SNP pairs was 0.02. We observed an initial pattern of decreasing Ne and estimated values were closer to 41 at 1 ~ 5 generations ago. The effective population size (41 heads) estimated in this study seems to be large considering Jeju horse's population size (about 2,000 heads), but it should be interpreted with caution because of the technical limitations of the methods and sample size.

Peeling back the evolutionary layers of molecular mechanisms responsive to exercise-stress in the skeletal muscle of the racing horse.

The modern horse (Equus caballus) is the product of over 50 million yrs of evolution. The athletic abilities of the horse have been enhanced during the past 6000 yrs under domestication. Therefore, the horse serves as a valuable model to understand the physiology and molecular mechanisms of adaptive responses to exercise. The structure and function of skeletal muscle show remarkable plasticity to the physical and metabolic challenges following exercise. Here, we reveal an evolutionary layer of responsiveness to exercise-stress in the skeletal muscle of the racing horse. We analysed differentially expressed genes and their co-expression networks in a large-scale RNA-sequence dataset comparing expression before and after exercise. By estimating genome-wide dN/dS ratios using six mammalian genomes, and FST and iHS using re-sequencing data derived from 20 horses, we were able to peel back the evolutionary layers of adaptations to exercise-stress in the horse. We found that the oldest and thickest layer (dN/dS) consists of system-wide tissue and organ adaptations. We further find that, during the period of horse domestication, the older layer (FST) is mainly responsible for adaptations to inflammation and energy metabolism, and the most recent layer (iHS) for neurological system process, cell adhesion, and proteolysis.