Discovery of novel genetic networks associated with 19 economically important traits in beef cattle

A network-based approach to understanding gene-biological processes affecting economically important traits of Nelore cattle

This study aimed to build gene-biological process networks with differentially expressed genes associated with economically important traits of Nelore cattle from 17 previous studies. The genes were clustered into three groups by evaluated traits: group 1, production traits; group 2, carcass traits; and group 3, meat quality traits. For each group, a gene-biological process network analysis was performed with the differentially expressed genes in common. For production traits, 37 genes were found in common, of which 13 genes were enriched for six Gene Ontology (GO) terms; these terms were not functionally grouped. However, the enriched GO terms were related to homeostasis, the development of muscles and the immune system. For carcass traits, four genes were found in common. Thus, it was not possible to functionally group these genes into a network. For meat quality traits, the analysis revealed 222 genes in common. CSRP3 was the only gene differentially expressed in all three groups. Non-redundant biological terms for clusters of genes were functionally grouped networks, reflecting the cross-talk between all biological processes and genes involved. Many biological processes and pathways related to muscles, the immune system and lipid metabolism were enriched, such as striated muscle cell development and triglyceride metabolic processes. This study provides insights into the genetic mechanisms of production, carcass and meat quality traits of Nelore cattle. This information is fundamental for a better understanding of the complex traits and could help in planning strategies for the production and selection systems of Nelore cattle.

Meta-analysis of RNA-Seq datasets highlights novel genes/pathways involved in fat deposition in fat-tail of sheep

Introduction

Fat-tail in sheep is considered as an important energy reservoir to provide energy as a survival buffer during harsh challenges. However, fat-tail is losing its importance in modern sheep industry systems and thin-tailed breeds are more desirable. Using comparative transcriptome analysis to compare fat-tail tissue between fat- and thin-tailed sheep breeds provides a valuable approach to study the complex genetic factors associated with fat-tail development. However, transcriptomic studies often suffer from issues with reproducibility, which can be improved by integrating multiple studies based on a meta-analysis.

Methods

Hence, for the first time, an RNA-Seq meta-analysis on sheep fat-tail transcriptomes was performed using six publicly available datasets.

Results and discussion

A total of 500 genes (221 up-regulated, 279 down-regulated) were identified as differentially expressed genes (DEGs). A jackknife sensitivity analysis confirmed the robustness of the DEGs. Moreover, QTL and functional enrichment analysis reinforced the importance of the DEGs in the underlying molecular mechanisms of fat deposition. Protein-protein interactions (PPIs) network analysis revealed the functional interactions among the DEGs and the subsequent sub-network analysis led to identify six functional sub-networks. According to the results of the network analysis, down-regulated DEGs in green and pink sub-networks (like collagen subunits IV, V, and VI, integrins 1 and 2, SCD, SCD5, ELOVL6, ACLY, SLC27A2, and LPIN1) may impair lipolysis or fatty acid oxidation and cause fat accumulation in tail. On the other hand, up-regulated DEGs, especially those are presented in green and pink sub-networks (like IL6, RBP4, LEPR, PAI-1, EPHX1, HSD11B1, and FMO2), might contribute to a network controlling fat accumulation in the tail of sheep breed through mediating adipogenesis and fatty acid biosynthesis. Our results highlighted a set of known and novel genes/pathways associated with fat-tail development, which could improve the understanding of molecular mechanisms behind fat deposition in sheep fat-tail.

Multi-locus genome-wide association studies reveal genomic regions and putative candidate genes associated with leaf spot diseases in African groundnut (Arachis hypogaea L.) germplasm

Early leaf spot (ELS) and late leaf spot (LLS) diseases are the two most destructive groundnut diseases in Ghana resulting in ≤ 70% yield losses which is controlled largely by chemical method. To develop leaf spot resistant varieties, the present study was undertaken to identify single nucleotide polymorphism (SNP) markers and putative candidate genes underlying both ELS and LLS. In this study, six multi-locus models of genome-wide association study were conducted with the best linear unbiased predictor obtained from 294 African groundnut germplasm screened for ELS and LLS as well as image-based indices of leaf spot diseases severity in 2020 and 2021 and 8,772 high-quality SNPs from a 48 K SNP array Axiom platform. Ninety-seven SNPs associated with ELS, LLS and five image-based indices across the chromosomes in the 2 two sub-genomes. From these, twenty-nine unique SNPs were detected by at least two models for one or more traits across 16 chromosomes with explained phenotypic variation ranging from 0.01 - 62.76%, with exception of chromosome (Chr) 08 (Chr08), Chr10, Chr11, and Chr19. Seventeen potential candidate genes were predicted at ± 300 kbp of the stable/prominent SNP positions (12 and 5, down- and upstream, respectively). The results from this study provide a basis for understanding the genetic architecture of ELS and LLS diseases in African groundnut germplasm, and the associated SNPs and predicted candidate genes would be valuable for breeding leaf spot diseases resistant varieties upon further validation.

Genome-Wide Association Analysis Reveals Novel Loci Related with Visual Score Traits in Nellore Cattle Raised in Pasture-Based Systems

Body conformation traits assessed based on visual scores are widely used in Zebu cattle breeding programs. The aim of this study was to identify genomic regions and biological pathways associated with body conformation (CONF), finishing precocity (PREC), and muscling (MUSC) in Nellore cattle. The measurements based on visual scores were collected in 20,807 animals raised in pasture-based systems in Brazil. In addition, 2775 animals were genotyped using a 35 K SNP chip, which contained 31,737 single nucleotide polymorphisms after quality control. Single-step GWAS was performed using the BLUPF90 software while candidate genes were identified based on the Ensembl Genes 69. PANTHER and REVIGO platforms were used to identify key biological pathways and STRING to create gene networks. Novel candidate genes were revealed associated with CONF, including ALDH9A1, RXRG, RAB2A, and CYP7A1, involved in lipid metabolism. The genes associated with PREC were ELOVL5, PID1, DNER, TRIP12, and PLCB4, which are related to the synthesis of long-chain fatty acids, lipid metabolism, and muscle differentiation. For MUSC, the most important genes associated with muscle development were SEMA6A, TIAM2, UNC5A, and UIMC1. The polymorphisms identified in this study can be incorporated in commercial genotyping panels to improve the accuracy of genomic evaluations for visual scores in beef cattle.

Effect of genetic variants in the SMAD1 and SMAD5 genes promoter on growth and beef quality traits in cattle

