Genome-wide association study for carcass traits in Hanwoo cattle using additional relatives' information of non-genotyped animals.

The aim of this study was to find significant genomic regions associated with carcass traits in Hanwoo cattle and to compare the benefit of using additional information from non-genotyped animals. Imputed whole-genome sequence data were used along with phenotypic data on 13 715 genotyped animals as well as phenotypes of 440 284 non-genotyped animals that were offspring of 454 genotyped sires. For carcass weight, 15 083 SNPs in 33 QTL regions and 313 candidate genes were identified. We found 410 SNPs in 17 QTL regions containing 122 candidate genes for back fat thickness. In total, 656 SNPs in 19 QTLs with 137 candidate genes for eye muscle area and 79 SNPs in 12 QTL regions with 77 candidate genes were identified for marbling score. The most important candidate genes included ZFAT, TG, PLAG1, CHCHD7, and TOX for carcass weight and eye muscle area, NOG for back fat thickness, and EVOVL5 for marbling score. This study showed that the use of phenotypic records on non-genotyped progeny along with imputed whole-genome sequence data increased the power of detecting new significant genomic regions.

Brachygnathia Inferior in Cloned Dogs Is Possibly Correlated with Variants of Wnt Signaling Pathway Initiators.

Abnormalities in animals cloned via somatic cell nuclear transfer (SCNT) have been reported. In this study, to produce bomb-sniffing dogs, we successfully cloned four healthy dogs through SCNT using the same donor genome from the skin of a male German shepherd old dog. Veterinary diagnosis (X-ray/3D-CT imaging) revealed that two cloned dogs showed normal phenotypes, whereas the others showed abnormal shortening of the mandible (brachygnathia inferior) at 1 month after birth, even though they were cloned under the same conditions except for the oocyte source. Therefore, we aimed to determine the genetic cause of brachygnathia inferior in these cloned dogs. To determine the genetic defects related to brachygnathia inferior, we performed karyotyping and whole-genome sequencing (WGS) for identifying small genetic alterations in the genome, such as single-nucleotide variations or frameshifts. There were no chromosomal numerical abnormalities in all cloned dogs. However, WGS analysis revealed variants of Wnt signaling pathway initiators (WNT5B, DVL2, DACT1, ARRB2, FZD 4/8) and cadherin (CDH11, CDH1like) in cloned dogs with brachygnathia inferior. In conclusion, this study proposes that brachygnathia inferior in cloned dogs may be associated with variants in initiators and/or regulators of the Wnt/cadherin signaling pathway.

The transcriptomic blueprint of molt in rooster using various tissues from Ginkkoridak (Korean long-tailed chicken).

BACKGROUND: Annual molt is a critical stage in the life cycle of birds. Although the most extensively documented aspects of molt are the renewing of plumage and the remodeling of the reproductive tract in laying hens, in chicken, molt deeply affects various tissues and physiological functions. However, with exception of the reproductive tract, the effect of molt on gene expression across the tissues known to be affected by molt has to date never been investigated. The present study aimed to decipher the transcriptomic effects of molt in Ginkkoridak, a Korean long-tailed chicken. Messenger RNA data available across 24 types of tissue samples (9 males) and a combination of mRNA and miRNA data on 10 males and 10 females blood were used. RESULTS: The impact of molt on gene expression and gene transcript usage appeared to vary substantially across tissues types in terms of histological entities or physiological functions particularly related to nervous system. Blood was the tissue most affected by molt in terms of differentially expressed genes in both sexes, closely followed by meninges, bone marrow and heart. The effect of molt in blood appeared to differ between males and females, with a more than fivefold difference in the number of down-regulated genes between both sexes. The blueprint of molt in roosters appeared to be specific to tissues or group of tissues, with relatively few genes replicating extensively across tissues, excepted for the spliceosome genes (U1, U4) and the ribosomal proteins (RPL21, RPL23). By integrating miRNA and mRNA data, when chickens molt, potential roles of miRNA were discovered such as regulation of neurogenesis, regulation of immunity and development of various organs. Furthermore, reliable candidate biomarkers of molt were found, which are related to cell dynamics, nervous system or immunity, processes or functions that have been shown to be extensively modulated in response to molt. CONCLUSIONS: Our results provide a comprehensive description at the scale of the whole organism deciphering the effects of molt on the transcriptome in chicken. Also, the conclusion of this study can be used as a valuable resource in transcriptome analyses of chicken in the future and provide new insights related to molt.

