Construction of a circRNA- lincRNA-lncRNA-miRNA-mRNA ceRNA regulatory network identifies genes and pathways linked to goat fertility.

Background: There is growing interest in the genetic improvement of fertility traits in female goats. With high-throughput genotyping, single-cell RNA sequencing (scRNA-seq) is a powerful tool for measuring gene expression profiles. The primary objective was to investigate comparative transcriptome profiling of granulosa cells (GCs) of high- and low-fertility goats, using scRNA-seq. Methods: Thirty samples from Ji'ning Gray goats (n = 15 for high fertility and n = 15 for low fertility) were retrieved from publicly available scRNA-seq data. Functional enrichment analysis and a literature mining approach were applied to explore modules and hub genes related to fertility. Then, interactions between types of RNAs identified were predicted, and the ceRNA regulatory network was constructed by integrating these interactions with other gene regulatory networks (GRNs). Results and discussion: Comparative transcriptomics-related analyses identified 150 differentially expressed genes (DEGs) between high- and low-fertility groups, based on the fold change (≥5 and ≤-5) and false discovery rate (FDR <0.05). Among these genes, 80 were upregulated and 70 were downregulated. In addition, 81 mRNAs, 58 circRNAs, 8 lincRNAs, 19 lncRNAs, and 55 miRNAs were identified by literature mining. Furthermore, we identified 18 hub genes (SMAD1, SMAD2, SMAD3, SMAD4, TIMP1, ERBB2, BMP15, TGFB1, MAPK3, CTNNB1, BMPR2, AMHR2, TGFBR2, BMP4, ESR1, BMPR1B, AR, and TGFB2) involved in goat fertility. Identified biological networks and modules were mainly associated with ovary signature pathways. In addition, KEGG enrichment analysis identified regulating pluripotency of stem cells, cytokine-cytokine receptor interactions, ovarian steroidogenesis, oocyte meiosis, progesterone-mediated oocyte maturation, parathyroid and growth hormone synthesis, cortisol synthesis and secretion, and signaling pathways for prolactin, TGF-beta, Hippo, MAPK, PI3K-Akt, and FoxO. Functional annotation of identified DEGs implicated important biological pathways. These findings provided insights into the genetic basis of fertility in female goats and are an impetus to elucidate molecular ceRNA regulatory networks and functions of DEGs underlying ovarian follicular development.

Massilia agrisoli sp. nov., isolated from rhizospheric soil of banana.

A Gram-stain-negative, aerobic, short rod-shaped and motile novel bacterial strain, designated MAHUQ-52T, was isolated from the rhizospheric soil of a banana plant. Colonies grew at 10-35 °C (optimum, 28 °C), pH 6.0-9.5 (optimum, pH 7.0-7.5), and in the presence of 0-1.0 % NaCl (optimum 0 %). The strain was positive for catalase and oxidase tests, as well as hydrolysis of gelatin, casein, starch and Tween 20. Based on the results of phylogenetic analysis using 16S rRNA gene and genome sequences, strain MAHUQ-52T clustered together within the genus Massilia. Strain MAHUQ-52T was closely related to Massilia soli R798T (98.6 %) and Massilia polaris RP-1-19T (98.3 %). The novel strain MAHUQ-52T has a draft genome size of 4 677 454 bp (25 contigs), annotated with 4193 protein-coding genes, 64 tRNA and 19 rRNA genes. The genomic DNA G+C content was 63.0 %. The average nucleotide identity (ANI) and digital DNA-DNA hybridization (dDDH) values between strain MAHUQ-52T and closely related type strains were ≤88.4 and 35.8 %, respectively. The only respiratory quinone was ubiquinone-8. The major fatty acids were identified as C16 : 0 and summed feature 3 (C15 : 0 iso 2-OH and/or C16 : 1 ω7c). Strain MAHUQ-52T contained phosphatidylethanolamine, diphosphatidylglycerol and phosphatidylglycerol as the major polar lipids. On the basis of dDDH and ANI values, as well as genotypic, chemotaxonomic and physiological data, strain MAHUQ-52T represents a novel species within the genus Massilia, for which the name Massilia agrisoli sp. nov. is proposed, with MAHUQ-52T (=KACC 21999T=CGMCC 1.18577T) as the type strain.