The SMAD1 and SMAD5 genes belong to mothers against decapentaplegic proteins family, which participate in the BMP pathway to control skeletal myogenesis and growth. In the present study, we analyzed the associations between polymorphisms of SMAD1 and SMAD5 genes promoter and important economical traits in Qinchuan cattle. Four SNPs in the SMAD1 gene promoter and three SNPs in the SMAD5 promoter were identified by sequencing of 448 Qinchuan cattles. Allelic and frequency analyses of these SNPs resulted in eight haplotypes both in the promoters of the two genes promoter and identified potential cis-regulatory transcription factor (TF) components. In addition, correlation analysis showed that cattle SMAD1 promoter activity of individuals with Hap4 (P < 0.01) was stronger than that of individuals with Hap2. while the transcriptional activity of individuals with Hap3 within SMAD5 gene promoter was significantly (P < 0.01) higher followed by H2. Uniformly, diplotypes H4-H6 of SMAD1 gene and H1-H3 of SMAD5 gene performed significant (P < 0.01) associations with body measurement and improved carcass quality traits. All these results have indicated that polymorphisms in SMAD1 and SMAD5 genes promoter could impact the transcriptional regulation and then affect muscle content in beef cattle. Moreover, both the SMAD1 and SMAD5 genes were expressed ubiquitously in 10 tissues and had higher expression in the longissimus thoracis tissue from 6-month-old and 12-month-old cattle than in cattle of other ages. We can conclude that SMAD1 and SMAD5 genes may play an important role in muscle growth and development, and the variants mapped within SMAD1 and SMAD5 genes can be utilized in molecular marker-assisted selection for cattle carcass quality and body measurement traits in breed improvement programs of Qinchuan cattle.

Identification and Validation of Marketing Weight-Related SNP Markers Using SLAF Sequencing in Male Yangzhou Geese

Growth performance is a complex economic trait for avian production. The swan goose (Anser cygnoides) has never been exploited genetically like chickens or other waterfowl species such as ducks. Traditional phenotypic selection is still the main method for genetic improvement of geese body weight. In this study, specific locus amplified fragment sequencing (SLAF-seq) with bulked segregant analysis (BSA) was conducted for discovering and genotyping single nucleotide polymorphisms (SNPs) associated with marketing weight trait in male geese. A total of 149,045 SNPs were obtained from 427,093 SLAF tags with an average sequencing depth of 44.97-fold and a Q30 value of 93.26%. After SNPs' filtering, a total of 12,917 SNPs were included in the study. The 31 highest significant SNPs-which had different allelic frequencies-were further validated by individual-based AS-PCR genotyping in two populations. The association between 10 novel SNPs and the marketing weight of male geese was confirmed. The 10 significant SNPs were involved in linear regression model analysis, which confirmed single-SNP associations and revealed three types of SNP networks for marketing weight. The 10 significant SNPs were located within or close to 10 novel genes, which were identified. The qPCR analysis showed significant difference between genotypes of each SNP in seven genes. Developed SLAF-seq and identified genes will enrich growth performance studies, promoting molecular breeding applications to boost the marketing weight of Chinese geese.

Introgression contributes to distribution of structural variations in cattle

Structural variations (SVs) are an important source of phenotypic diversity in cattle. Here, 72 whole genome sequences representing taurine and zebu cattle were used to identify SVs. Applying multiple approaches, 16,738 SVs were identified. A comparison against the Database of Genomic Variants archives revealed that 1575 SVs were novel in our data. A novel duplication covering the entire GALNT15 gene, was observed only in N'Dama. A duplication, which was previously reported only in zebu and associated with navel length, was also observed in N'Dama. Investigation of a novel deletion located upstream of CAST13 gene and identified only in Italian cattle and zebu, revealed its introgressed origin in the former. Overall, our data highlights how the SVs distribution in cattle is also shaped by forces such as demographical differences and gene flow. The cattle SVs of this study and its meta-data can be visualized on an interactive genome browser at https://tinyurl.com/svCowArs.

Single-nucleotide polymorphism detecting of some candidate genes related to lipid metabolism in Booroola Merino-Afshari sheep by Bayesian model averaging

This study was designed to identify single-nucleotide polymorphisms (SNPs) of some candidate genes related to lipid metabolism and their association with carcass fat in male crossbred lambs. Hence, 96 of almost 11-month-old Booroola Merino-Afshari crossbred lambs (first-generation backcross) were used by considering their phenotypic carcass traits. Then, DNA was extracted and DNA targets were amplified using designed specific primers by PCR procedure. Identification of potential SNPs was done by a direct sequencing method for LEP, FABP4, DGAT1, GH, and TRIB3 genes using the sequencing-RFLP procedure. Then, the most probable statistical models based on additive and genotypic effects of identified SNPs in each trait were obtained by the Bayesian model averaging (BMA) approach of R software (Ver. 3.3.1) to assess the association of SNPs with traits. Detected SNPs in this study included two SNPs in exon 3 of LEP, one SNP in exon 2 of TRIB3, one SNP in intron 2 of FABP4, one SNP in 5' UTR of DGAT1, and two SNPs in 3' UTR of GH genes. For carcass weight trait, one of the identified SNP genotypes in the LEP (c.587G > A) had a higher probability in the model. Carcass weight of lambs with GA genotype was 2.46 kg heavier than GG genotype. Also, two genes of TRIB3 and GH2 had the highest probability in the models of fat tail and waste weight, respectively. Based on the results, these polymorphisms can be used in the marker-assisted selection of breeding programs and designing DNA chips for genomic selection.

Allelic and genotypic frequencies for loci associated with meat quality in Mexican Braunvieh cattle

The objective was to estimate allelic and genotypic frequencies for loci associated with meat quality in a Mexican population of Braunvieh cattle. Information was obtained from 300 animals genotyped with the Genomic Profile Bovine LD chip of 30K and 50K SNPs. After the final edition, including quality control, the data contained information for 12 loci of the CAPN1, CAPN3, CAPN5, CAPN14, DGAT1, DGAT2, TG, ANK1, and MADH3 genes. Allelic and genotypic frequencies and Hardy-Weinberg equilibrium were estimated with the Cervus 3.0.7 software. The studied population markers were in Hardy-Weinberg equilibrium, except for those associated with CAPN1, DGAT1, and MADH3. Frequencies higher than those reported for other breeds were found for genotypes associated with meat softness, higher marbling score, lower quantity of saturated fatty acids, and lower shear force (CAPN1 and DGAT2). There were similarities with frequencies reported for Bos taurus breeds for the CAPN3 and TG genes. For the DGAT1 and ANK1 genes, the frequencies of the desired genotypes were low. A marker for DGAT1 and another for MADH3 were monomorphic. The results of this study are encouraging in terms of the potential of the Braunvieh population studied for breeding programs aiming to increase meat quality. The breed has strengths that could be used either by crossbreeding to generate heterozygous animals or by selection to increase frequencies of valuable alleles.

Genome-wide association studies detects candidate genes for wool traits by re-sequencing in Chinese fine-wool sheep

BACKGROUND

The quality and yield of wool determine the economic value of the fine-wool sheep. Therefore, discovering markers or genes relevant to wool traits is the cornerstone for the breeding of fine-wool sheep. In this study, we used the Illumina HiSeq X Ten platform to re-sequence 460 sheep belonging to four different fine-wool sheep breeds, namely, Alpine Merino sheep (AMS), Chinese Merino sheep (CMS), Aohan fine-wool sheep (AHS) and Qinghai fine-wool sheep (QHS). Eight wool traits, including fiber diameter (FD), fiber diameter coefficient of variance (FDCV), fiber diameter standard deviation (FDSD), staple length (SL), greasy fleece weight (GFW), clean wool rate (CWR), staple strength (SS) and staple elongation (SE) were examined. A genome-wide association study (GWAS) was performed to detect the candidate genes for the eight wool traits.

