Genome-wide analysis reveals the effects of artificial selection on production and meat quality traits in Qinchuan cattle

Genetic diversity of United States Rambouillet, Katahdin and Dorper sheep

BACKGROUND

Managing genetic diversity is critically important for maintaining species fitness. Excessive homozygosity caused by the loss of genetic diversity can have detrimental effects on the reproduction and production performance of a breed. Analysis of genetic diversity can facilitate the identification of signatures of selection which may contribute to the specific characteristics regarding the health, production and physical appearance of a breed or population. In this study, breeds with well-characterized traits such as fine wool production (Rambouillet, N = 745), parasite resistance (Katahdin, N = 581) and environmental hardiness (Dorper, N = 265) were evaluated for inbreeding, effective population size (Ne), runs of homozygosity (ROH) and Wright's fixation index (FST) outlier approach to identify differential signatures of selection at 36,113 autosomal single nucleotide polymorphisms (SNPs).

RESULTS

Katahdin sheep had the largest current Ne at the most recent generation estimated with both the GONe and NeEstimator software. The most highly conserved ROH Island was identified in Rambouillet with a signature of selection on chromosome 6 containing 202 SNPs called in an ROH in 50 to 94% of the individuals. This region contained the DCAF16, LCORL and NCAPG genes that have been previously reported to be under selection and have biological roles related to milk production and growth traits. The outlier regions identified through the FST comparisons of Katahdin with Rambouillet and Dorper contained genes with known roles in milk production and mastitis resistance or susceptibility, and the FST comparisons of Rambouillet with Katahdin and Dorper identified genes related to wool growth, suggesting these traits have been under natural or artificial selection pressure in these populations. Genes involved in the cytokine-cytokine receptor interaction pathways were identified in all FST breed comparisons, which indicates the presence of allelic diversity between these breeds in genomic regions controlling cytokine signaling mechanisms.

CONCLUSIONS

In this paper, we describe signatures of selection within diverse and economically important U.S. sheep breeds. The genes contained within these signatures are proposed for further study to understand their relevance to biological traits and improve understanding of breed diversity.

Genome-wide analysis reveals genomic diversity and signatures of selection in Qinchuan beef cattle

BACKGROUND

Indigenous Chinese cattle have abundant genetic diversity and a long history of artificial selection, giving local breeds advantages in adaptability, forage tolerance and resistance. The detection of selective sweeps and comparative genome analysis of selected breeds and ancestral populations provide a basis for understanding differences among breeds and for the identification and utilization of candidate genes. We investigated genetic diversity, population structure, and signatures of selection using genome-wide sequencing data for a new breed of Qinchuan cattle (QNC, n = 21), ancestral Qinchuan cattle (QCC, n = 20), and Zaosheng cattle (ZSC, n = 19).

RESULTS

A population structure analysis showed that the ancestry components of QNC and ZSC were similar. In addition, the QNC and ZSC groups showed higher proportions of European taurine ancestry than that of QCC, and this may explain the larger body size of QNC, approaching that of European cattle under long-term domestication and selection. A neighbor-joining tree revealed that QCC individuals were closely related, whereas QNC formed a distinct group. To search for signatures of selection in the QNC genome, we evaluated nucleotide diversity (θπ), the fixation index (FST) and Tajima's D. Overlapping selective sweeps were enriched for one KEGG pathway, the apelin signaling pathway, and included five candidate genes (MEF2A, SMAD2, CAMK4, RPS6, and PIK3CG). We performed a comprehensive review of genomic variants in QNC, QCC, and ZSC using whole-genome sequencing data. QCC was rich in novel genetic diversity, while diversity in QNC and ZSC cattle was reduced due to strong artificial selection, with divergence from the original cattle.

CONCLUSIONS

We identified candidate genes associated with production traits. These results support the success of selective breeding and can guide further breeding and resource conservation of Qinchuan cattle.

RNA-Seq explores the functional role of the fibroblast growth factor 10 gene in bovine adipocytes differentiation

OBJECTIVE

The present study was executed to explore the molecular mechanism of fibroblast growth factor 10 (FGF10) gene in bovine adipogenesis.

METHODS

The bovine FGF10 gene was overexpressed through Ad-FGF10 or inhibited through siFGF10 and their negative control (NC) in bovine adipocytes, and the multiplicity of infection, transfection efficiency, interference efficiency were evaluated through quantitative real-time polymerase chain reaction, western blotting and fluorescence microscopy. The lipid droplets, triglycerides (TG) content and the expression levels of adipogenic marker genes were measured during preadipocytes differentiation. The differentially expressed genes were explored through deep RNA sequencing.