Exploring and Identifying Candidate Genes and Genomic Regions Related to Economically Important Traits in Hanwoo Cattle.

The purpose of the current review was to explore and summarize different studies concerning the detection and characterization of candidate genes and genomic regions associated with economically important traits in Hanwoo beef cattle. Hanwoo cattle, the indigenous premium beef cattle of Korea, were introduced for their marbled fat, tenderness, characteristic flavor, and juiciness. To date, there has been a strong emphasis on the genetic improvement of meat quality and yields, such as backfat thickness (BFT), marbling score (MS), carcass weight (CW), eye muscle area (EMA), and yearling weight (YW), as major selection criteria in Hanwoo breeding programs. Hence, an understanding of the genetics controlling these traits along with precise knowledge of the biological mechanisms underlying the traits would increase the ability of the industry to improve cattle to better meet consumer demands. With the development of high-throughput genotyping, genomewide association studies (GWAS) have allowed the detection of chromosomal regions and candidate genes linked to phenotypes of interest. This is an effective and useful tool for accelerating the efficiency of animal breeding and selection. The GWAS results obtained from the literature review showed that most positional genes associated with carcass and growth traits in Hanwoo are located on chromosomes 6 and 14, among which LCORL, NCAPG, PPARGC1A, ABCG2, FAM110B, FABP4, DGAT1, PLAG1, and TOX are well known. In conclusion, this review study attempted to provide comprehensive information on the identified candidate genes associated with the studied traits and genes enriched in the functional terms and pathways that could serve as a valuable resource for future research in Hanwoo breeding programs.

Multi-Omics Integration and Network Analysis Reveal Potential Hub Genes and Genetic Mechanisms Regulating Bovine Mastitis.

Mastitis, inflammation of the mammary gland, is the most prevalent disease in dairy cattle that has a potential impact on profitability and animal welfare. Specifically designed multi-omics studies can be used to prioritize candidate genes and identify biomarkers and the molecular mechanisms underlying mastitis in dairy cattle. Hence, the present study aimed to explore the genetic basis of bovine mastitis by integrating microarray and RNA-Seq data containing healthy and mastitic samples in comparative transcriptome analysis with the results of published genome-wide association studies (GWAS) using a literature mining approach. The integration of different information sources resulted in the identification of 33 common and relevant genes associated with bovine mastitis. Among these, seven genes-CXCR1, HCK, IL1RN, MMP9, S100A9, GRO1, and SOCS3-were identified as the hub genes (highly connected genes) for mastitis susceptibility and resistance, and were subjected to protein-protein interaction (PPI) network and gene regulatory network construction. Gene ontology annotation and enrichment analysis revealed 23, 7, and 4 GO terms related to mastitis in the biological process, molecular function, and cellular component categories, respectively. Moreover, the main metabolic-signalling pathways responsible for the regulation of immune or inflammatory responses were significantly enriched in cytokine-cytokine-receptor interaction, the IL-17 signaling pathway, viral protein interaction with cytokines and cytokine receptors, and the chemokine signaling pathway. Consequently, the identification of these genes, pathways, and their respective functions could contribute to a better understanding of the genetics and mechanisms regulating mastitis and can be considered a starting point for future studies on bovine mastitis.

Pinibacter aurantiacus gen. nov., sp. nov., isolated from rhizospheric soil of a pine tree.