RESULTS

A total of 8.222 Tb of raw data was generated, with an average of approximately 8.59X sequencing depth. After quality control, 12,561,225 SNPs were available for analysis. And a total of 57 genome-wide significant SNPs and 30 candidate genes were detected for the desired wool traits. Among them, 7 SNPs and 6 genes are related to wool fineness indicators (FD, FDCV and FDSD), 10 SNPs and 7 genes are related to staple length, 13 SNPs and 7 genes are related to wool production indicators (GFW and CWR), 27 SNPs and 10 genes associated with staple elongation. Among these candidate genes, UBE2E3 and RHPN2 associated with fiber diameter, were found to play an important role in keratinocyte differentiation and cell proliferation. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment results, revealed that multitude significant pathways are related to keratin and cell proliferation and differentiation, such as positive regulation of canonical Wnt signaling pathway (GO:0090263).

CONCLUSION

This is the first GWAS on the wool traits by using re-sequencing data in Chinese fine-wool sheep. The newly detected significant SNPs in this study can be used in genome-selective breeding for the fine-wool sheep. And the new candidate genes would provide a good theoretical basis for the fine-wool sheep breeding.

Selection signatures of Fuzhong Buffalo based on whole-genome sequences

BACKGROUND

Fuzhong buffalo, a native breed of Guangxi Zhuang Autonomous Region, is traditionally used as a draft animal to provide farm power in the rice cultivation. In addition, the Fuzhong buffalo also prepared for the bullfighting festival organized by the locals. The detection of the selective signatures in its genome can help in elucidating the selection mechanisms in its stamina and muscle development of a draft animal.

RESULTS

In this study, we analyzed 27 whole genomes of buffalo (including 15 Fuzhong buffalo genomes and 12 published buffalo genomes from Upper Yangtze region). The ZHp, ZFst, π-Ratio, and XP-EHH statistics were used to identify the candidate signatures of positive selection in Fuzhong buffalo. Our results detected a set of candidate genes involving in the pathways and GO terms associated with the response to exercise (e.g., ALDOA, STAT3, AKT2, EIF4E2, CACNA2D2, TCF4, CDH2), immunity (e.g., PTPN22, NKX2-3, PIK3R1, ITK, TMEM173), nervous system (e.g., PTPN21, ROBO1, HOMER1, MAGI2, SLC1A3, NRG3, SNAP47, CTNNA2, ADGRL3). In addition, we also identified several genes related to production and growth traits (e.g., PHLPP1, PRKN, MACF1, UCN3, RALGAPA1, PHKB, PKD1L). Our results depicted several pathways, GO terms, and candidate genes to be associated with response to exercise, immunity, nervous system, and growth traits.

CONCLUSIONS

The selective sweep analysis of the Fuzhong buffalo demonstrated positive selection pressure on potential target genes involved in behavior, immunity, and growth traits, etc. Our findings provided a valuable resource for future research on buffalo breeding and an insight into the mechanisms of artificial selection.

Genome Regulation and Gene Interaction Networks Inferred From Muscle Transcriptome Underlying Feed Efficiency in Pigs

Improvement of feed efficiency (FE) is key for Sustainability and cost reduction in pig production. Our aim was to characterize the muscle transcriptomic profiles in Danbred Duroc (Duroc; n = 13) and Danbred Landrace (Landrace; n = 28), in relation to FE for identifying potential biomarkers. RNA-seq data on the 41 pigs was analyzed employing differential gene expression methods, gene-gene interaction and network analysis, including pathway and functional analysis. We also compared the results with genome regulation in human exercise data, hypothesizing that increased FE mimics processes triggered in exercised muscle. In the differential expression analysis, 13 genes were differentially expressed, including: MRPS11, MTRF1, TRIM63, MGAT4A, KLH30. Based on a novel gene selection method, the divergent count, we performed pathway enrichment analysis. We found five significantly enriched pathways related to feed conversion ratio (FCR). These pathways were mainly related to mitochondria, and summarized in the mitochondrial translation elongation (MTR) pathway. In the gene interaction analysis, the most interesting genes included the mitochondrial genes: PPIF, MRPL35, NDUFS4 and the fat metabolism and obesity genes: AACS, SMPDL3B, CTNNBL1, NDUFS4, and LIMD2. In the network analysis, we identified two modules significantly correlated with FCR. Pathway enrichment of module genes identified MTR, electron transport chain and DNA repair as enriched pathways. The network analysis revealed the mitochondrial gene group NDUF as key network hub genes, showing their potential as biomarkers. Results show that genes related to human exercise were enriched in identified FCR related genes. We conclude that mitochondrial activity is a key driver for FCR in muscle tissue, and mitochondrial genes could be potential biomarkers for FCR in pigs.

Multi-breed genome-wide association studies across countries for electronically recorded behavior traits in local dual-purpose cows

Basic bovine behavior is a crucial parameter influencing cattle domestication. In addition, behavior has an impact on cattle productivity, welfare and adaptation. The aim of the present study was to infer quantitative genetic and genomic mechanisms contributing to natural dual-purpose cow behavior in grazing systems. In this regard, we genotyped five dual-purpose breeds for a dense SNP marker panel from four different European countries. All cows from the across-country study were equipped with the same electronic recording devices. In this regard, we analyzed 97,049 longitudinal sensor behavior observations from 319 local dual-purpose cows for rumination, feeding, basic activity, high active, not active and ear temperature. According to the specific sensor behaviors and following a welfare protocol, we computed two different welfare indices. For genomic breed characterizations and multi-breed genome-wide association studies, sensor traits and test-day production records were merged with 35,826 SNP markers per cow. For the estimation of variance components, we used the pedigree relationship matrix and a combined similarity matrix that simultaneously included both pedigree and genotypes. Heritabilities for feeding, high active and not active were in a moderate range from 0.16 to 0.20. Estimates were very similar from both relationship matrix-modeling approaches and had quite small standard errors. Heritabilities for the remaining sensor traits (feeding, basic activity, ear temperature) and welfare indices were lower than 0.09. Five significant SNPs on chromosomes 11, 17, 27 and 29 were associated with rumination, and two different SNPs significantly influenced the sensor traits “not active” (chromosome 13) and “feeding” (chromosome 23). Gene annotation analyses inferred 22 potential candidate genes with a false discovery rate lower than 20%, mostly associated with rumination (13 genes) and feeding (8 genes). Mendelian randomization based on genomic variants (i.e., the instrumental variables) was used to infer causal inference between an exposure and an outcome. Significant regression coefficients among behavior traits indicate that all specific behavioral mechanisms contribute to similar physiological processes. The regression coefficients of rumination and feeding on milk yield were 0.10 kg/% and 0.12 kg/%, respectively, indicating their positive influence on dual-purpose cow productivity. Genomically, an improved welfare behavior of grazing cattle, i.e., a higher score for welfare indices, was significantly associated with increased fat and protein percentages.