RESULTS

The highest mRNA level was found in omasum, subcutaneous fat, and intramuscular fat. Moreover, the highest mRNA level was found in adipocytes at day 4 of differentiation. The results of red-oil o staining showed that overexpression (Ad-FGF10) of the FGF10 gene significantly (p<0.05) reduced the lipid droplets and TG content, and their downregulation (siFGF10) increased the measurement of lipid droplets and TG in differentiated bovine adipocytes. Furthermore, the overexpression of the FGF10 gene down regulated the mRNA levels of adipogenic marker genes such as CCAAT enhancer binding protein alpha (C/EBPα), fatty acid binding protein (FABP4), peroxisome proliferator-activated receptor-γ (PPARγ), lipoprotein lipase (LPL), and Fas cell surface death receptor (FAS), similarly, down-regulation of the FGF10 gene enriched the mRNA levels of C/EBPα, PPARγ, FABP4, and LPL genes (p<0.01). Additionally, the protein levels of PPARγ and FABP4 were reduced (p<0.05) in adipocytes infected with Ad-FGF10 gene and enriched in adipocytes transfected with siFGF10. Moreover, a total of 1,774 differentially expressed genes (DEGs) including 157 up regulated and 1,617 down regulated genes were explored in adipocytes infected with Ad-FGF10 or Ad-NC through deep RNA-sequencing. The top Kyoto encyclopedia of genes and genomes pathways regulated through DEGs were the PPAR signaling pathway, cell cycle, base excision repair, DNA replication, apoptosis, and regulation of lipolysis in adipocytes.

CONCLUSION

Therefore, we can conclude that the FGF10 gene is a negative regulator of bovine adipogenesis and could be used as a candidate gene in marker-assisted selection.

Novel polymorphism at KLF15 gene and its association with growth traits in Hu sheep

Growth traits are important economic characteristics of livestock and poultry. In the present study, the expression features of KLF15 and the relationship between KLF15 gene polymorphisms and growth traits in Hu sheep were investigated by using real-time quantitative PCR technology (qPCR), Sanger sequencing, and Kaspar genotyping technology. The qPCR results showed that the KLF15 gene is expressed widely in the tested tissues of Hu sheep, and the expression level of the KLF15 gene in the heart and the muscle was significantly higher than in other tissues (p < 0.05). Missense mutation c.62565119 A > G was found in KLF15, and an association analysis showed that it was correlated with the growth traits (body weight, body height, and body length) of Hu sheep (p < 0.05). The body weight, body height, and body length of the sheep carrying the AA genotype were remarkably higher than those of the GG and AG genotypes (p < 0.05). These results showed that novel polymorphisms at the KLF15 gene can be used as a genetic marker of growth traits of Hu sheep.

Roles of MEF2A and HOXA5 in the transcriptional regulation of the bovine FoxO1 gene

The Forkhead box factor 1 (FoxO1) gene plays a vital role in the growth and development of skeletal muscle. In the present study, expression analysis of the bovine FoxO1 gene exhibited the highest expression in longissimus dorsi muscle followed by its expression in adipose tissue. Moreover, high mRNA expression of FoxO1 gene was found in differentiated bovine myoblasts and adipocytes at day 6 of induced differentiation (p < 0.05). The regulatory pattern of the bovine FoxO1 gene was investigated through screening and dual-luciferase activity of the 1.7 kb 5'UTR (untranslated region) within pGL3-basic vector and a core promoter region was explored at (-285/-27) upstream of the transcription start site. The transcription factors (TFs) MEF2A and HOXA5 within the core promoter region (-285/-27) were found as the regulatory cis-acting element. The siRNA interference of the TFs, chromatin immunoprecipitation (ChIP) assay, and site-directed mutation validated that MEF2A and HOXA5 binding occurs in the region -285/-27 bp and performs an essential role in the transcriptional regulation of bovine FoxO1 gene. These findings explored the regulatory network mechanism of the FoxO1 gene in skeletal muscle development and adipogenesis for the bovine breed improvement program.

Genetic variants of TORC1 gene promoter and their association with carcass quality and body measurement traits in Qinchuan beef cattle