A Gram-stain-negative, aerobic and rod-shaped novel bacterial strain, designated MAH-26T, was isolated from rhizospheric soil of a pine tree. The colonies were orange coloured, smooth, spherical and 0.7-1.8 mm in diameter when grown on Reasoner's 2A (R2A) agar for 2 days. Strain MAH-26T was able to grow at 10-40 °C, at pH 6.0-9.0 and with 0-1.0 % NaCl. Cell growth occurred on nutrient agar, R2A agar, tryptone soya agar and Luria-Bertani agar. The strain gave positive results in oxidase and catalase tests. Strain MAH-26T was closely related to Flavihumibacter sediminis CJ663T and Parasegetibacter terrae SGM2-10T with a low 16S rRNA gene sequence similarity (92.8 and 92.9 %, respectively) and phylogenetic analysis indicated that the strain formed a distinct phylogenetic lineage from the members of the closely related genera of the family Chitinophagaceae. Strain MAH-26T has a draft genome size of 6 857 405 bp, annotated with 5173 protein-coding genes, 50 tRNA and two rRNA genes. The genomic DNA G+C content was 41.5 mol%. The predominant isoprenoid quinone was menaquinone 7. The major fatty acids were identified as iso-C15:0, iso-C15:1 G and iso-C17:0 3OH. On the basis of phylogenetic inference and phenotypic, chemotaxonomic and molecular properties, strain MAH-26T represents a novel species of a novel genus of the family Chitinophagaceae, for which the name Pinibacter aurantiacus gen. nov., sp. nov. is proposed. The type strain of Pinibacter aurantiacus is MAH-26T (=KACC 19749T=CGMCC 1.13701T).

Evaluation of Genome-Enabled Prediction for Carcass Primal Cut Yields Using Single-Step Genomic Best Linear Unbiased Prediction in Hanwoo Cattle.

There is a growing interest worldwide in genetically selecting high-value cut carcass weights, which allows for increased profitability in the beef cattle industry. Primal cut yields have been proposed as a potential indicator of cutability and overall carcass merit, and it is worthwhile to assess the prediction accuracies of genomic selection for these traits. This study was performed to compare the prediction accuracy obtained from a conventional pedigree-based BLUP (PBLUP) and a single-step genomic BLUP (ssGBLUP) method for 10 primal cut traits-bottom round, brisket, chuck, flank, rib, shank, sirloin, striploin, tenderloin, and top round-in Hanwoo cattle with the estimators of the linear regression method. The dataset comprised 3467 phenotypic observations for the studied traits and 3745 genotyped individuals with 43,987 single-nucleotide polymorphisms. In the partial dataset, the accuracies ranged from 0.22 to 0.30 and from 0.37 to 0.54 as evaluated using the PBLUP and ssGBLUP models, respectively. The accuracies of PBLUP and ssGBLUP with the whole dataset varied from 0.45 to 0.75 (average 0.62) and from 0.52 to 0.83 (average 0.71), respectively. The results demonstrate that ssGBLUP performed better than PBLUP averaged over the 10 traits, in terms of prediction accuracy, regardless of considering a partial or whole dataset. Moreover, ssGBLUP generally showed less biased prediction and a value of dispersion closer to 1 than PBLUP across the studied traits. Thus, the ssGBLUP seems to be more suitable for improving the accuracy of predictions for primal cut yields, which can be considered a starting point in future genomic evaluation for these traits in Hanwoo breeding practice.

Genome-wide Association Study for Carcass Primal Cut Yields Using Single-step Bayesian Approach in Hanwoo Cattle.