Single Nucleotide Polymorphisms in SLC5A1, CCNA1, and ABCC1 and the Association with Litter Size in Small-Tail Han Sheep

Simple Summary

Litter size is one of the most important reproductive traits in sheep. Four single nucleotide polymorphisms (SNPs), g.70067210 T > C in SLC5A1, g.25350431 C > T and g.25360220 T > C in CCNA1, and g.14413132 C > T in ABCC1, were identified by mass spectrometry and may be associated with litter size in sheep. Four SNPs were genotyped in Small-Tail Han, Hu, Cele Black, Suffolk, Sunite, Prairie Tibetan, and Tan sheep, and the expression patterns of SLC5A1, CCNA1, and ABCC1 were determined in Small-Tail Han sheep with different fecundities. Furthermore, we also studied the FecB mutation’s association with litter size in Small-Tail Han sheep. The results indicated that all genes included in this study were differentially expressed in the ovary and uterus of polytocous and monotocous Small-Tail Han sheep. Furthermore, association analysis indicated that both g.70067210 T > C in SLC5A1 and the FecB mutation in BMPR-IB were significantly associated with litter size in Small-Tail Han sheep. Linear regression analysis of the association of multiple markers (FecB and g.70067210 T > C in SCL5A1) with litter size indicated that homozygous ewes carrying the BB/TT genotype had a larger litter size than any other genotype.

Abstract

SLC5A1, CCNA1, and ABCC1 have been extensively studied as candidate genes because of their great influence on the reproductive traits of animals. However, little is known about the association between polymorphisms of the SLC5A1, CCNA1, and ABCC1 genes and litter size in Small-Tail Han sheep. In this study, the expression levels of SLC5A1, CCNA1, and ABCC1 in HPG (hypothalamic–pituitary–gonadal) axis tissues of polytocous and monotocous Small-Tail Han sheep were analyzed by qPCR. To better understand the effects of four single nucleotide polymorphisms (SNPs) comprising of g.70067210 T > C in SLC5A1, g.25350431 C > T and g.25360220 T > C in CCNA1, and g.14413132 C > T in ABCC1, a population genetic analysis was conducted using data obtained from genotyping in 728 sheep from seven breeds. The results indicated that all genes included in this study were differentially expressed in the pituitary and uterus of polytocous and monotocous Small-Tail Han sheep (p < 0.05). The associations of these four SNPs and the FecB mutation with litter size in 384 Small-Tail Han sheep were analyzed, therefore, and it was found that both g.70067210T > C and the FecB mutation were significantly associated with litter size (p < 0.05). The linear regression analysis of the association of multiple markers (FecB and g.70067210 T > C in SCL5A1) with litter size indicated that homozygous ewes carrying the BB/TT genotype had larger litter size than any ewes with any other genotype. In conclusion, the SLC5A1 SNPs significantly affect litter size in sheep and are useful as genetic marker for litter size.

Cloning and Expression Levels of TFAM and TFB2M Genes and their Correlation with Meat and Carcass Quality Traits in Jiaxing Black Pig

The coding sequences (CDS) of TFAM and TFB2M genes from Jiaxing Black Pig (JBP) were first obtained by RT-PCR and DNA-seq in the present study. Sequence analyses showed that the TFAM gene contains a 741-bp CDS region encoding 246 amino acids sharing a 100% homology with the sequence on NCBI, while TFB2M gene contains a CDS region of 1176 bp encoding 391 amino acids with two missense mutations. The results of quantitative Real-Time PCR for TFAM and TFB2M revealed that transcripts of the genes were both presented at the highest levels in spleen tissue followed by liver tissue, while the least levels in longissimus dorsi muscle (LDM), and obviously the higher levels in two adipose tissues than those in LDM tissue (P<0.01). Meanwhile, a total of forty-two JBPs were employed in this experiment to investigate the effect of these two genes on the carcass, meat quality traits and flavor substances such as fatty acids, intramuscular fat (IMF) in LDM. As expected, some strong correlations of gene expression abundance of TFAM and TFB2M mRNA in particular tissues such as liver and LDM with carcass and meat quality traits including marbling score, as well as the content of saturated fatty acid (SFA), in JBP were found.

A copy number variant scan in the autochthonous Valdostana Red Pied cattle breed and comparison with specialized dairy populations

Copy number variants (CNVs) are an important source of genomic structural variation, recognized to influence phenotypic variation in many species. Many studies have focused on identifying CNVs within and between human and livestock populations alike, but only few have explored population-genetic properties in cattle based on CNVs derived from a high-density SNP array. We report a high-resolution CNV scan using Illumina’s 777k BovineHD Beadchip for Valdostana Red Pied (VRP), an autochthonous Italian dual-purpose cattle population reared in the Alps that did not undergo strong selection for production traits. After stringent quality control and filtering, CNVs were called across 108 bulls using the PennCNV software. A total of 6,784 CNVs were identified, summarized to 1,723 CNV regions (CNVRs) on 29 autosomes covering a total of ~59 Mb of the UMD3.1 assembly. Among the mapped CNVRs, there were 812 losses, 832 gains and 79 complexes. We subsequently performed a comparison of CNVs detected in the VRP and those available from published studies in the Italian Brown Swiss (IBS) and Mexican Holstein (HOL). A total of 171 CNVRs were common to all three breeds. Between VRP and IBS, 474 regions overlapped, while only 313 overlapped between VRP and HOL, indicating a more similar genetic background among populations with common origins, i.e. the Alps. The principal component, clustering and admixture analyses showed a clear separation of the three breeds into three distinct clusters. In order to describe the distribution of CNVs within and among breeds we used the pair V_ST statistic, considering only the CNVRs shared to more than 5 individuals (within breed). We identified unique and highly differentiated CNVs (n = 33), some of which could be due to specific breed selection and adaptation. Genes and QTL within these regions were characterized.

Genome-wide association study of carcass weight in commercial Hanwoo cattle

Objective

The objective of the present study was to validate genes and genomic regions associated with carcass weight using a low-density single nucleotide polymorphism (SNP) Chip in Hanwoo cattle breed.

Methods

Commercial Hanwoo steers (n = 220) were genotyped with 20K GeneSeek genomic profiler BeadChip. After applying the quality control of criteria of a call rate ≥90% and minor allele frequency (MAF) ≥0.01, a total of 15,235 autosomal SNPs were left for genome-wide association (GWA) analysis. The GWA tests were performed using single-locus mixed linear model. Age at slaughter was fitted as fixed effect and sire included as a covariate. The level of genome-wide significance was set at 3.28×10⁻⁶ (0.05/15,235), corresponding to Bonferroni correction for 15,235 multiple independent tests.

Results

By employing EMMAX approach which is based on a mixed linear model and accounts for population stratification and relatedness, we identified 17 and 16 loci significantly (p<0.001) associated with carcass weight for the additive and dominant models, respectively. The second most significant (p = 0.000049) SNP (ARS-BFGL-NGS-28234) on bovine chromosome 4 (BTA4) at 21 Mb had an allele substitution effect of 43.45 kg. Some of the identified regions on BTA2, 6, 14, 22, and 24 were previously reported to be associated with quantitative trait loci for carcass weight in several beef cattle breeds.