In the present study, sequencing of TORC1 prompter region explored three SNPs at loci g.80G>T, g.93A>T, and g.1253G>A. The SNP1 produced GG, GT and TT, SNP2 AA, AT and TT, and SNP3 produced GG, GA and AA genotypes. Allelic and genotypic frequencies analysis exhibited that SNP1 is within Hardy-Weinberg equilibrium (HWE). All three SNPs were found highly polymorphic as PIC value (0.25 < PIC < 0.50). At loci g.80G>T the cattle with genotype GG showed significantly (P <0.01) larger body length (BL), Wither height (WH), Hip height (HH), Rump length (RL), Hip width (HW), Chest depth (CD), and Chest circumference (CC). The genotype AA at g.93A>T showed significantly (P< 0.01 and 0.05) Larger body length (BL), Wither height (WH), Hip height, Rump length (RL), Hip width (HW), Chest depth (CD), and Chest circumference (CC). Interestingly, the carcass quality parameters such as Ultrasound loin area (ULA) and Intramuscular fat percentage (IF%) was highest in genotype GG at loci g.1253G>A. These findings conclude that genotype GG at loci g.80 G>T and AA at loci g.93A>T could be used as genetic markers for body measurement and genotype GG at loci g.1253G>A for carcass quality traits of TORC1 gene in Qinchuan beef cattle.

The body weight heritability and the effect of non-genetic factors on the body measurement traits in Qinchuan cattle

This study aimed to understand the influence of the non-genetic factors that include breeding year, season, and sex of growth and development traits of Qinchuan cattle and to estimate the heritability of body weight at different growth stages. The Qinchuan cattle measurement records were by the Experiment farm of the National Beef Cattle Improvement Center (Yangling, China) from 2000 to 2017. SPSS and R software were used to analyze the influence of non-genetic factors on body size traits that include body weight (BW), withers height (WH), hip height (HH), body length (BL), chest circumference (CC), abdominal girth (AG), and calf girth (CG), at birth, 6, 12, 18, and 24 months of age. Meanwhile, the single-trait animal model of DMU software was used to estimate the variance component and the heritability of body weight. The results of GLM analysis showed as follows: sex, birth year, and birth season had effects on the body size traits of Qinchuan cattle at different growth stages. Respectively, the heritability of body weight at birth, 6, 12, 18, and 24 months of age were 0.43, 0.32, 0.37, 0.32, and 0.38.

The growth curve determination and economic trait correlation for Qinchuan bull population

Carcass weight, as a measure of meat yield, and body measurements are directly correlated traits in livestock. However, longitudinally collected phenotype records of local breeds are not comprehensive. The research was performed on Qinchuan bull population to understand their growth and development, and data from Qinchuan bull that was weighed and measured at birth, 6, 12, 18, and 24 months of age was analyzed. Furthermore, Logistic, Brody, Gompertz, and Bertallanffy were used to fit the growth curves for weight and body size traits. The results showed that the four curve models have good fitting degrees for the weight and body size (R2 > 0.99), and the Bertallanffy model exhibited a good fit to the measured data of body weight, and the model estimated the inflection point of body weight as (5.43 months of age, 122.01 kg). Particularly, the limited mature body weight can reach 557.8 kg by the Brody model. Body weight was significantly positively correlated with body height, hip height, body length, chest circumference, abdominal girth, and calf girth (p < 0.0001), and the correlation between body weight and body length was the highest (r = 0.975). The regression equation predicting body weight was Y = -275.691 + 3.28 X3 + 1.311 X4 - 0.397 X5.

Selection signatures of Qinchuan cattle based on whole-genome sequences

Qinchuan cattle has gradually improved in body shape and growth rate in the long-term breeding process from the draft cattle to beef cattle. As the head of the five local yellow cattle in China, the Qinchuan cattle has been designated as a specialized beef cattle breed. We investigated the selection signatures using whole genome sequencing data in Qinchuan cattle. Based on Fst, we detected hundreds of candidate genes under selection across Qinchuan, Red Angus, and Japanese Black cattle. Through protein-protein interaction analysis and functional annotation of candidate genes, the results revealed that KMT2E, LTBP1 and NIPBL were related to brain size, body characteristics, and limb development, respectively, suggesting that these potential genes may affect the growth and development traits in Qinchuan cattle. ARIH2, DACT1 and DNM2, et al. are related to meat quality. Meanwhile, TBXA2R can be used as a gene associated with reproductive function, and USH2A affect coat color. This provided a glimpse into the formation of breeds and molecular genetic breeding. Our findings will promote genome-assisted breeding to improve animal production and health.

Bioinformatics and genetic variants analysis of FGF10 gene promoter with their association at carcass quality and body measurement traits in Qinchuan beef cattle