The importance of meat and carcass quality is growing in beef cattle production to meet both producer and consumer demands. Primal cut yields, which reflect the body compositions of carcass, could determine the carcass grade and, consequently, command premium prices. Despite its importance, there have been few genome-wide association studies on these traits. This study aimed to identify genomic regions and putative candidate genes related to 10 primal cut traits, including tenderloin, sirloin, striploin, chuck, brisket, top round, bottom round, shank, flank, and rib in Hanwoo cattle using a single-step Bayesian regression (ssBR) approach. After genomic data quality control, 43,987 SNPs from 3,745 genotyped animals were available, of which 3,467 had phenotypic records for the analyzed traits. A total of 16 significant genomic regions (1-Mb window) were identified, of which five large-effect quantitative trait loci (QTLs) located on chromosomes 6 at 38-39 Mb, 11 at 21-22 Mb, 14 at 6-7 Mb and 26-27 Mb, and 19 at 26-27 Mb were associated with more than one trait, while the remaining 11 QTLs were trait-specific. These significant regions were harbored by 154 genes, among which TOX, FAM184B, SPP1, IBSP, PKD2, SDCBP, PIGY, LCORL, NCAPG, and ABCG2 were noteworthy. Enrichment analysis revealed biological processes and functional terms involved in growth and lipid metabolism, such as growth (GO:0040007), muscle structure development (GO:0061061), skeletal system development (GO:0001501), animal organ development (GO:0048513), lipid metabolic process (GO:0006629), response to lipid (GO:0033993), metabolic pathways (bta01100), focal adhesion (bta04510), ECM-receptor interaction (bta04512), fat digestion and absorption (bta04975), and Rap1 signaling pathway (bta04015) being the most significant for the carcass primal cut traits. Thus, identification of quantitative trait loci regions and plausible candidate genes will aid in a better understanding of the genetic and biological mechanisms regulating carcass primal cut yields.

Estimation of Genetic Correlations of Primal Cut Yields with Carcass Traits in Hanwoo Beef Cattle.

This study was carried out to estimate the variance components, heritability, and genetic correlations between the carcass traits and primal cut yields in Hanwoo cattle. Carcass traits comprising 5622 records included back fat thickness (BFT), carcass weight (CW), eye muscle area (EMA), and marbling score (MS). The 10 primal cut yields from 3467 Hanwoo steers included the tenderloin (TLN), sirloin (SLN), striploin (STLN), chuck (CHK), brisket (BSK), top round (TRD), bottom round (BRD), rib (RB), shank (SK), and flank (FK). In addition, three composite traits were formed by combining primal cut yields as novel traits according to consumer preferences and market price: high-value cuts (HVC), medium-value cuts (MVC), and low-value cuts (LVC). Heritability estimates for the interest of traits were moderate to high, ranging from 0.21 ± 0.04 for CHK to 0.59 ± 0.05 for MS. Except genetic correlations between RB and other primal cut traits, favorable and moderate to high correlations were observed among the yields of primal cut that ranged from 0.38 ± 0.14 (CHK and FK) to 0.93 ± 0.01 (TRD and BRD). Moreover, the estimated genetic correlations of CW and EMA with primal cut yields and three composite traits were positive and moderate to strong, except for BFT, which was negative. These results indicate that genetic progress can be achieved for all traits, and selection to increase the yields of primal cuts can lead to considerable profitability in the Hanwoo beef industry.

Multi-Trait Single-Step GBLUP Improves Accuracy of Genomic Prediction for Carcass Traits Using Yearling Weight and Ultrasound Traits in Hanwoo.

There has been a growing interest in the genetic improvement of carcass traits as an important and primary breeding goal in the beef cattle industry over the last few decades. The use of correlated traits and molecular information can aid in obtaining more accurate estimates of breeding values. This study aimed to assess the improvement in the accuracy of genetic predictions for carcass traits by using ultrasound measurements and yearling weight along with genomic information in Hanwoo beef cattle by comparing four evaluation models using the estimators of the recently developed linear regression method. We compared the performance of single-trait pedigree best linear unbiased prediction [ST-BLUP and single-step genomic (ST-ssGBLUP)], as well as multi-trait (MT-BLUP and MT-ssGBLUP) models for the studied traits at birth and yearling date of steers. The data comprised of 15,796 phenotypic records for yearling weight and ultrasound traits as well as 5,622 records for carcass traits (backfat thickness, carcass weight, eye muscle area, and marbling score), resulting in 43,949 single-nucleotide polymorphisms from 4,284 steers and 2,332 bulls. Our results indicated that averaged across all traits, the accuracy of ssGBLUP models (0.52) was higher than that of pedigree-based BLUP (0.34), regardless of the use of single- or multi-trait models. On average, the accuracy of prediction can be further improved by implementing yearling weight and ultrasound data in the MT-ssGBLUP model (0.56) for the corresponding carcass traits compared to the ST-ssGBLUP model (0.49). Moreover, this study has shown the impact of genomic information and correlated traits on predictions at the yearling date (0.61) using MT-ssGBLUP models, which was advantageous compared to predictions at birth date (0.51) in terms of accuracy. Thus, using genomic information and high genetically correlated traits in the multi-trait model is a promising approach for practical genomic selection in Hanwoo cattle, especially for traits that are difficult to measure.