Conclusion

This is the first genome-wide association study using SNP chips on commercial Hanwoo steers, and some of the loci newly identified in this study may help to better DNA markers that determine increased beef production in commercial Hanwoo cattle. Further studies using a larger sample size will allow confirmation of the candidates identified in this study.

Determination of Genetic Structure and Signatures of Selection in Three Strains of Tanzania Shorthorn Zebu, Boran and Friesian Cattle by Genome-Wide SNP Analyses

Background

More than 90 percent of cattle in Tanzania belong to the indigenous Tanzania Short Horn Zebu (TSZ) population which has been classified into 12 strains based on historical evidence, morphological characteristics, and geographic distribution. However, specific genetic information of each TSZ population has been lacking and has caused difficulties in designing programs such as selection, crossbreeding, breed improvement or conservation. This study was designed to evaluate the genetic structure, assess genetic relationships, and to identify signatures of selection among cattle of Tanzania with the main goal of understanding genetic relationship, variation and uniqueness among them.

Methodology/Principal findings

The Illumina Bos indicus SNP 80K BeadChip was used to genotype genome wide SNPs in 168 DNA samples obtained from three strains of TSZ cattle namely Maasai, Tarime and Sukuma as well as two comparative breeds; Boran and Friesian. Population structure and signatures of selection were examined using principal component analysis (PCA), admixture analysis, pairwise distances (F_ST), integrated haplotype score (iHS), identical by state (IBS) and runs of homozygosity (ROH). There was a low level of inbreeding (F~0.01) in the TSZ population compared to the Boran and Friesian breeds. The analyses of F_ST, IBS and admixture identified no considerable differentiation between TSZ trains. Importantly, common ancestry in Boran and TSZ were revealed based on admixture and IBD, implying gene flow between two populations. In addition, Friesian ancestry was found in Boran. A few common significant iHS were detected, which may reflect influence of recent selection in each breed or strain.

Conclusions

Population admixture and selection signatures could be applied to develop conservation plan of TSZ cattle as well as future breeding programs in East African cattle.

Iron Content Affects Lipogenic Gene Expression in the Muscle of Nelore Beef Cattle

Iron (Fe) is an essential mineral for metabolism and plays a central role in a range of biochemical processes. Therefore, this study aimed to identify differentially expressed (DE) genes and metabolic pathways in Longissimus dorsi (LD) muscle from cattle with divergent iron content, as well as to investigate the likely role of these DE genes in biological processes underlying beef quality parameters. Samples for RNA extraction for sequencing and iron, copper, manganese, and zinc determination were collected from LD muscles at slaughter. Eight Nelore steers, with extreme genomic estimated breeding values for iron content (Fe-GEBV), were selected from a reference population of 373 animals. From the 49 annotated DE genes (FDR<0.05) found between the two groups, 18 were up-regulated and 31 down-regulated for the animals in the low Fe-GEBV group. The functional enrichment analyses identified several biological processes, such as lipid transport and metabolism, and cell growth. Lipid metabolism was the main pathway observed in the analysis of metabolic and canonical signaling pathways for the genes identified as DE, including the genes FASN, FABP4, and THRSP, which are functional candidates for beef quality, suggesting reduced lipogenic activities with lower iron content. Our results indicate metabolic pathways that are partially influenced by iron, contributing to a better understanding of its participation in skeletal muscle physiology.

Genome-Wide Association Study of Meat Quality Traits in Nellore Cattle

The objective of this study was to identify genomic regions that are associated with meat quality traits in the Nellore breed. Nellore steers were finished in feedlots and slaughtered at a commercial slaughterhouse. This analysis included 1,822 phenotypic records of tenderness and 1,873 marbling records. After quality control, 1,630 animals genotyped for tenderness, 1,633 animals genotyped for marbling, and 369,722 SNPs remained. The results are reported as the proportion of variance explained by windows of 150 adjacent SNPs. Only windows with largest effects were considered. The genomic regions were located on chromosomes 5, 15, 16 and 25 for marbling and on chromosomes 5, 7, 10, 14 and 21 for tenderness. These windows explained 3,89% and 3,80% of the additive genetic variance for marbling and tenderness, respectively. The genes associated with the traits are related to growth, muscle development and lipid metabolism. The study of these genes in Nellore cattle is the first step in the identification of causal mutations that will contribute to the genetic evaluation of the breed.

Genome Wide Sampling Sequencing for SNP Genotyping: Methods, Challenges and Future Development

Genetic polymorphisms, particularly single nucleotide polymorphisms (SNPs), have been widely used to advance quantitative, functional and evolutionary genomics. Ideally, all genetic variants among individuals should be discovered when next generation sequencing (NGS) technologies and platforms are used for whole genome sequencing or resequencing. In order to improve the cost-effectiveness of the process, however, the research community has mainly focused on developing genome-wide sampling sequencing (GWSS) methods, a collection of reduced genome complexity sequencing, reduced genome representation sequencing and selective genome target sequencing. Here we review the major steps involved in library preparation, the types of adapters used for ligation and the primers designed for amplification of ligated products for sequencing. Unfortunately, currently available GWSS methods have their drawbacks, such as inconsistency in the number of reads per sample library, the number of sites/targets per individual, and the number of reads per site/target, all of which result in missing data. Suggestions are proposed here to improve library construction, genotype calling accuracy, genome-wide marker density and read mapping rate. In brief, optimized GWSS library preparation should generate a unique set of target sites with dense distribution along chromosomes and even coverage per site across all individuals.

Identification of Laying-Related SNP Markers in Geese Using RAD Sequencing

Laying performance is an important economical trait of goose production. As laying performance is of low heritability, it is of significance to develop a marker-assisted selection (MAS) strategy for this trait. Definition of sequence variation related to the target trait is a prerequisite of quantitating MAS, but little is presently known about the goose genome, which greatly hinders the identification of genetic markers for the laying traits of geese. Recently developed restriction site-associated DNA (RAD) sequencing is a possible approach for discerning large-scale single nucleotide polymorphism (SNP) and reducing the complexity of a genome without having reference genomic information available. In the present study, we developed a pooled RAD sequencing strategy for detecting geese laying-related SNP. Two DNA pools were constructed, each consisting of equal amounts of genomic DNA from 10 individuals with either high estimated breeding value (HEBV) or low estimated breeding value (LEBV). A total of 139,013 SNP were obtained from 42,291,356 sequences, of which 18,771,943 were for LEBV and 23,519,413 were for HEBV cohorts. Fifty-five SNP which had different allelic frequencies in the two DNA pools were further validated by individual-based AS-PCR genotyping in the LEBV and HEBV cohorts. Ten out of 55 SNP exhibited distinct allele distributions in these two cohorts. These 10 SNP were further genotyped in a goose population of 492 geese to verify the association with egg numbers. The result showed that 8 of 10 SNP were associated with egg numbers. Additionally, liner regression analysis revealed that SNP Record-111407, 106975 and 112359 were involved in a multiplegene network affecting laying performance. We used IPCR to extend the unknown regions flanking the candidate RAD tags. The obtained sequences were subjected to BLAST to retrieve the orthologous genes in either ducks or chickens. Five novel genes were cloned for geese which harbored the candidate laying-related SNP, including membrane associated guanylate kinase (MAGI-1), KIAA1462, Rho GTPase activating protein 21 (ARHGAP21), acyl-CoA synthetase family member 2 (ACSF2), astrotactin 2 (ASTN2). Collectively, our data suggests that 8 SNP and 5 genes might be promising candidate markers or targets for marker-assisted selection of egg numbers in geese.

