Han-Ha Chai

Whole-Genome Resequencing Analysis of Hanwoo and Yanbian Cattle to Identify Genome-Wide SNPs and Signatures of Selection.

Over the last 30 years, Hanwoo has been selectively bred to improve economically important traits. Hanwoo is currently the representative Korean native beef cattle breed, and it is believed that it shared an ancestor with a Chinese breed, Yanbian cattle, until the last century. However, these two breeds have experienced different selection pressures during recent decades. Here, we whole-genome sequenced 10 animals each of Hanwoo and Yanbian cattle (20 total) using the Illumina HiSeq 2000 sequencer. A total of approximately 3.12 and 3.07 billion sequence reads were mapped to the bovine reference sequence assembly (UMD 3.1) at an average of approximately 10.71- and 10.53-fold coverage for Hanwoo and Yanbian cattle, respectively. A total of 17,936,399 single nucleotide polymorphisms (SNPs) were yielded, of which 22.3% were found to be novel. By annotating the SNPs, we further retrieved numerous nonsynonymous SNPs that may be associated with traits of interest in cattle. Furthermore, we performed whole-genome screening to detect signatures of selection throughout the genome. We located several promising selective sweeps that are potentially responsible for economically important traits in cattle; the PPP1R12A gene is an example of a gene that potentially affects intramuscular fat content. These discoveries provide valuable genomic information regarding potential genomic markers that could predict traits of interest for breeding programs of these cattle breeds.

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.

Characterization of Genes for Beef Marbling Based on Applying Gene Coexpression Network

Marbling is an important trait in characterization beef quality and a major factor for determining the price of beef in the Korean beef market. In particular, marbling is a complex trait and needs a system-level approach for identifying candidate genes related to the trait. To find the candidate gene associated with marbling, we used a weighted gene coexpression network analysis from the expression value of bovine genes. Hub genes were identified; they were topologically centered with large degree and BC values in the global network. We performed gene expression analysis to detect candidate genes in M. longissimus with divergent marbling phenotype (marbling scores 2 to 7) using qRT-PCR. The results demonstrate that transmembrane protein 60 (TMEM60) and dihydropyrimidine dehydrogenase (DPYD) are associated with increasing marbling fat. We suggest that the network-based approach in livestock may be an important method for analyzing the complex effects of candidate genes associated with complex traits like marbling or tenderness.

Copy number variations in Hanwoo and Yanbian cattle genomes using the massively parallel sequencing data.

Hanwoo is an indigenous Korean beef cattle breed, and it shared an ancestor with Yanbian cattle that are found in the Northeast provinces in China until the last century. During recent decades, those cattle breeds experienced different selection pressures. Here, we present genome-wide copy number variations (CNVs) by comparing Hanwoo and Yanbian cattle sequencing data. We used ~3.12 and ~3.07 billion sequence reads from Hanwoo and Yanbian cattle, respectively. A total of 901 putative CNV regions (CNVRs) were identified throughout the genome, representing 5,513,340bp. This is a smaller number than has been reported in previous studies, indicating that Hanwoo are genetically close to Yanbian cattle. Of the CNVRs, 53.2% and 46.8% were found to be gains and losses in Hanwoo. Potential functional roles of each CNVR were assessed by annotating all CNVRs and gene ontology (GO) enrichment analysis. We found that 278 CNVRs overlapped with cattle gene-sets (genic-CNVRs) that could be promising candidates to account for economically important traits in cattle. The enrichment analysis indicated that genes were significantly over-represented in GO terms, including developmental process, multicellular organismal process, reproduction, and response to stimulus. These results provide a valuable genomic resource for determining how CNVs are associated with cattle traits.

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.

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 Ca2+-binding to the enzyme is thought to include activation, as well as the dissociation, aggregation, and autolysis of small regular subunits. Unfortunately, the Ca2+-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 Ca2+-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 Ca2+-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.

Demographic Trends in Korean Native Cattle Explained Using Bovine SNP50 Beadchip

Linkage disequilibrium (LD) is the non-random association between the loci and it could give us a preliminary insight into the genetic history of the population. In the present study LD patterns and effective population size (Ne) of three Korean cattle breeds along with Chinese, Japanese and Mongolian cattle were compared using the bovine Illumina SNP50 panel. The effective population size (Ne) is the number of breeding individuals in a population and is particularly important as it determines the rate at which genetic variation is lost. The genotype data in our study comprised a total of 129 samples, varying from 4 to 39 samples. After quality control there were ~29,000 single nucleotide polymorphisms (SNPs) for which r2 value was calculated. Average distance between SNP pairs was 1.14 Mb across all breeds. Average r2 between adjacent SNP pairs ranged between was 0.1 for Yanbian to 0.3 for Qinchuan. Effective population size of the breeds based on r2 varied from 16 in Hainan to 226 in Yanbian. Amongst the Korean native breeds effective population size of Brindle Hanwoo was the least with Ne = 59 and Brown Hanwoo was the highest with Ne = 83. The effective population size of the Korean cattle breeds has been decreasing alarmingly over the past generations. We suggest appropriate measures to be taken to prevent these local breeds in their native tracts.

HSF2 Co-regulates Protein-coding and Long Non-coding RNA Genes Specific to Black Tissues of the Black Chicken, Yeonsan Ogye