Comparative methylation and RNA-seq expression analysis in CpG context to identify genes involved in Backfat vs. Liver diversification in Nanchukmacdon Pig.

BACKGROUND: DNA methylation and demethylation at CpG islands is one of the main regulatory factors that allow cells to respond to different stimuli. These regulatory mechanisms help in developing tissue without affecting the genomic composition or undergoing selection. Liver and backfat play important roles in regulating lipid metabolism and control various pathways involved in reproductive performance, meat quality, and immunity. Genes inside these tissue store a plethora of information and an understanding of these genes is required to enhance tissue characteristics in the future generation. RESULTS: A total of 16 CpG islands were identified, and they were involved in differentially methylation regions (DMRs) as well as differentially expressed genes (DEGs) of liver and backfat tissue samples. The genes C7orf50, ACTB and MLC1 in backfat and TNNT3, SIX2, SDK1, CLSTN3, LTBP4, CFAP74, SLC22A23, FOXC1, GMDS, GSC, GATA4, SEMA5A and HOXA5 in the liver, were categorized as differentially-methylated. Subsequently, Motif analysis for DMRs was performed to understand the role of the methylated motif for tissue-specific differentiation. Gene ontology studies revealed association with collagen fibril organization, the Bone Morphogenetic Proteins (BMP) signaling pathway in backfat and cholesterol biosynthesis, bile acid and bile salt transport, and immunity-related pathways in methylated genes expressed in the liver. CONCLUSIONS: In this study, to understand the role of genes in the differentiation process, we have performed whole-genome bisulfite sequencing (WGBS) and RNA-seq analysis of Nanchukmacdon pigs. Methylation and motif analysis reveals the critical role of CpG islands and transcriptional factors binding site (TFBS) in guiding the differential patterns. Our findings could help in understanding how methylation of certain genes plays an important role and can be used as biomarkers to study tissue specific characteristics.

Whole metagenome sequencing of cecum microbiomes in Ethiopian indigenous chickens from two different altitudes reveals antibiotic resistance genes.

We analyzed the whole genomes of cecum microbiomes of Ethiopian indigenous chickens from two distinct geographical zones: Afar (AF) district (Dulecha, 730 m above sea level) and Amhara (AM) district (Menz Gera Midir, 3300 m). Through metagenomic analysis we found that microbial populations were mainly dominated by Bacteroidetes and Firmicutes. We identified 2210 common genes in the two groups. LEfSe showed that the distribution of Coprobacter, Geobacter, Cronobacter, Alloprevotella, and Dysgonomonas were more abundant in AF than AM. Analyses using KEGG, eggNOG, and CAZy databases indicated that the pathways of metabolism, genetic information processing, environmental information processing, and cellular process were significantly enriched. Functional abundance was found to be associated with the nutrient absorption and microbial localization of indigenous chickens. We also investigated antibiotic resistant genes and found antibiotics like LSM, cephalosporin, and tetracycline were significantly more abundant in AF than AM.

Population analysis of the Korean native duck using whole-genome sequencing data.

BACKGROUND: Advances in next-generation sequencing technologies have provided an opportunity to perform population-level comparative genomic analysis to discover unique genomic characteristics of domesticated animals. Duck is one of the most popular domesticated waterfowls, which is economically important as a source of meat, eggs, and feathers. The objective of this study is to perform population and functional analyses of Korean native duck, which has a distinct meat flavor and texture phenotype, using whole-genome sequencing data. To study the distinct genomic features of Korean native duck, we conducted population-level genomic analysis of 20 Korean native ducks together with 15 other duck breeds. RESULTS: A total of 15.56 million single nucleotide polymorphisms were detected in Korean native duck. Based on the unique existence of non-synonymous single nucleotide polymorphisms in Korean native duck, a total of 103 genes related to the unique genomic characteristics of Korean native duck were identified in comparison with 15 other duck breeds, and their functions were investigated. The nucleotide diversity and population structures among the used duck breeds were then compared, and their phylogenetic relationship was analyzed. Finally, highly differentiated genomic regions among Korean native duck and other duck breeds were identified, and functions of genes in those regions were examined. CONCLUSIONS: This is the first study to compare the population of Korean native duck with those of other duck breeds by using whole-genome sequencing data. Our findings can be used to expand our knowledge of genomic characteristics of Korean native duck, and broaden our understanding of duck breeds.