Anti-diabetic effect of purple corn extract on C57BL/KsJ db/db mice

BACKGROUND/OBJECTIVES

Recently, anthocyanins have been reported to have various biological activities. Furthermore, anthocyanin-rich purple corn extract (PCE) ameliorated insulin resistance and reduced diabetes-associated mesanginal fibrosis and inflammation, suggesting that it may have benefits for the prevention of diabetes and diabetes complications. In this study, we determined the anthocyanins and non-anthocyanin component of PCE by HPLC-ESI-MS and investigated its anti-diabetic activity and mechanisms using C57BL/KsJ db/db mice.

MATERIALS/METHODS

The db/db mice were divided into four groups: diabetic control group (DC), 10 or 50 mg/kg PCE (PCE 10 or PCE 50), or 10 mg/kg pinitol (pinitol 10) and treated with drugs once per day for 8 weeks. During the experiment, body weight and blood glucose levels were measured every week. At the end of treatment, we measured several diabetic parameters.

RESULTS

Compared to the DC group, Fasting blood glucose levels were 68% lower in PCE 50 group and 51% lower in the pinitol 10 group. Furthermore, the PCE 50 group showed 2- fold increased C-peptide and adiponectin levels and 20% decreased HbA1c levels, than in the DC group. In pancreatic islets morphology, the PCE- or pinitol-treated mice showed significant prevention of pancreatic β-cell damage and higher insulin content. Microarray analyses results indicating that gene and protein expressions associated with glycolysis and fatty acid metabolism in liver and fat tissues. In addition, purple corn extract increased the phosphorylation of AMP-activated protein kinase (AMPK) and decreased phosphoenolpyruvate carboxykinase (PEPCK), glucose 6-phosphatase (G6pase) genes in liver, and also increased glucose transporter 4 (GLUT4) expressions in skeletal muscle.

CONCLUSIONS

Our results suggested that PCE exerted anti-diabetic effects through protection of pancreatic β-cells, increase of insulin secretion and AMPK activation in the liver of C57BL/KsJ db/db mice.

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.

An approach to identify SNPs in the gene encoding acetyl-CoA acetyltransferase-2 (ACAT-2) and their proposed role in metabolic processes in pig

The novel liver protein acetyl-CoA acetyltransferase-2 (ACAT2) is involved in the beta-oxidation and lipid metabolism. Its comprehensive relative expression, in silico non-synonymous single nucleotide polymorphism (nsSNP) analysis, as well as its annotation in terms of metabolic process with another protein from the same family, namely, acetyl-CoA acyltransferase-2 (ACAA2) was performed in Sus scrofa. This investigation was conducted to understand the most important nsSNPs of ACAT2 in terms of their effects on metabolic activities and protein conformation. The two most deleterious mutations at residues 122 (I to V) and 281 (R to H) were found in ACAT2. Validation of expression of genes in the laboratory also supported the idea of differential expression of ACAT2 and ACAA2 conceived through the in silico analysis. Analysis of the relative expression of ACAT2 and ACAA2 in the liver tissue of Jeju native pig showed that the former expressed significantly higher (P<0.05). Overall, the computational prediction supported by wet laboratory analysis suggests that ACAT2 might contribute more to metabolic processes than ACAA2 in swine. Further associations of SNPs in ACAT2 with production traits might guide efforts to improve growth performance in Jeju native pigs.

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.

Detecting novel SNPs and breed-specific haplotypes at calpastatin gene in Iranian fat- and thin-tailed sheep breeds and their effects on protein structure

Calpastatin has been introduced as a potential candidate gene for growth and meat quality traits. In this study, genetic variability was investigated in the exon 6 and its intron boundaries of ovine CAST gene by PCR-SSCP analysis and DNA sequencing. Also a protein sequence and structural analysis were performed to predict the possible impact of amino acid substitutions on physicochemical properties and structure of the CAST protein. A total of 487 animals belonging to four ancient Iranian sheep breeds with different fat metabolisms, Lori-Bakhtiari and Chall (fat-tailed), Zel-Atabay cross-bred (medium fat-tailed) and Zel (thin-tailed), were analyzed. Eight unique SSCP patterns, representing eight different sequences or haplotypes, CAST-1, CAST-2 and CAST-6 to CAST-11, were identified. Haplotypes CAST-1 and CAST-2 were most common with frequency of 0.365 and 0.295. The novel haplotype CAST-8 had considerable frequency in Iranian sheep breeds (0.129). All the consensus sequences showed 98-99%, 94-98%, 92-93% and 82-83% similarity to the published ovine, caprine, bovine and porcine CAST locus sequences, respectively. Sequence analysis revealed four SNPs in intron 5 (C24T, G62A, G65T and T69-) and three SNPs in exon 6 (c.197A>T, c.282G>T and c.296C>G). All three SNPs in exon 6 were missense mutations which would result in p.Gln 66 Leu, p.Glu 94 Asp and p.Pro 99 Arg substitutions, respectively, in CAST protein. All three amino acid substitutions affected the physicochemical properties of ovine CAST protein including hydrophobicity, amphiphilicity and net charge and subsequently might influence its structure and effect on the activity of Ca2+ channels; hence, they might regulate calpain activity and afterwards meat tenderness and growth rate. The Lori-Bakhtiari population showed the highest heterozygosity in the ovine CAST locus (0.802). Frequency difference of haplotypes CAST-10 and CAST-8 between Lori-Bakhtiari (fat-tailed) and Zel (thin-tailed) breeds was highly significant (P<0.001), indicating that these two haplotypes might be breed-specific haplotypes that distinguish between fat-tailed and thin-tailed sheep breeds.

Obesity gene atlas in mammals

Obesity in humans has increased at an alarming rate over the past two decades and has become one of the leading public health problems worldwide. Studies have revealed a large number of genes/markers that are associated with obesity and/or obesity-related phenotypes, indicating an urgent need to develop a central database for helping the community understand the genetic complexity of obesity. In the present study, we collected a total of 1,736 obesity associated loci and created a freely available obesity database, including 1,515 protein-coding genes and 221 microRNAs (miRNAs) collected from four mammalian species: human, cattle, rat, and mouse. These loci were integrated as orthologs on comparative genomic views in human, cattle, and mouse. The database and genomic views are freely available online at: http://www.integratomics-time.com/fat_deposition. Bioinformatics analyses of the collected data revealed some potential novel obesity related molecular markers which represent focal points for testing more targeted hypotheses and designing experiments for further studies. We believe that this centralized database on obesity and adipogenesis will facilitate development of comparative systems biology approaches to address this important health issue in human and their potential applications in animals.