The fibroblast growth factor 10 (FGF10) gene regulates adipogenesis and myogensis. In this study, sequencing of FGF10 prompter region identified three SNPs at loci g.78G > A, g.116C > T and g.201A > T. Each SNP yields three genotypes as GG, GA and AA at loci g.78G > A, CC, CT and TT at loci g.116C > T and AA, AT and TT at loci g.201A > T. Allelic and genotypic frequencies of all three SNPs deviated from the Hardy-Weinberg equilibrium (HWE) (P < 0.05) and were found highly polymorphic as PIC (0.25 < PIC < 0.50). Moreover, we found highest LD (D'/γ2) between SNP2 and SNP3 (0.989/0.909), followed by SNP1 and SNP3 (0.944/0.796). Moreover, three variants of FGF10 gene promoter exhibited significant (P < 0.05) association with body measurement and carcass quality traits in Qinchuan beef cattle. At loci g.78G > A, the genotype GG showed significantly (P < 0.01) larger body length (BL), rump length (RL), chest depth (CD), chest circumference (CC) and ultrasound loin area (ULA). The genotype TC at loci g.116C > T showed significantly (P < 0.01 and 0.05) larger body measurement and intramuscular fat, and ultrasound loin area (ULA). In addition to that, at loci g.201A > T, genotype TT showed significantly (P < 0.01 and P < 0.05) larger body length (BL), rump length (RL), hip width (HW), chest circumference (CC) and ultrasound loin area (ULA). Additionally, screening of promoter sequence of FGF10 gene explored loss of four TFs binding sites (KLF3, ZNF37α, GLIS2 and BCL11A) at g.116C > T because of SNP2. However, a single TF binding site was lost at g.202A > T due to SNP3. Interestingly, none of TF binding site was lost at g.78G > A in SNP1; however, one new TF binding site was gained at this location due to SNP1. These findings conclude that genotype GG, TC and TT could be used as genetic markers of FGF10 gene for body measurement and carcass quality traits in Qinchuan beef cattle.

Genome-wide association study for primal cut lean traits in Canadian beef cattle

This study identified genomic variants and underlying candidate genes related to the whole carcass and individual primal cut lean content in Canadian commercial crossbred beef cattle. Genotyping information of 1035 crossbred beef cattle were available alongside estimated and actual carcass lean meat yield and individual primal cut lean content in all carcasses. Significant fixed effects and covariates were identified and included in the animal model. Genome-wide association analysis were implemented using the weighted single-step genomic best linear unbiased prediction (WssGBLUP). A number of candidate genes identified linked to lean tissue production were unrelated to estimated lean meat yield and were specific to the actual lean traits. Among these, 41 genes were common for actual lean traits, on specific regions of BTA4, BTA13 and BTA25 indicating potential involvement in lean mass synthesis. Therefore, the results suggested the inclusion of primal cut lean traits as a selection objective in breeding programs with consideration of further functional studies of the identified genes could help in optimizing lean yield for maximal carcass value.

Identification of Missense Variants Affecting Carcass Traits for Hanwoo Precision Breeding

This study aimed to identify causal variants associated with important carcass traits such as weight and meat quality in Hanwoo cattle. We analyzed missense mutations extracted from imputed sequence data (ARS-UCD1.2) and performed an exon-specific association test on the carcass traits of 16,970 commercial Hanwoo. We found 33, 2, 1, and 3 significant SNPs associated with carcass weight (CW), backfat thickness (BFT), eye muscle area (EMA), and marbling score (MS), respectively. In CW and EMA, the most significant missense SNP was identified at 19,524,263 on BTA14 and involved the PRKDC. A missense SNP in the ZFAND2B, located at 107,160,304 on BTA2 was identified as being involved in BFT. For MS, missense SNP in the ACVR2B gene, located at 11,849,704 in BTA22 was identified as the most significant marker. The contribution of the most significant missense SNPs to genetic variance was confirmed to be 8.47%, 2.08%, 1.73%, and 1.19% in CW, BFT, EMA, and MS, respectively. We generated favorable and unfavorable haplotype combinations based on the significant SNPs for CW. Significant differences in GEBV (Genomic Estimated Breeding Values) were observed between groups with each favorable and unfavorable haplotype combination. In particular, the missense SNPs in PRKDC, MRPL9, and ANKFN1 appear to significantly affect the protein's function and structure, making them strong candidates as causal mutations. These missense SNPs have the potential to serve as valuable markers for improving carcass traits in Hanwoo commercial farms.

SLC16A7 Promotes Triglyceride Deposition by De Novo Lipogenesis in Chicken Muscle Tissue

Triglyceride (TG) content in chicken muscle tissue signifies intramuscular fat (IMF) content, which is important for improving meat quality. However, the genetic basis of TG deposition in chicken is still unclear. Using 520 chickens from an artificially selected line with significantly increased IMF content and a control line, a genome-wide association study (GWAS) with TG content reports a region of 802 Kb located in chromosome 1. The XP-EHH and gene expression analysis together reveal that the solute carrier family 16 member A7 (SLC16A7) gene is the key candidate gene associated with TG content in chicken muscle tissue. Furthermore, the weighted gene co-expression network analysis (WGCNA) confirmed the regulatory effects of SLC16A7 on promoting TG deposition by de novo lipogenesis (DNL). Functional verification of SLC16A7 in vitro also supports this view, and reveals that this effect mainly occurs in myocytes. Our data highlight a potential IMF deposition pathway by DNL, induced by SLC16A7 in chicken myocytes. These findings will improve the understanding of IMF regulation in chicken and guide the formulation of breeding strategies for high-quality chicken.

Weighted Single-Step Genome-Wide Association Study Uncovers Known and Novel Candidate Genomic Regions for Milk Production Traits and Somatic Cell Score in Valle del Belice Dairy Sheep

Simple Summary

Milk production is the most economically crucial dairy sheep trait and constitutes the major genetic enhancement purpose via selective breeding. Also, mastitis is one of the most frequently encountered diseases, having a significant impact on animal welfare, milk yield, and quality. The aim of this study was to identify genomic region(s) associated with the milk production traits and somatic cell score (SCS) in Valle del Belice sheep using single-step genome-wide association (ssGWA) and genotyping data from medium density SNP panels. We identified several genomic regions (OAR1, OAR2, OAR3, OAR4, OAR6, OAR9, and OAR25) and candidate genes implicated in milk production traits and SCS. Our findings offer new insights into the genetic basis of milk production traits and SCS in dairy sheep.

Abstract

The objective of this study was to uncover genomic regions explaining a substantial proportion of the genetic variance in milk production traits and somatic cell score in a Valle del Belice dairy sheep. Weighted single-step genome-wide association studies (WssGWAS) were conducted for milk yield (MY), fat yield (FY), fat percentage (FAT%), protein yield (PY), protein percentage (PROT%), and somatic cell score (SCS). In addition, our aim was also to identify candidate genes within genomic regions that explained the highest proportions of genetic variance. Overall, the full pedigree consists of 5534 animals, of which 1813 ewes had milk data (15,008 records), and 481 ewes were genotyped with a 50 K single nucleotide polymorphism (SNP) array. The effects of markers and the genomic estimated breeding values (GEBV) of the animals were obtained by five iterations of WssGBLUP. We considered the top 10 genomic regions in terms of their explained genomic variants as candidate window regions for each trait. The results showed that top ranked genomic windows (1 Mb windows) explained 3.49, 4.04, 5.37, 4.09, 3.80, and 5.24% of the genetic variances for MY, FY, FAT%, PY, PROT%, and total SCS, respectively. Among the candidate genes found, some known associations were confirmed, while several novel candidate genes were also revealed, including PPARGC1A, LYPLA1, LEP, and MYH9 for MY; CACNA1C, PTPN1, ROBO2, CHRM3, and ERCC6 for FY and FAT%; PCSK5 and ANGPT1 for PY and PROT%; and IL26, IFNG, PEX26, NEGR1, LAP3, and MED28 for SCS. These findings increase our understanding of the genetic architecture of six examined traits and provide guidance for subsequent genetic improvement through genome selection.

Comparative Enhancer Map of Cattle Muscle Genome Annotated by ATAC-Seq

Annotating regulatory elements could benefit the interpretation of the molecular mechanism of genome-wide association study (GWAS) hits. In this work, we performed transposase-accessible chromatin with sequencing (ATAC-seq) to annotate the cattle muscle genome's functional elements. A total of 10,023 and 11,360 peaks were revealed in muscle genomes of adult and embryo cattle, respectively. The two peak sets produced 8,850 differentially accessible regions (DARs), including 2,515 promoters and 4,319 putative enhancers. These functional elements were associated with the cell cycle, muscle development, and lipid metabolism. A total of 15 putative enhancers were selected for a dual-luciferase reporter assay, and 12 of them showed enhancer activity in cattle myoblasts. Interestingly, the GeneHancer database has annotated the interactions of eight active enhancers with gene promoters, such as embryo-specific peak1053 (log2FC = 1.81, embryo/adult, E/A) with ligand-dependent nuclear receptor corepressor-like protein (LCORL) and embryo-specific peak4218 (log2FC = 1.81) with FERM domain-containing 8 (FRMD8). A total of 295 GWAS loci from the animal QTL database were mapped to 183 putative enhancers, including rs109554838 (associated with cattle body weight and average daily gain) to peak1053 and rs110294629 (associated with beef shear force and tenderness score) to peak4218. Notably, peak4218 has been found to be involved in mouse embryo development. Deleting peak4218 clearly reduced luciferase activity (P = 3.30E-04). Our comparative enhancer map is expected to benefit the area of beef cattle breeding.

Integration of selection signatures and multi-trait GWAS reveals polygenic genetic architecture of carcass traits in beef cattle

Carcass merits are widely considered as economically important traits affecting beef production in the beef cattle industry. However, the genetic basis of carcass traits remains to be well understood. Here, we applied multiple methods, including the Composite of Likelihood Ratio (CLR) and Genome-wide Association Study (GWAS), to explore the selection signatures and candidate variants affecting carcass traits. We identified 11,600 selected regions overlapping with 2214 candidate genes, and most of those were enriched in binding and gene regulation. Notably, we identified 66 and 110 potential variants significantly associated with carcass traits using single-trait and multi-traits analyses, respectively. By integrating selection signatures with single and multi-traits associations, we identified 12 and 27 putative genes, respectively. Several highly conserved missense variants were identified in OR5M13D, NCAPG, and TEX2. Our study supported polygenic genetic architecture of carcass traits and provided novel insights into the genetic basis of complex traits in beef cattle.

Genomics of Adaptations in Ungulates

Ungulates are a group of hoofed animals that have long interacted with humans as essential sources of food, labor, clothing, and transportation. These consist of domesticated, feral, and wild species raised in a wide range of habitats and biomes. Given the diverse and extreme environments inhabited by ungulates, unique adaptive traits are fundamental for fitness. The documentation of genes that underlie their genomic signatures of selection is crucial in this regard. The increasing availability of advanced sequencing technologies has seen the rapid growth of ungulate genomic resources, which offers an exceptional opportunity to understand their adaptive evolution. Here, we summarize the current knowledge on evolutionary genetic signatures underlying the adaptations of ungulates to different habitats.

Opportunities and Challenges for Improving the Productivity of Swamp Buffaloes in Southeastern Asia

The swamp buffalo is a domesticated animal commonly found in Southeast Asia. It is a highly valued agricultural animal for smallholders, but the production of this species has unfortunately declined in recent decades due to rising farm mechanization. While swamp buffalo still plays a role in farmland cultivation, this species' purposes has shifted from draft power to meat, milk, and hide production. The current status of swamp buffaloes in Southeast Asia is still understudied compared to its counterparts such as the riverine buffaloes and cattle. This review discusses the background of swamp buffalo, with an emphasis on recent work on this species in Southeast Asia, and associated genetics and genomics work such as cytogenetic studies, phylogeny, domestication and migration, genetic sequences and resources. Recent challenges to realize the potential of this species in the agriculture industry are also discussed. Limited genetic resource for swamp buffalo has called for more genomics work to be done on this species including decoding its genome. As the economy progresses and farm mechanization increases, research and development for swamp buffaloes are focused on enhancing its productivity through understanding the genetics of agriculturally important traits. The use of genomic markers is a powerful tool to efficiently utilize the potential of this animal for food security and animal conservation. Understanding its genetics and retaining and maximizing its adaptability to harsher environments are a strategic move for food security in poorer nations in Southeast Asia in the face of climate change.

Bta-miR-376a Targeting KLF15 Interferes with Adipogenesis Signaling Pathway to Promote Differentiation of Qinchuan Beef Cattle Preadipocytes

Intramuscular fat (IMF) is a quality index associated with the taste and juiciness of meat. The deposition of IMF is affected by genetic and non-genetic factors, such as age, slaughter location, gender of the animal, and diet. Micro-ribonucleic acids (miRNA) are transcriptional regulators involved in adipogenesis, but the specific role of miR-376a in regulation of bovine adipocytes remains unknown. Our findings indicated that miR-376a was a potential negative regulator of bovine adipocyte differentiation. A bta-miR-376a mimic inhibited mRNA and protein expression of the marker genes, CDK1, CDK2, PCNA, C/EBPα, FAS, and PPAR γ, and significantly reduced ratios (%) of S-phase cells, the number of cells stained with 5-ethynyl-2'-deoxyuridine, and adipocyte proliferation. Oil red O staining and triglyceride content analysis also confirmed that bta-miR-376a was involved in adipocyte differentiation. Luciferase activities confirmed that Krüppel-like transcription factor 15 (KLF15) was a direct target gene of bta-miR-376a, and that KLF15 was a key transcription factor in adipogenesis. Therefore, bta-miR-376a might be a target for increasing beef IMF.

The effect of Acot2 overexpression or downregulation on the preadipocyte differentiation in Chinese Red Steppe cattle

The quality and nutritional value of beef is closely linked to its content of intramuscular fat (IMF). The differentiation of preadipocytes and the deposition of lipid droplets in the adipocytes are the key to regulate the IMF content. The differentiation of adipocytes is regulated by a series of transcription factors and genes. Acyl-CoA thioesterase 2 (Acot2) hydrolyzes the acyl-coenzyme A (CoA) into free fatty acids and CoA and has the potential to maintain the free fatty acids and acyl CoA at the cellular level. In this work, we detected the expression of the Acot2 gene during the adipocyte differentiation in Chinese Red Steppe cattle, and then interfered and overexpressed the Acot2 gene in the preadipocytes to explore its regulatory role in the adipocyte differentiation. The results showed that the expression and regulation of Acot2 mainly occurred at the later stage of the adipocyte differentiation. The interference with the Acot2 gene significantly inhibited the lipid droplets accumulation and triglyceride content, while its overexpression significantly promoted both of them. The results of this study show that the Acot2 gene is a positive regulator of the adipocyte differentiation and may become a new target to improve the beef quality.

Transcriptional regulation of adipogenic marker genes for the improvement of intramuscular fat in Qinchuan beef cattle

The intramuscular fat content plays a crucial role in meat quality traits. Increasing the degree of adipogenesis in beef cattle leads to an increase in the content of intramuscular fat. Adipogenesis a complex biochemical process which is under firm genetic control. Over the last three decades, the Qinchuan beef cattle have been extensively studied for the improvement of meat production and quality traits. In this study, we reviewed the literature regarding adipogenesis and intramuscular fat deposition. Then, we summarized the research conducted on the transcriptional regulation of key adipogenic marker genes, and also reviewed the roles of adipogenic marker genes in adipogenesis of Qinchuan beef cattle. This review will elaborate our understanding regarding transcriptional regulation which is a vital physiological process regulated by a cascade of transcription factors (TFs), key target marker genes, and regulatory proteins. This synergistic action of TFs and target genes ensures the accurate and diverse transmission of the genetic information for the accomplishment of central physiological processes. This information will provide an insight into the transcriptional regulation of the adipogenic marker genes and its role in bovine adipogenesis for the breed improvement programs especially for the trait of intramuscular fat deposition.

Genetic variants in MYF5 affected growth traits and beef quality traits in Chinese Qinchuan cattle

Myogenic factor 5 plays actively roles in the regulation of myogenesis. The aims of this study are to identify the evolution information of MYF5 protein among 10 domestic and mammalian animals, to uncover the expression patterns of MYF5 gene in calves and adults of Qinchuan cattle, and to expose the genetic variants of the MYF5 gene and explore its effect on cattle growth traits and beef quality traits in Qinchuan cattle. The bioinformatics results showed that the MYF5 proteins highly conserved in different mammalian or domestic animals apart from chicken. The expression level of MYF5 gene in the heart, muscle, lung, large intestine and liver was greater than that of other tissues. PCR amplicons sequencing identified four novel SNPs at g.5738A>G, g.5785C>T and g.5816A>G in the 3rd exon region and g.6535A>G in the 3' UTR. Genotypic frequencies of g.5785C>T was harshly deviated from the HWE (P < .05). Genetic diversity was low or intermediate for the four SNPs and those SNPs were in the weak linkage disequilibrium. Association analysis results indicated g.5785C>T, g.5816A>G and g.6535A>G significant effect on growth performance and beef quality traits of Qinchuan cattle. H1H3 diplotype had greater body size and better beef quality. All the results implicate that the MYF5 gene might be applied as a promising candidate gene in Qinchuan cattle breeding.

Genomic analyses reveal distinct genetic architectures and selective pressures in buffaloes

BACKGROUND

The domestic buffalo (Bubalus bubalis) is an essential farm animal in tropical and subtropical regions, whose genomic diversity is yet to be fully discovered.

RESULTS

In this study, we describe the demographic events and selective pressures of buffalo by analyzing 121 whole genomes (98 newly reported) from 25 swamp and river buffalo breeds. Both uniparental and biparental markers were investigated to provide the final scenario. The ancestors of swamp and river buffalo diverged ∼0.23 million years ago and then experienced independent demographic histories. They were domesticated in different regions, the swamp buffalo at the border between southwest China and southeast Asia, while the river buffalo in south Asia. The domestic stocks migrated to other regions and further differentiated, as testified by (at least) 2 ancestral components identified in each subspecies. Different signals of selective pressures were also detected in these 2 types of buffalo. The swamp buffalo, historically used as a draft animal, shows selection signatures in genes associated with the nervous system, while in river dairy breeds, genes under selection are related to heat stress and immunity.

CONCLUSIONS

Our findings substantially expand the catalogue of genetic variants in buffalo and reveal new insights into the evolutionary history and distinct selective pressures in river and swamp buffalo.

Function and Transcriptional Regulation of Bovine TORC2 Gene in Adipocytes: Roles of C/EBP, XBP1, INSM1 and ZNF263

The TORC2 gene is a member of the transducer of the regulated cyclic adenosine monophosphate (cAMP) response element binding protein gene family, which plays a key role in metabolism and adipogenesis. In the present study, we confirmed the role of TORC2 in bovine preadipocyte proliferation through cell cycle staining flow cytometry, cell counting assay, 5-ethynyl-2′-deoxyuridine staining (EdU), and mRNA and protein expression analysis of proliferation-related marker genes. In addition, Oil red O staining analysis, immunofluorescence of adiponectin, mRNA and protein level expression of lipid related marker genes confirmed the role of TORC2 in the regulation of bovine adipocyte differentiation. Furthermore, the transcription start site and sub-cellular localization of the TORC2 gene was identified in bovine adipocytes. To investigate the underlying regulatory mechanism of the bovine TORC2, we cloned a 1990 bp of the 5’ untranslated region (5′UTR) promoter region into a luciferase reporter vector and seven vector fragments were constructed through serial deletion of the 5′UTR flanking region. The core promoter region of the TORC2 gene was identified at location −314 to −69 bp upstream of the transcription start site. Based on the results of the transcriptional activities of the promoter vector fragments, luciferase activities of mutated fragments and siRNAs interference, four transcription factors (CCAAT/enhancer-binding protein C/BEPγ, X-box binding protein 1 XBP1, Insulinoma-associated 1 INSM1, and Zinc finger protein 263 ZNF263) were identified as the transcriptional regulators of TORC2 gene. These findings were further confirmed through Electrophoretic Mobility Shift Assay (EMSA) within nuclear extracts of bovine adipocytes. Furthermore, we also identified that C/EBPγ, XBP1, INSM1 and ZNF263 regulate TORC2 gene as activators in the promoter region. We can conclude that TORC2 gene is potentially a positive regulator of adipogenesis. These findings will not only provide an insight for the improvement of intramuscular fat in cattle, but will enhance our understanding regarding therapeutic intervention of metabolic syndrome and obesity in public health as well.

Identification of genetic markers associated with milk production traits in Chinese Holstein cattle based on post genome-wide association studies

With the rapid development of dairy industry, the breeding process of dairy cows has been accelerated. In previous genome-wide association studies (GWAS), a large number of genetic markers have been reported which may contribute to the selection of Holstein populations with superior milk-producing traits, but they remain to be further verified before practical application. In this study, 90 single nucleotide polymorphisms (SNPs) were selected, which were reported to be significantly associated with five milk production traits, including 305-day milk yield (305MY), 305-day milk fat percent (305FC), 305-day milk protein percent (305PC), 305-day milk fat yield (305FY) and 305-day milk protein yield (305PY). Effective 305-day data and fresh DNA samples were obtained from 295 healthy cows with gestational age of 1-4. Matrix-assisted laser desorption/ionization-time of flight (MALDI-TOF) was used to perform precise genotyping of these loci, followed by site association and haplotype analysis. Results showed that 36 out of 90 loci were supported to be used as genetic markers. In particular, several novel and effective haplotypes were also presented. Overall, our results verified tens of useful markers and provided a basis for further development of breeding strategies.