Serial gene co-expression network approach to mine biological meanings from integrated transcriptomes of the porcine endometrium during estrous cycle.

The estrous cycle is a complex process regulated by several hormones. To understand the dynamic changes in gene expression that takes place in the swine endometrium during the estrous cycle relative to the day of estrus onset, we performed RNA-sequencing analysis on days 0, 3, 6, 9, 12, 15, and 18, resulting in the identification of 4495 differentially expressed genes (DEGs; Q ≤ 0.05 and |log2FC| ≥ 1) at various phases in the estrous cycle. These DEGs were integrated into multiple gene co-expression networks based on different fold changes and correlation coefficient (R2) thresholds and a suitable network, which included 899 genes (|log2FC| ≥ 2 and R2 ≥ 0.99), was identified for downstream analyses based on the biological relevance of the Gene Ontology (GO) terms enriched. The genes in this network were partitioned into 6 clusters based on the expression pattern. Several GO terms including cell cycle, apoptosis, hormone signaling, and lipid biosynthetic process were found to be enriched. Furthermore, we found 15 significant KEGG pathways, including cell adhesion molecules, cytokine-cytokine receptor signaling, steroid biosynthesis, and estrogen signaling pathways. We identified several genes and GO terms to be stage-specific. Moreover, the identified genes and pathways extend our understanding of porcine endometrial regulation during estrous cycle and will serve as a good resource for future studies.

Use of gene expression and whole-genome sequence information to improve the accuracy of genomic prediction for carcass traits in Hanwoo cattle.

BACKGROUND: In this study, we assessed the accuracy of genomic prediction for carcass weight (CWT), marbling score (MS), eye muscle area (EMA) and back fat thickness (BFT) in Hanwoo cattle when using genomic best linear unbiased prediction (GBLUP), weighted GBLUP (wGBLUP), and a BayesR model. For these models, we investigated the potential gain from using pre-selected single nucleotide polymorphisms (SNPs) from a genome-wide association study (GWAS) on imputed sequence data and from gene expression information. We used data on 13,717 animals with carcass phenotypes and imputed sequence genotypes that were split in an independent GWAS discovery set of varying size and a remaining set for validation of prediction. Expression data were used from a Hanwoo gene expression experiment based on 45 animals. RESULTS: Using a larger number of animals in the reference set increased the accuracy of genomic prediction whereas a larger independent GWAS discovery dataset improved identification of predictive SNPs. Using pre-selected SNPs from GWAS in GBLUP improved accuracy of prediction by 0.02 for EMA and up to 0.05 for BFT, CWT, and MS, compared to a 50 k standard SNP array that gave accuracies of 0.50, 0.47, 0.58, and 0.47, respectively. Accuracy of prediction of BFT and CWT increased when BayesR was applied with the 50 k SNP array (0.02 and 0.03, respectively) and was further improved by combining the 50 k array with the top-SNPs (0.06 and 0.04, respectively). By contrast, using BayesR resulted in limited improvement for EMA and MS. wGBLUP did not improve accuracy but increased prediction bias. Based on the RNA-seq experiment, we identified informative expression quantitative trait loci, which, when used in GBLUP, improved the accuracy of prediction slightly, i.e. between 0.01 and 0.02. SNPs that were located in genes, the expression of which was associated with differences in trait phenotype, did not contribute to a higher prediction accuracy. CONCLUSIONS: Our results show that, in Hanwoo beef cattle, when SNPs are pre-selected from GWAS on imputed sequence data, the accuracy of prediction improves only slightly whereas the contribution of SNPs that are selected based on gene expression is not significant. The benefit of statistical models to prioritize selected SNPs for estimating genomic breeding values is trait-specific and depends on the genetic architecture of each trait.

Genome sequencing and protein domain annotations of Korean Hanwoo cattle identify Hanwoo-specific immunity-related and other novel genes.