Investigation of four candidate genes (IGF2, JHDM1A, COPB1 and TEF1) for growth rate and backfat thickness traits on SSC2q in Large White pigs

As important quantitative traits, the growth rate and backfat thickness are controlled by multiple genes. The aim of this investigation was to evaluate the effect of the single and multiple SNPs of four candidate genes (IGF2, JHDM1A, COPB1 and TEF-1) on growth rate and backfat thickness. The four candidate genes were mapped on the p arm of SSC 2, and there are several QTLs, such as average daily gain, backfat thickness, an imprinted QTLs affecting muscle mass and fat deposition have been reported in this region. The polymorphisms of these genes were detected using PCR-RFLP methods, mixed procedure was used to analyze the single marker association with the growth and backfat thickness traits, and the gene-gene combination was investigated using multiple-markers analysis. The single marker association analysis indicated that the IGF2 intron-3 g.3072G > A and the substitution g.93G > A of TEF-1 gene were significantly associated with the age at 100 kg (P < 0.05). The JHDM1A 3′UTR g.224C > G, the c.3096C > T polymorphism of COPB1 gene and the substitution g.93G > A of TEF-1 gene were all significantly associated with the backfat at the shoulder (P < 0.05), backfat at the last rib, backfat at the lumbar, and the average backfat thickness, respectively. The multiple-markers analysis indicated that IGF2 and TEF-1 integrated gene networks for the age at 100 kg. Therefore, we can suggest that the polymorphism of IGF2 and TEF-1 gene could be used in marker-assisted selection for the age at 100 kg in Large White pigs.

Bovine dedifferentiated adipose tissue (DFAT) cells: DFAT cell isolation

Dedifferentiated fat cells (DFAT cells) are derived from lipid-containing (mature) adipocytes, which possess the ability to symmetrically or asymmetrically proliferate, replicate, and redifferentiate/transdifferentiate. Robust cell isolation and downstream culture methods are needed to isolate large numbers of DFAT cells from any (one) adipose depot in order to establish population dynamics and regulation of the cells within and across laboratories. In order to establish more consistent/repeatable methodology here we report on two different methods to establish viable DFAT cell cultures: both traditional cell culture flasks and non-traditional (flat) cell culture plates were used for ceiling culture establishment. Adipocytes (maternal cells of the DFAT cells) were easier to remove from flat culture plates than flasks and the flat plates also allowed cloning rings to be utilized for cell/cell population isolation. While additional aspects of usage of flat-bottomed cell culture plates may yet need to be optimized by definition of optimum bio-coating to enhance cell attachment, utilization of flat plate approaches will allow more efficient study of the dedifferentiation process or the DFAT progeny cells. To extend our preliminary observations, dedifferentiation of Wagyu intramuscular fat (IMF)-derived mature adipocytes and redifferentiation ability of DFAT cells utilizing the aforementioned isolation protocols were examined in traditional basal media/differentiation induction media (DMI) containing adipogenic inducement reagents. In the absence of treatment approximately 10% isolated Wagyu IMF-mature adipocytes dedifferentiated spontaneously and 70% DFAT cells displayed protracted adipogenesis 12 d after confluence in vitro. Lipid-free intracellular vesicles in the cytoplasm (vesicles possessing an intact membrane but with no any observable or stainable lipid inside) were observed during redifferentiation. One to 30% DFAT cells redifferentiated into lipid-assimilating adipocytes in the DMI media, with distinct lipid-droplets in the cytoplasm and with no observable lipid-free vesicles inside. Moreover, a high confluence level promoted the redifferentiation efficiency of DFAT cells. Wagyu IMF dedifferentiated DFAT cells exhibited unique adipogenesis modes in vitro, revealing a useful cell model for studying adipogenesis and lipid metabolism.

Use of a bovine genome chip to identify new biological pathways for beef quality in cattle

The accumulation of muscle is largely influenced by the genetic background of cattle. Muscle tissue was collected from the longissimus muscle of Lilu beef cattle at 12, 18, 24 and 30 months old. Using meat quality analysis, we found that the Lilu beef cattle have good production and slaughter performance, the performance meets the criterion of beef cattle. Microarray analysis was able to identify a total of 4,219 genes that are differentially expressed (P ≤ 0.01) between the two groups of cattle (12 vs 18; 18 vs 24; 24 vs 30). Bioinformatics analysis results suggested that most of the differentially expressed genes are involved in the metabolic pathways and neuroactive ligand-receptor interaction pathways. In the future study that aims to look for genes relating to growth and meat quality, we will focus on the genes that have been shown to have a significant variation between groups and are involved in the two pathways.

Quantitative genomics of 30 complex phenotypes in Wagyu x Angus F₁ progeny