Estimation of Genetic Parameters and Correlation between Yearling Ultrasound Measurements and Carcass Traits in Hanwoo Cattle.

Genetic parameters have a significant role in designing a breeding program and are required to evaluate economically important traits. The objective of this study was to estimate heritability and genetic correlation between yearling ultrasound measurements, such as backfat thickness (UBFT), eye muscle area (UEMA), intramuscular fat content (UIMF), and carcass traits, such as backfat thickness (BFT), carcass weight (CW), eye muscle area (EMA), marbling score (MS) at approximately 24 months of age, as well as yearling weight (YW) in Hanwoo bulls (15,796) and steers (5682). The (co) variance components were estimated using a multi-trait animal model. Moderate to high heritability estimates were obtained and were 0.42, 0.50, 0.56, and 0.59 for CW, EMA, BFT, and MS, respectively. Heritability estimates for yearling measurements of YW, UEMA, UBFT, and UIMF were 0.31, 0.32, 0.30, and 0.19, respectively. Favorable and strong genetic correlations were observed between UIMF and MS (0.78), UBFT and BFT (0.63), and UEMA and EMA (0.65). Moreover, the estimated genetic correlation between YW and CW was high (0.84) and relatively moderate between YW and EMA (0.43). These results suggest that genetic improvement can be achieved for carcass traits when using yearling ultrasound measurements as selection criteria in ongoing Hanwoo breeding programs.

Genomic Prediction Using Alternative Strategies of Weighted Single-Step Genomic BLUP for Yearling Weight and Carcass Traits in Hanwoo Beef Cattle.

The weighted single-step genomic best linear unbiased prediction (GBLUP) method has been proposed to exploit information from genotyped and non-genotyped relatives, allowing the use of weights for single-nucleotide polymorphism in the construction of the genomic relationship matrix. The purpose of this study was to investigate the accuracy of genetic prediction using the following single-trait best linear unbiased prediction methods in Hanwoo beef cattle: pedigree-based (PBLUP), un-weighted (ssGBLUP), and weighted (WssGBLUP) single-step genomic methods. We also assessed the impact of alternative single and window weighting methods according to their effects on the traits of interest. The data was comprised of 15,796 phenotypic records for yearling weight (YW) and 5622 records for carcass traits (backfat thickness: BFT, carcass weight: CW, eye muscle area: EMA, and marbling score: MS). Also, the genotypic data included 6616 animals for YW and 5134 for carcass traits on the 43,950 single-nucleotide polymorphisms. The ssGBLUP showed significant improvement in genomic prediction accuracy for carcass traits (71%) and yearling weight (99%) compared to the pedigree-based method. The window weighting procedures performed better than single SNP weighting for CW (11%), EMA (11%), MS (3%), and YW (6%), whereas no gain in accuracy was observed for BFT. Besides, the improvement in accuracy between window WssGBLUP and the un-weighted method was low for BFT and MS, while for CW, EMA, and YW resulted in a gain of 22%, 15%, and 20%, respectively, which indicates the presence of relevant quantitative trait loci for these traits. These findings indicate that WssGBLUP is an appropriate method for traits with a large quantitative trait loci effect.

Genetic parameters and correlations of related feed efficiency, growth, and carcass traits in Hanwoo beef cattle.