BACKGROUND: Identification of genetic mechanisms and idiosyncrasies at the breed-level can provide valuable information for potential use in evolutionary studies, medical applications, and breeding of selective traits. Here, we analyzed genomic data collected from 136 Korean Native cattle, known as Hanwoo, using advanced statistical methods. RESULTS: Results revealed Hanwoo-specific protein domains which were largely characterized by immunoglobulin function. Furthermore, domain interactions of novel Hanwoo-specific genes reveal additional links to immunity. Novel Hanwoo-specific genes linked to muscle and other functions were identified, including protein domains with functions related to energy, fat storage, and muscle function that may provide insight into the mechanisms behind Hanwoo cattle's uniquely high percentage of intramuscular fat and fat marbling. CONCLUSION: The identification of Hanwoo-specific genes linked to immunity are potentially useful for future medical research and selective breeding. The significant genomic variations identified here can crucially identify genetic novelties that are arising from useful adaptations. These results will allow future researchers to compare and classify breeds, identify important genetic markers, and develop breeding strategies to further improve significant traits.

Non-Coding Transcriptome Maps across Twenty Tissues of the Korean Black Chicken, Yeonsan Ogye.

Yeonsan Ogye is a rare Korean domestic chicken breed whose entire body, including feathers and skin, has a unique black coloring. Although some protein-coding genes related to this unique feature have been examined, non-coding elements have not been widely investigated. Thus, we evaluated coding and non-coding transcriptome expression and identified long non-coding RNAs functionally linked to protein-coding genes in Ogye. High-throughput RNA sequencing and DNA methylation sequencing were performed to profile the expression of 14,264 Ogye protein-coding and 6900 long non-coding RNA (lncRNA) genes and detect DNA methylation in 20 different tissues of an individual Ogye. Approximately 75% of Ogye lncRNAs and 45% of protein-coding genes showed tissue-specific expression. For some genes, tissue-specific expression levels were inversely correlated with DNA methylation levels in their promoters. Approximately 39% of tissue-specific lncRNAs displayed functional associations with proximal or distal protein-coding genes. Heat shock transcription factor 2-associated lncRNAs appeared to be functionally linked to protein-coding genes specifically expressed in black skin tissues, more syntenically conserved in mammals, and differentially expressed in black relative to in white tissues. Pending experimental validation, our findings increase the understanding of how the non-coding genome regulates unique phenotypes and can be used for future genomic breeding of chickens.

Cattle genome-wide analysis reveals genetic signatures in trypanotolerant N'Dama.

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.

Accuracy of genomic breeding value prediction for intramuscular fat using different genomic relationship matrices in Hanwoo (Korean cattle).

OBJECTIVE: Intramuscular fat is one of the meat quality traits that is considered in the selection strategies for Hanwoo (Korean cattle). Different methods are used to estimate the breeding value of selection candidates. In the present work we focused on accuracy of different genotype relationship matrices as described by forni and pedigree based relationship matrix. METHODS: The data set included a total of 778 animals that were genotyped for BovineSNP50 BeadChip. Among these 778 animals, 72 animals were sires for 706 reference animals and were used as a validation dataset. Single trait animal model (best linear unbiased prediction and genomic best linear unbiased prediction) was used to estimate the breeding values from genomic and pedigree information. RESULTS: The diagonal elements for the pedigree based coefficients were slightly higher for the genomic relationship matrices (GRM) based coefficients while off diagonal elements were considerably low for GRM based coefficients. The accuracy of breeding value for the pedigree based relationship matrix (A) was 13% while for GRM (GOF, G05, and Yang) it was 0.37, 0.45, and 0.38, respectively. CONCLUSION: Accuracy of GRM was 1.5 times higher than A in this study. Therefore, genomic information will be more beneficial than pedigree information in the Hanwoo breeding program.

A genome-wide assessment of genetic diversity and population structure of Korean native cattle breeds.