In the present study, a total of 91 genes involved in various pathways were investigated for their associations with six carcass traits and twenty-four fatty acid composition phenotypes in a Wagyu×Angus reference population, including 43 Wagyu bulls and their potential 791 F₁ progeny. Of the 182 SNPs evaluated, 102 SNPs that were in Hardy-Weinberg equilibrium with minor allele frequencies (MAF>0.15) were selected for parentage assignment and association studies with these quantitative traits. The parentage assignment revealed that 40 of 43 Wagyu sires produced over 96.71% of the calves in the population. Linkage disequilibrium analysis identified 75 of 102 SNPs derived from 54 genes as tagged SNPs. After Bonferroni correction, single-marker analysis revealed a total of 113 significant associations between 44 genes and 29 phenotypes (adjusted P<0.05). Multiple-marker analysis confirmed single-gene associations for 10 traits, but revealed two-gene networks for 9 traits and three-gene networks for 8 traits. Particularly, we observed that TNF (tumor necrosis factor) gene is significantly associated with both beef marbling score (P=0.0016) and palmitic acid (C16:0) (P=0.0043), RCAN1 (regulator of calcineurin 1) with rib-eye area (P=0.0103), ASB3 (ankyrin repeat and SOCS box-containing 3) with backfat (P=0.0392), ABCA1 (ATP-binding cassette A1) with both palmitic acid (C16:0) (P=0.0025) and oleic acid (C18:1n9) (P=0.0114), SLC27A1(solute carrier family 27 A1) with oleic acid (C18:1n9) (P=0.0155), CRH (corticotropin releasing hormone) with both linolenic acid (OMEGA-3) (P=0.0200) and OMEGA 6:3 RATIO (P=0.0054), SLC27A2 (solute carrier family 27 A2) with both linoleic acid (OMEGA-6) (P=0.0121) and FAT (P=0.0333), GNG3 (guanine nucleotide binding protein gamma 3 with desaturase 9 (P=0.0115), and EFEMP1 (EGF containing fibulin-like extracellular matrix protein 1), PLTP (phospholipid transfer protein) and DSEL (dermatan sulfate epimerase-like) with conjugated linoleic acid (P=0.0042-0.0044), respectively, in the Wagyu x Angus F₁ population. In addition, we observed an interesting phenomenon that crossbreeding of different breeds might change gene actions to dominant and overdominant modes, thus explaining the origin of heterosis. The present study confirmed that these important families or pathway-based genes are useful targets for improving meat quality traits and healthful beef products in cattle.

Investigation of four porcine candidate genes (H-FABP, MYOD1, UCP3 and MASTR) for meat quality traits in Large White pigs

Meat quality traits are economically important traits of swine, and are controlled by multiple genes as complex quantitative traits. In the present study four genes, H-FABP (heart fatty acid-binding protein), MASTR (MEF2 activating motif and SAP domain containing transcriptional regulator), UCP3 (uncoupling protein 3) and MYOD1 (myogenic differentiation 1) were researched in Large White pigs. The polymorphisms H-FABP T/C of 5'UTR, MYOD1 g.257 A>C, UCP3 g.1406 G>A in exon 3 and MASTR c.187 C>T have been reported to be associated with meat quality traits in pigs. The aim of this study was to analyze the effect of single and multiple markers for single traits in Large White pigs. The single marker association analysis showed that the H-FABP and MASTR genes were associated with IMF (intramuscular fat content) (P < 0.05), and that the g.257 A>C of MYOD1 gene was most significantly related to muscle pH value (P < 0.01). The multiple markers for IMF were analyzed by combining the markers and quantitative trait modes into the linear regression. The results revealed that H-FABP and MASTR integrate gene networks for IMF. Thus, our study results suggested that H-FABP and MASTR polymorphisms could be used as genetic markers in the marker-assisted selection towards the improvement of IMF in Large White pigs.

Identification of candidate genes related to bovine marbling using protein-protein interaction networks

Complex traits are determined by the combined effects of many loci and are affected by gene networks or biological pathways. Systems biology approaches have an important role in the identification of candidate genes related to complex diseases or traits at the system level. The present study systemically analyzed genes associated with bovine marbling score and identified their relationships. The candidate nodes were obtained using MedScan text-mining tools and linked by protein-protein interaction (PPI) from the Human Protein Reference Database (HPRD). To determine key node of marbling, the degree and betweenness centrality (BC) were used. The hub nodes and biological pathways of our network are consistent with the previous reports about marbling traits, and also suggest unknown candidate genes associated with intramuscular fat. Five nodes were identified as hub genes, which was consistent with the network analysis using quantitative reverse-transcription PCR (qRT-PCR). Key nodes of the PPI network have positive roles (PPARγ, C/EBPα, and RUNX1T1) and negative roles (RXRA, CAMK2A) in the development of intramuscular fat by several adipogenesis-related pathways. This study provides genetic information for identifying candidate genes for the marbling trait in bovine.

The heparan and heparin metabolism pathway is involved in regulation of fatty acid composition

Six genes involved in the heparan sulfate and heparin metabolism pathway, DSEL (dermatan sulfate epimerase-like), EXTL1 (exostoses (multiple)-like 1), HS6ST1 (heparan sulfate 6-O-sulfotransferase 1), HS6ST3 (heparan sulfate 6-O-sulfotransferase 3), NDST3 (N-deacetylase/N-sulfotransferase (heparan glucosaminyl) 3), and SULT1A1 (sulfotransferase family, cytosolic, 1A, phenol-preferring, member 1), were investigated for their associations with muscle lipid composition using cattle as a model organism. Nineteen single nucleotide polymorphisms (SNPs)/multiple nucleotide length polymorphisms (MNLPs) were identified in five of these six genes. Six of these mutations were then genotyped on 246 Wagyu x Limousin F₂ animals, which were measured for 5 carcass, 6 eating quality and 8 fatty acid composition traits. Association analysis revealed that DSEL, EXTL1 and HS6ST1 significantly affected two stearoyl-CoA desaturase activity indices, the amount of conjugated linoleic acid (CLA), and the relative amount of saturated fatty acids (SFA) and monounsaturated fatty acids (MUFA) in skeletal muscle (P<0.05). In particular, HS6ST1 joined our previously reported SCD1 and UQCRC1 genes to form a three gene network for one of the stearoyl-CoA desaturase activity indices. These results provide evidence that genes involved in heparan sulfate and heparin metabolism are also involved in regulation of lipid metabolism in bovine muscle. Whether the SNPs affected heparan sulfate proteoglycan structure is unknown and warrants further investigation.

The reverse cholesterol transport pathway improves understanding of genetic networks for fat deposition and muscle growth in beef cattle

In the present study, thirteen genes involved in the reverse cholesterol transport (RCT) pathway were investigated for their associations with three fat depositions, eight fatty acid compositions and two growth-related phenotypes in a Wagyu x Limousin reference population, including 6 F₁ bulls, 113 F₁ dams, and 246 F₂ progeny. A total of 37 amplicons were used to screen single nucleotide polymorphisms (SNPs) on 6 F₁ bulls. Among 36 SNPs detected in 11 of these 13 genes, 19 were selected for genotyping by the Sequenom assay design on all F₂ progeny. Single-marker analysis revealed seven SNPs in ATP binding cassette A1, apolipoproteins A1, B and E, phospholipid transfer protein and paraoxinase 1 genes significantly associated with nine phenotypes (P<0.05). Previously, we reported genetic networks associated with 19 complex phenotypes based on a total of 138 genetic polymorphisms derived from 71 known functional genes. Therefore, after Bonferroni correction, these significant (adjusted P<0.05) and suggestive (adjusted P<0.10) associations were then used to identify genetic networks related to the RCT pathway. Multiple-marker analysis suggested possible genetic networks involving the RCT pathway for kidney-pelvic-heart fat percentage, rib-eye area, and subcutaneous fat depth phenotypes with markers derived from paraoxinase 1, apolipoproteins A1 and E, respectively. The present study confirmed that genes involved in cholesterol homeostasis are useful targets for investigating obesity in humans as well as for improving meat quality phenotypes in a livestock production.

Lipid metabolism, adipocyte depot physiology and utilization of meat animals as experimental models for metabolic research

Meat animals are unique as experimental models for both lipid metabolism and adipocyte studies because of their direct economic value for animal production. This paper discusses the principles that regulate adipogenesis in major meat animals (beef cattle, dairy cattle, and pigs), the definition of adipose depot-specific regulation of lipid metabolism or adipogenesis, and introduces the potential value of these animals as models for metabolic research including mammary biology and the ontogeny of fatty livers.

Allied industry approaches to alter intramuscular fat content and composition in beef animals

Biochemical and biophysical research tools are used to define the developmental dynamics of numerous cell lineages from a variety of tissues relevant to meat quality. With respect to the adipose cell lineage, much of our present understanding of adipogenesis and lipid metabolism was initially determined through the use of these methods, even though the in vitro or molecular environments are far removed from the tissues of meat animals. This concise review focuses on recent cellular and molecular biology-related research with adipocytes, and how the research might be extended to the endpoint of altering red meat quality. Moreover, economic and policy impacts of such in animal production regimens is discussed. These issues are important, not only with respect to palatability, but also to offer enhanced health benefits to the consumer by altering content of bioactive components in adipocytes.