OBJECTIVE: This study aimed to estimate the genetic parameters and genetic correlations for related feed efficiency, growth, and carcass traits in Hanwoo cattle. METHODS: Phenotypic data from 15,279 animals born between 1989 and 2015 were considered. The related feed efficiency traits considered were Kleiber ratio (KR) and relative growth rate (RGR). Carcass traits analyzed were backfat thickness (BT), carcass weight, eye muscle area, and marbling score. Growth traits were assessed by the average daily gain (ADG), metabolic body weight (MBW) at mid-test age from 6 to 24 months, and yearling weight (YW). Variance and covariance components were estimated using restricted maximum likelihood using nine multi-trait animal models. RESULTS: The heritability estimates for related feed efficiency (0.28±0.04 for KR and RGR) and growth traits (0.26±0.02 to 0.33±0.04) were moderate, but the carcass traits tended to be higher (0.38±0.04 to 0.61±0.06). The related feed efficiency traits were positively genetically correlated with all the carcass traits (0.37±0.09 to 0.47±0.07 for KR, and 0.14±0.09 to 0.37±0.09 for RGR), except for BT, which showed null to weak correlation. Conversely, the genetic correlations of RGR with MBW (-0.36±0.08) and YW (-0.30±0.08) were negative, and those of KR with MBW and YW were close to zero, whereas the genetic correlations of ADG with RGR (0.40±0.08) and KR (0.70±0.05) were positive and relatively moderate to high. The genetic (0.92±0.02) correlations between KR and RGR were very high. CONCLUSION: Sufficient genetic variability and heritability were observed for traits of interest. Moreover, the inclusion of KR and/or RGR in Hanwoo cattle breeding programs could improve the feed efficiency without producing any unfavorable effects on the carcass traits.

Improving the accuracy of genomic evaluation for linear body measurement traits using single-step genomic best linear unbiased prediction in Hanwoo beef cattle.

BACKGROUND: Recently, there has been a growing interest in the genetic improvement of body measurement traits in farm animals. They are widely used as predictors of performance, longevity, and production traits, and it is worthwhile to investigate the prediction accuracies of genomic selection for these traits. In genomic prediction, the single-step genomic best linear unbiased prediction (ssGBLUP) method allows the inclusion of information from genotyped and non-genotyped relatives in the analysis. Hence, we aimed to compare the prediction accuracy obtained from a pedigree-based BLUP only on genotyped animals (PBLUP-G), a traditional pedigree-based BLUP (PBLUP), a genomic BLUP (GBLUP), and a single-step genomic BLUP (ssGBLUP) method for the following 10 body measurement traits at yearling age of Hanwoo cattle: body height (BH), body length (BL), chest depth (CD), chest girth (CG), chest width (CW), hip height (HH), hip width (HW), rump length (RL), rump width (RW), and thurl width (TW). The data set comprised 13,067 phenotypic records for body measurement traits and 1523 genotyped animals with 34,460 single-nucleotide polymorphisms. The accuracy for each trait and model was estimated only for genotyped animals using five-fold cross-validations. RESULTS: The accuracies ranged from 0.02 to 0.19, 0.22 to 0.42, 0.21 to 0.44, and from 0.36 to 0.55 as assessed using the PBLUP-G, PBLUP, GBLUP, and ssGBLUP methods, respectively. The average predictive accuracies across traits were 0.13 for PBLUP-G, 0.34 for PBLUP, 0.33 for GBLUP, and 0.45 for ssGBLUP methods. Our results demonstrated that averaged across all traits, ssGBLUP outperformed PBLUP and GBLUP by 33 and 43%, respectively, in terms of prediction accuracy. Moreover, the least root of mean square error was obtained by ssGBLUP method. CONCLUSIONS: Our findings suggest that considering the ssGBLUP model may be a promising way to ensure acceptable accuracy of predictions for body measurement traits, especially for improving the prediction accuracy of selection candidates in ongoing Hanwoo breeding programs.

Integrating Single-Step GWAS and Bipartite Networks Reconstruction Provides Novel Insights into Yearling Weight and Carcass Traits in Hanwoo Beef Cattle.