The Yeonsan Ogye (Ogye) is a rare Korean domestic chicken breed, the entire body of which, including its 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 globally investigated. In this study, high-throughput RNA sequencing and DNA methylation sequencing were performed to dissect the expression landscape of 14,264 Ogye protein-coding and 6900 long non-coding RNA (lncRNA) genes along with DNA methylation landscape in twenty different Ogye tissues. About 75% of Ogye lncRNAs showed tissue-specific expression whereas about 45% of protein-coding genes did. For some genes, the tissue-specific expression levels were inversely correlated with DNA methylation levels in their promoters. About 39% of the tissue-specific lncRNAs displayed functional association with proximal or distal protein-coding genes. In particular, heat shock transcription factor 2 (HSF2)-associated lncRNAs were discovered to be functionally linked to protein-coding genes that are specifically expressed in black skin tissues, tended to be more syntenically conserved in mammals, and were differentially expressed in black tissues relative to white tissues. Our results not only facilitate understanding how the non-coding genome regulates unique phenotypes but also should be of use for future genomic breeding of chickens.The Yeonsan Ogye (Ogye) is a rare chicken breed that populates the Korean peninsula. The entire body of this bird, including its 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 globally investigated. In this study, high-throughput RNA sequencing (RNA-seq) and reduced representation bisulfite sequencing (RRBS) were performed to construct whole non-coding transcriptome maps across twenty different Ogye tissues. The resulting maps included 6900 long non-coding RNA (lncRNA) genes comprising 1290 known and 5610 novel lncRNA genes. Compared to lncRNAs previously annotated in the Galllus gallus red junglefowl, a considerable number were either fragments of protein-coding genes or not expressed in Ogye tissues. Newly annotated Ogye lncRNA genes showed tissue-specific expression and simple gene structures containing 2 or 3 exons. Systematic analyses of sequencing data and other genomic data demonstrated that about 39% of the tissue-specific lncRNAs displayed evidence of function. In particular, heat shock transcription factor 2 (HSF2)-associated lncRNAs were discovered to be functionally linked to protein-coding genes that are specifically expressed in black skin tissues, tended to be more syntenically conserved in mammals, and were differentially expressed in black tissues relative to white tissues. Our findings and resulting maps provide not only a comprehensive catalogue of lncRNAs but also a set of functional lncRNAs that will facilitate understanding how the non-coding genome regulates unique phenotypes. Furthermore, our results should be of use for future genomic breeding of chickens.

Effects of a Novel p.A41P Mutation in the Swine Myogenic factor 5 (MYF5) Gene on Protein Stabilizing, Muscle Fiber Characteristics and Meat Quality

Abstract Myogenic factor 5 (MYF5) plays an important role in regulating skeletal muscle fiber characteristics, consequently affecting meat production and quality. We identified a novel p.A41P mutation in exon1 of the porcine MYF5 gene by direct sequencing. The mutation was predicted to be destabilizing in protein structure based on the resultant amino acid substitution. We estimated the significant substitution effect of p.A41P on the energy stabilization of Myf5 protein structure. Then, we demonstrated that the mutation in Yorkshire population significantly affected muscle fiber type I composition (p<0.05), loin-eye area of lean meat content (p<0.05) and filter-fluid uptake of meat quality (p<0.01). Furthermore, dominant effects significantly influenced total muscle fiber number (p<0.05). This study suggests that the novel p.A41P mutation in porcine MYF5 may be a valuable genetic marker to affect the muscle fiber characteristics and consequently improve meat production quality and quantity.

Structural analyses of putative porcine cluster of differentiation (CD) 24 gene

24 (HSA) molecule is a mucin-type glycoprotein attached to the cell surface by glycosylphosphatidylinositol (GPI)anchor, promoting adhesive interactions between cells or in extracellular matrix. The aim of this study was to determine the not yet fully identified porcine CD24 gene and protein structure using computational analysis and to validate variants reported in exons of CD24 gene using direct sequencing. A total of 59 samples belonging to Yorkshire, Landrace, Berkshire, Jeju black pig and wild boar were used in the study. Human CD24 mRNA sequences were used as a reference and subjected to BLAST searches to retrieve the orthologous expressed sequence tags (ESTs) or cDNA sequences against NCBI and Ensemble databases. Assembled ESTs and retrieved cDNA sequences for the porcine CD24 gene were used for specific BLAST search to determine its genomic structure. We found porcine CD24 gene to consist of two exons and a relatively long intron. Second exon of porcine CD24 gene had a long 3’ untranslated region (UTR) and was very similar to that of human, mouse, rat, and sheep. The sequence homology of porcine CD24 protein was 65.38-84.62%, when analyzed with amino acid sequences of rat, mouse, human, cattle, and sheep CD24 protein. Nterminal signal sequence, O-glycosylation sites and GPI-anchoring signal sites were also predicted in pig, which showed these motifs to be evolutionary conserved across the species. Variant analysis in exonic regions of porcine CD24 among the multiple breeds showed that only second exon contained eight SNPs and three insertions in a 3’ UTR. Taken together, this study reports putative porcine CD24 gene and its protein structure using in silico approaches, which will be helpful for any further functional studies. Original Article