BACKGROUND: The native cattle breeds are an important genetic resource for meat and milk production throughout Asia. In Asia cattle were domesticated around 10,000 years ago and in Korea cattle are being raised since 2000 B.C. There are three native breeds of cattle in Korea viz. Brown Hanwoo, Brindle Hanwoo and Jeju Black. While one of these breeds, Brown Hanwoo, is a part of a Food and Agricultural Organization and national genetic evaluation plans, others get little attention. This study is an effort to understand and provide a detailed insight into the population structure and genetic variability of the Korean cattle breeds along with other Asian breeds using various methods. In this study we report the genetic variation and structure of the Korean cattle breeds and their comparison with five other Asian cattle breeds along with a panel of animals from European taurine, African taurine and indicine cattle breeds. RESULTS: Asian cattle were found to be least differentiated which reflects their recent history. Amongst the Asian breeds Hainan, which is an indicine breed, had the lowest gene diversity while Yanbian had the highest followed by Mongolian and Korean cattle. Amongst the Korean breeds Brown Hanwoo had the highest diversity followed by Brindle Hanwoo and Jeju Black. The genetic diversity in Asian cattle breeds was found comparable to the European taurines and more than the African taurines and Zebu cattle. Korean cattle breed, Brown Hanwoo was consistently found to be closer to Yanbian, a Chinese cattle breed. We found low divergence and moderate levels of genetic diversity among the native Korean breeds. Indicine introgression from Hainan was seen in other Asian breeds. From Europe, Limousin, Holstein and Hereford introgression was found in Asian breeds. CONCLUSIONS: In this study we provide a genome-wide insight into the genetic history of the native cattle breeds of Korea. The outcomes of this study will help in prioritization and designing of the conservation plans.

Gene Expression Patterns Associated with Peroxisome Proliferator-activated Receptor (PPAR) Signaling in the Longissimus dorsi of Hanwoo (Korean Cattle).

Adipose tissue deposited within muscle fibers, known as intramuscular fat (IMF or marbling), is a major determinant of meat quality and thereby affects its economic value. The biological mechanisms that determine IMF content are therefore of interest. In this study, 48 genes involved in the bovine peroxisome proliferator-activated receptor signaling pathway, which is involved in lipid metabolism, were investigated to identify candidate genes associated with IMF in the longissimus dorsi of Hanwoo (Korean cattle). Ten genes, retinoid X receptor alpha, peroxisome proliferator-activated receptor gamma (PPARG), phospholipid transfer protein, stearoyl-CoA desaturase, nuclear receptor subfamily 1 group H member 3, fatty acid binding protein 3 (FABP3), carnitine palmitoyltransferase II, acyl-Coenzyme A dehydrogenase long chain (ACADL), acyl-Coenzyme A oxidase 2 branched chain, and fatty acid binding protein 4, showed significant effects with regard to IMF and were differentially expressed between the low- and high-marbled groups (p<0.05). Analysis of the gene co-expression network based on Pearson's correlation coefficients identified 10 up-regulated genes in the high-marbled group that formed a major cluster. Among these genes, the PPARG-FABP4 gene pair exhibited the strongest correlation in the network. Glycerol kinase was found to play a role in mediating activation of the differentially expressed genes. We categorized the 10 significantly differentially expressed genes into the corresponding downstream pathways and investigated the direct interactive relationships among these genes. We suggest that fatty acid oxidation is the major downstream pathway affecting IMF content. The PPARG/RXRA complex triggers activation of target genes involved in fatty acid oxidation resulting in increased triglyceride formation by ATP production. Our findings highlight candidate genes associated with the IMF content of the loin muscle of Korean cattle and provide insight into the biological mechanisms that determine adipose deposition within muscle.

Homology modeling study of bovine µ-calpain inhibitor-binding domains.

The activated mammalian CAPN-structures, the CAPN/CAST complex in particular, have become an invaluable target model using the structure-based virtual screening of drug candidates from the discovery phase to development for over-activated CAPN linked to several diseases, such as post-ischemic injury and cataract formation. The effect of Ca²âº-binding to the enzyme is thought to include activation, as well as the dissociation, aggregation, and autolysis of small regular subunits. Unfortunately, the Ca²âº-activated enzyme tends to aggregate when provided as a divalent ion at the high-concentration required for the protease crystallization. This is also makes it very difficult to crystallize the whole-length enzyme itself, as well as the enzyme-inhibitor complex. Several parameters that influence CAPN activity have been investigated to determine its roles in Ca²âº-modulation, autoproteolysis, phosphorylation, and intracellular distribution and inhibition by its endogenous inhibitor CAST. CAST binds and inhibits CAPN via its CAPN-inhibitor domains (four repeating domains 1-4; CAST1-4) when CAPN is activated by Ca²âº-binding. An important key to understanding CAPN1 inhibition by CAST is to determine how CAST interacts at the molecular level with CAPN1 to inhibit its protease activity. In this study, a 3D structure model of a CAPN1 bound bovine CAST4 complex was built by comparative modeling based on the only known template structure of a rat CAPN2/CAST4 complex. The complex model suggests certain residues of bovine CAST4, notably, the TIPPKYQ motif sequence, and the structural elements of these residues, which are important for CAPN1 inhibition. In particular, as CAST4 docks near the flexible active site of CAPN1, conformational changes at the interaction site after binding could be directly related to CAST4 inhibitory activity. These functional interfaces can serve as a guide to the site-mutagenesis in research on bovine CAPN1 structure-function relationships for the design of small molecules inhibitors to prevent uncontrolled and unspecific degradation in the proteolysis of key protease substrates.