In recent years, studies on the biological mechanisms underlying complex traits have been facilitated by innovations in high-throughput genotyping technology. We conducted a weighted single-step genome-wide association study (WssGWAS) to evaluate backfat thickness, carcass weight, eye muscle area, marbling score, and yearling weight in a cohort of 1540 Hanwoo beef cattle using BovineSNP50 BeadChip. The WssGWAS uncovered thirty-three genomic regions that explained more than 1% of the additive genetic variance, mostly located on chromosomes 6 and 14. Among the identified window regions, seven quantitative trait loci (QTL) had pleiotropic effects and twenty-six QTL were trait-specific. Significant pathways implicated in the measured traits through Gene Ontology (GO) term enrichment analysis included the following: lipid biosynthetic process, regulation of lipid metabolic process, transport or localization of lipid, regulation of growth, developmental growth, and multicellular organism growth. Integration of GWAS results of the studied traits with pathway and network analyses facilitated the exploration of the respective candidate genes involved in several biological functions, particularly lipid and growth metabolism. This study provides novel insight into the genetic bases underlying complex traits and could be useful in developing breeding schemes aimed at improving growth and carcass traits in Hanwoo beef cattle.

Comparison of conventional BLUP and single-step genomic BLUP evaluations for yearling weight and carcass traits in Hanwoo beef cattle using single trait and multi-trait models.

Hanwoo, an important indigenous and popular breed of beef cattle in Korea, shows rapid growth and has high meat quality. Its yearling weight (YW) and carcass traits (backfat thickness, carcass weight- CW, eye muscle area, and marbling score) are economically important for selection of young and proven bulls. However, measuring carcass traits is difficult and expensive, and can only be performed postmortem. Genomic selection has become an appealing procedure for genetic evaluation of these traits (by inclusion of the genomic data) along with the possibility of multi-trait analysis. The aim of this study was to compare conventional best linear unbiased prediction (BLUP) and single-step genomic BLUP (ssGBLUP) methods, using both single-trait (ST-BLUP, ST-ssGBLUP) and multi-trait (MT-BLUP, MT-ssGBLUP) models to investigate the improvement of breeding-value accuracy for carcass traits and YW. The data comprised of 15,279 phenotypic records for YW and 5,824 records for carcass traits, and 1,541 genotyped animals for 34,479 single-nucleotide polymorphisms. Accuracy for each trait and model was estimated only for genotyped animals by five-fold cross-validation. ssGBLUP models (ST-ssGBLUP and MT-ssGBLUP) showed ~19% and ~36% greater accuracy than conventional BLUP models (ST-BLUP and MT-BLUP) for YW and carcass traits, respectively. Within ssGBLUP models, the accuracy of the genomically estimated breeding value for CW increased (19%) when ST-ssGBLUP was replaced with the MT-ssGBLUP model, as the inclusion of YW in the analysis led to a strong genetic correlation with CW (0.76). For backfat thickness, eye muscle area, and marbling score, ST- and MT-ssGBLUP models yielded similar accuracy. Thus, combining pedigree and genomic data via the ssGBLUP model may be a promising way to ensure acceptable accuracy of predictions, especially among young animals, for ongoing Hanwoo cattle breeding programs. MT-ssGBLUP is highly recommended when phenotypic records are limited for one of the two highly correlated genetic traits.

Exploring novel single nucleotide polymorphisms and haplotypes of the diacylglycerol O-acyltransferase 1 (DGAT1) gene and their effects on protein structure in Iranian buffalo.