Understanding the interaction determinants of CAPN1 inhibition by CAST4 from bovines using molecular modeling techniques.

HCV-induced CAPN activation and its effects on virus-infected cells in a host-immune system have been studied recently. It has been shown that the HCV-nonstructural 5A protein acts as both an inducer and a substrate for host CAPN protease; it participates in suppressing the TNF-α-induced apoptosis response and downstream IFN-induced antiviral processes. However, little is known regarding the disturbance of antiviral responses generated by bovine CAPN activation by BVDV, which is a surrogate model of HCV and is one of the most destructive diseases leading to great economic losses in cattle herds worldwide. This is also thought to be associated with the effects of either small CAPN inhibitors or the natural inhibitor CAST. They mainly bind to the binding site of CAPN substrate proteins and competitively inhibit the binding of the enzyme substrates to possibly defend against the two viruses (HCV and BVDV) for anti-viral immunity. To devise a new stratagem to discover lead candidates for an anti-BVDV drug, we first attempted to understand the bovine CAPN-CAST interaction sites and the interaction constraints of local binding architectures, were well reflected in the geometry between the pharmacophore features and its shape constraints identified using our modeled bovine CAPN1/CAST4 complex structures. We propose a computer-aided molecular design of an anti-BVDV drug as a mimetic CAST inhibitor to develop a rule-based screening function for adjusting the puzzle of relationship between bovine CAPN1 and the BVDV nonstructural proteins from all of the data obtained in the study.

Transcriptome profiles of organ tissues from pigs experimentally infected with African swine fever virus in early phase of infection.

African swine fever, caused by African swine fever virus (ASFV), is a highly contagious and fatal disease that poses a significant threat to the global pig industry. The limited information on ASFV pathogenesis and ASFV-host interactions has recently prompted numerous transcriptomic studies. However, most of these studies have focused on elucidating the transcriptome profiles of ASFV-infected porcine alveolar macrophages in vitro. Here, we analyzed dynamic transcriptional patterns in vivo in nine organ tissues (spleen, submandibular lymph node, mesenteric lymph node, inguinal lymph node, tonsils, lungs, liver, kidneys, and heart) obtained from pigs in the early stages of ASFV infection (1 and 3 d after viremia). We observed rapid spread of ASFV to the spleen after viremia, followed by broad transmission to the liver and lungs and subsequently, the submandibular and inguinal lymph nodes. Profound variations in gene expression patterns were observed across all organs and at all time-points, providing an understanding of the distinct defence strategies employed by each organ against ASFV infection. All ASFV-infected organs exhibited a collaborative response, activating immune-associated genes such as S100A8, thereby triggering a pro-inflammatory cytokine storm and interferon activation. Functional analysis suggested that ASFV exploits the PI3K-Akt signalling pathway to evade the host immune system. Overall, our findings provide leads into the mechanisms underlying pathogenesis and host immune responses in different organs during the early stages of infection, which can guide further explorations, aid the development of efficacious antiviral strategies against ASFV, and identify valuable candidate gene targets for vaccine development.

Whole-genome resequencing of Hanwoo (Korean cattle) and insight into regions of homozygosity.

BACKGROUND: Hanwoo (Korean cattle), which originated from natural crossbreeding between taurine and zebu cattle, migrated to the Korean peninsula through North China. Hanwoo were raised as draft animals until the 1970s without the introduction of foreign germplasm. Since 1979, Hanwoo has been bred as beef cattle. Genetic variation was analyzed by whole-genome deep resequencing of a Hanwoo bull. The Hanwoo genome was compared to that of two other breeds, Black Angus and Holstein, and genes within regions of homozygosity were investigated to elucidate the genetic and genomic characteristics of Hanwoo. RESULTS: The Hanwoo bull genome was sequenced to 45.6-fold coverage using the ABI SOLiD system. In total, 4.7 million single-nucleotide polymorphisms and 0.4 million small indels were identified by comparison with the Btau4.0 reference assembly. Of the total number of SNPs and indels, 58% and 87%, respectively, were novel. The overall genotype concordance between the SNPs and BovineSNP50 BeadChip data was 96.4%. Of 1.6 million genetic differences in Hanwoo, approximately 25,000 non-synonymous SNPs, splice-site variants, and coding indels (NS/SS/Is) were detected in 8,360 genes. Among 1,045 genes containing reliable specific NS/SS/Is in Hanwoo, 109 genes contained more than one novel damaging NS/SS/I. Of the genes containing NS/SS/Is, 610 genes were assigned as trait-associated genes. Moreover, 16, 78, and 51 regions of homozygosity (ROHs) were detected in Hanwoo, Black Angus, and Holstein, respectively. 'Regulation of actin filament length' was revealed as a significant gene ontology term and 25 trait-associated genes for meat quality and disease resistance were found in 753 genes that resided in the ROHs of Hanwoo. In Hanwoo, 43 genes were located in common ROHs between whole-genome resequencing and SNP chips in BTA2, 10, and 13 coincided with quantitative trait loci for meat fat traits. In addition, the common ROHs in BTA2 and 16 were in agreement between Hanwoo and Black Angus. CONCLUSIONS: We identified 4.7 million SNPs and 0.4 million small indels by whole-genome resequencing of a Hanwoo bull. Approximately 25,000 non-synonymous SNPs, splice-site variants, and coding indels (NS/SS/Is) were detected in 8,360 genes. Additionally, we found 25 trait-associated genes for meat quality and disease resistance among 753 genes that resided in the ROHs of Hanwoo. These findings will provide useful genomic information for identifying genes or casual mutations associated with economically important traits in cattle.

Gene Co-expression Analysis to Characterize Genes Related to Marbling Trait in Hanwoo (Korean) Cattle.

Marbling (intramuscular fat) is an important trait that affects meat quality and is a casual factor determining the price of beef in the Korean beef market. It is a complex trait and has many biological pathways related to muscle and fat. There is a need to identify functional modules or genes related to marbling traits and investigate their relationships through a weighted gene co-expression network analysis based on the system level. Therefore, we investigated the co-expression relationships of genes related to the 'marbling score' trait and systemically analyzed the network topology in Hanwoo (Korean cattle). As a result, we determined 3 modules (gene groups) that showed statistically significant results for marbling score. In particular, one module (denoted as red) has a statistically significant result for marbling score (p = 0.008) and intramuscular fat (p = 0.02) and water capacity (p = 0.006). From functional enrichment and relationship analysis of the red module, the pathway hub genes (IL6, CHRNE, RB1, INHBA and NPPA) have a direct interaction relationship and share the biological functions related to fat or muscle, such as adipogenesis or muscle growth. This is the first gene network study with m.logissimus in Hanwoo to observe co-expression patterns in divergent marbling phenotypes. It may provide insights into the functional mechanisms of the marbling trait.

Identification of recently selected mutations driven by artificial selection in hanwoo (korean cattle).

Hanwoo have been subjected over the last seventy years to intensive artificial selection with the aim of improving meat production traits such as marbling and carcass weight. In this study, we performed a signature of selection analysis to identify recent positive selected regions driven by a long-term artificial selection process called a breeding program using whole genome SNP data. In order to investigate homozygous regions across the genome, we estimated iES (integrated Extended Haplotype Homozygosity SNP) for the each SNPs. As a result, we identified two highly homozygous regions that seem to be strong and/or recent positive selection. Five genes (DPH5, OLFM3, S1PR1, LRRN1 and CRBN) were included in this region. To go further in the interpretation of the observed signatures of selection, we subsequently concentrated on the annotation of differentiated genes defined according to the iES value of SNPs localized close or within them. We also described the detection of the adaptive evolution at the molecular level for the genes of interest. As a result, this analysis also led to the identification of OLFM3 as having a strong signal of selection in bovine lineage. The results of this study indicate that artificial selection which might have targeted most of these genes was mainly oriented towards improvement of meat production.