BACKGROUND: Diacylglycerol O-acyltransferase 1 (DGAT1) plays a key role in the synthesis of triglycerides. Recent studies have shown that a transition mutation resulting in substitutions of guanine by adenine in the DGAT1 gene in cattle has considerable effects on milk yield and composition. Currently, there is no systematic research reporting on the utilization of this gene segment in Iranian buffalo (Bubalus bubalis). OBJECTIVE: In this study, the genetic differentiation of three indigenous Iranian buffalo populations was investigated in the region spanning exon 3 to exon 17 of the DGAT1 gene. METHODS: A total of 200 buffaloes were genotyped, all the samples were sequenced directly in both directions with forward and reverse sequencing primers. RESULTS: Sequence analysis showed novel SNPs compared to the reference GenBank sequence (DQ886485) at nucleotide positions g.6097A>G, g.7036C>T, g.7338G>A, g.7710C>T, g.8087C>T, g.8259G>A, g.8275G>A, g.8367C>T, and g.8426C>T. No polymorphisms were found within exon 8. Therefore, the K232A position was thought to be a conserved and fixed region for high milk fat content (K allele) in Bos indicus and all buffalo breeds. Comparison with Indian buffalo revealed three exonic SNPs, one of which was nonsynonymous. A unique 22 bp insertion was observed in intron 10 of DGAT1. Linkage disequilibrium analysis allowed the identification of nine haplotypes among the sampled animals. To our knowledge, this is the first report of sequencing analysis of the DGAT1 gene in Iranian buffalo. CONCLUSION: Our results suggest that genetic diversity exists and could be useful in examining the association between the DGAT1 gene and milk production traits in buffalo.

Predictive performance of genomic selection methods for carcass traits in Hanwoo beef cattle: impacts of the genetic architecture.

BACKGROUND: Hanwoo beef is known for its marbled fat, tenderness, juiciness and characteristic flavor, as well as for its low cholesterol and high omega 3 fatty acid contents. As yet, there has been no comprehensive investigation to estimate genomic selection accuracy for carcass traits in Hanwoo cattle using dense markers. This study aimed at evaluating the accuracy of alternative statistical methods that differed in assumptions about the underlying genetic model for various carcass traits: backfat thickness (BT), carcass weight (CW), eye muscle area (EMA), and marbling score (MS). METHODS: Accuracies of direct genomic breeding values (DGV) for carcass traits were estimated by applying fivefold cross-validation to a dataset including 1183 animals and approximately 34,000 single nucleotide polymorphisms (SNPs). RESULTS: Accuracies of BayesC, Bayesian LASSO (BayesL) and genomic best linear unbiased prediction (GBLUP) methods were similar for BT, EMA and MS. However, for CW, DGV accuracy was 7% higher with BayesC than with BayesL and GBLUP. The increased accuracy of BayesC, compared to GBLUP and BayesL, was maintained for CW, regardless of the training sample size, but not for BT, EMA, and MS. Genome-wide association studies detected consistent large effects for SNPs on chromosomes 6 and 14 for CW. CONCLUSIONS: The predictive performance of the models depended on the trait analyzed. For CW, the results showed a clear superiority of BayesC compared to GBLUP and BayesL. These findings indicate the importance of using a proper variable selection method for genomic selection of traits and also suggest that the genetic architecture that underlies CW differs from that of the other carcass traits analyzed. Thus, our study provides significant new insights into the carcass traits of Hanwoo cattle.

Random Regression Models Using Legendre Polynomials to Estimate Genetic Parameters for Test-day Milk Protein Yields in Iranian Holstein Dairy Cattle.

The objective of this study was to estimate the genetic parameters of milk protein yields in Iranian Holstein dairy cattle. A total of 1,112,082 test-day milk protein yield records of 167,269 first lactation Holstein cows, calved from 1990 to 2010, were analyzed. Estimates of the variance components, heritability, and genetic correlations for milk protein yields were obtained using a random regression test-day model. Milking times, herd, age of recording, year, and month of recording were included as fixed effects in the model. Additive genetic and permanent environmental random effects for the lactation curve were taken into account by applying orthogonal Legendre polynomials of the fourth order in the model. The lowest and highest additive genetic variances were estimated at the beginning and end of lactation, respectively. Permanent environmental variance was higher at both extremes. Residual variance was lowest at the middle of the lactation and contrarily, heritability increased during this period. Maximum heritability was found during the 12th lactation stage (0.213±0.007). Genetic, permanent, and phenotypic correlations among test-days decreased as the interval between consecutive test-days increased. A relatively large data set was used in this study; therefore, the estimated (co)variance components for random regression coefficients could be used for national genetic evaluation of dairy cattle in Iran.