A marker-derived gene network reveals the regulatory role of PPARGC1A, HNF4G, and FOXP3 in intramuscular fat deposition of beef cattle

Regulatory impacts of PPARGC1A gene expression on milk production and cellular metabolism in buffalo mammary epithelial cells

The PPARGC1A gene plays a fundamental role in regulating cellular energy metabolism, including adaptive thermogenesis, mitochondrial biogenesis, adipogenesis, gluconeogenesis, and glucose/fatty acid metabolism. In a previous study, our group investigated seven SNPs in Mediterranean buffalo associated with milk production traits, and the current study builds on this research by exploring the regulatory influences of the PPARGC1A gene in buffalo mammary epithelial cells (BuMECs). Our findings revealed that knockdown of PPARGC1A gene expression significantly affected the growth of BuMECs, including proliferation, cell cycle, and apoptosis. Additionally, we observed downregulated triglyceride secretion after PPARGC1A knockdown. Furthermore, the critical genes related to milk production, including the STATS, BAD, P53, SREBF1, and XDH genes were upregulated after RNAi, while the FABP3 gene, was downregulated. Moreover, Silencing the PPARGC1A gene led to a significant downregulation of β-casein synthesis in BuMECs. Our study provides evidence of the importance of the PPARGC1A gene in regulating cell growth, lipid, and protein metabolism in the buffalo mammary gland. In light of our previous research, the current study underscores the potential of this gene for improving milk production efficiency and overall dairy productivity in buffalo populations.

An association weight matrix identified biological pathways associated with bull fertility traits in a multi-breed population

Using seven indicator traits, we investigated the genetic basis of bull fertility and predicted gene interactions from SNP associations. We used percent normal sperm as the key phenotype for the association weight matrix-partial correlation information theory (AWM-PCIT) approach. Beyond a simple list of candidate genes, AWM-PCIT predicts significant gene interactions and associations for the selected traits. These interactions formed a network of 537 genes: 38 genes were transcription cofactors, and 41 genes were transcription factors. The network displayed two distinct clusters, one with 294 genes and another with 243 genes. The network is enriched in fertility-associated pathways: steroid biosynthesis, p53 signalling, and the pentose phosphate pathway. Enrichment analysis also highlighted gene ontology terms associated with 'regulation of neurotransmitter secretion' and 'chromatin formation'. Our network recapitulates some genes previously implicated in another network built with lower-density genotypes. Sequence-level data also highlights additional candidate genes relevant to bull fertility, such as FOXO4, FOXP3, GATA1, CYP27B1, and EBP. A trio of regulatory genes-KDM5C, LRRK2, and PME-was deemed core to the network because of their overarching connections. This trio probably influences bull fertility through their interaction with genes, both known and unknown as to their role in male fertility. Future studies may target the trio and their target genes to enrich our understanding of male fertility further.

Population structure and identification of genomic regions associated with productive traits in five Italian beef cattle breeds

Italy has a long history in beef production, with local breeds such as Marchigiana, Chianina, Romagnola, Maremmana, and Podolica which produce high-quality meat. Selection has improved meat production, precocity, growth ability and muscle development, but the genetic determinism of such traits is mostly unknown. Using 33K SNPs-data from young bulls (N = 4064) belonging to these five Italian breeds, we demonstrated that the Maremmana and Podolica rustic breeds are closely related, while the specialised Marchigiana, Chianina, and Romagnola breeds are more differentiated. A genome-wide association study for growth and muscle development traits (average daily gain during the performance test, weight at 1 year old, muscularity) was conducted in the five Italian breeds. Results indicated a region on chromosome 2, containing the myostatin gene (MSTN), which displayed significant genome-wide associations with muscularity in Marchigiana cattle, a breed in which the muscle hypertrophy phenotype is segregating. Moreover, a significant SNP on chromosome 14 was associated, in the Chianina breed, to muscularity. The identification of diverse genomic regions associated with conformation traits might increase our knowledge about the genomic basis of such traits in Italian beef cattle and, eventually, such information could be used to implement marker-assisted selection of young bulls tested in the performance test.

Genome-wide scans identify biological and metabolic pathways regulating carcass and meat quality traits in beef cattle

Genome association studies (GWAS) provides knowledge about the genetic architecture of beef-related traits that allow linking the target phenotype to genomic information aiding breeding decision. Thus, the present study aims to uncover the genetic mechanism involved in carcass (REA: rib eye area, BF: backfat thickness, and HCW: hot carcass weight) and meat quality traits (SF: shear-force, MARB: marbling score, and IMF: intramuscular fat content) in Nellore cattle. For this, 6910 young bulls with phenotypic information and 23,859 animals genotyped with 435 k markers were used to perform the weighted single-step GBLUP (WssGBLUP) approach, considering two iterations. The top 10 genomic regions explained 8.13, 11.81, and 9.58% of the additive genetic variance, harboring a total of 119, 143, and 95 positional candidate genes for REA, BF, and HCW, respectively. For meat quality traits, the top 10 windows explained a large proportion of the total genetic variance for SF (14.95%), MARB (17.56%), and IMF (21.41%) surrounding 92, 155, and 111 candidate genes, respectively. Relevant candidate genes (CAST, PLAG1, XKR4, PLAGL2, AQP3/AQP7, MYLK2, WWOX, CARTPT, and PLA2G16) are related to physiological aspects affecting growth, carcass, meat quality, feed intake, and reproductive traits by signaling pathways controlling muscle control, key signal metabolic molecules INS / IGF-1 pathway, lipid metabolism, and adipose tissue development. The GWAS results provided insights into the genetic control of the traits studied and the genes found are potential candidates to be used in the improvement of carcass and meat quality traits.

RNA sequencing identifies key genes involved in intramuscular fat deposition in chickens at different developmental stages

BACKGROUND

Intramuscular fat (IMF) is an important factor in meat quality, and triglyceride (TG) and Phospholipids (PLIP), as the main components of IMF, are of great significance to the improvement of meat quality.

RESULTS

In this study, we used 30 RNA sequences generated from the transcriptome of chicken breast muscle tissues at different developmental stages to construct a gene expression matrix to map RNA sequence reads to the chicken genome and identify the transcript of origin. We used weighted gene co-expression network analysis (WGCNA) and identified 27 co-expression modules, 10 of which were related to TG and PLIP. We identified 150 highly-connected hub genes related to TG and PLIP, respectively, which were found to be mainly enriched in the adipocytokine signaling pathway, MAPK signaling pathway, mTOR signaling pathway, FoxO signaling pathway, and TGF-beta signaling pathway. Additionally, using the BioMart database, we identified 134 and 145 candidate genes related to fat development in the TG-related module and PLIP-related module, respectively. Among them, RPS6KB1, BRCA1, CDK1, RPS3, PPARGC1A, ACSL1, NDUFAB1, NDUFA9, ATP5B and PRKAG2 were identified as candidate genes related to fat development and highly-connected hub genes in the module, suggesting that these ten genes may be important candidate genes affecting IMF deposition.

CONCLUSIONS

RPS6KB1, BRCA1, CDK1, RPS3, PPARGC1A, ACSL1, NDUFAB1, NDUFA9, ATP5B and PRKAG2 may be important candidate genes affecting IMF deposition. The purpose of this study was to identify the co-expressed gene modules related to chicken IMF deposition using WGCNA and determine key genes related to IMF deposition, so as to lay a foundation for further research on the molecular regulation mechanism underlying chicken fat deposition.

Multitrait meta-analyses identify potential candidate genes for growth-related traits in Holstein heifers

Understanding the underlying pleiotropic relationships among growth and body size traits is important for refining breeding strategies in dairy cattle for optimal body size and growth rate. Therefore, we performed single-trait GWAS for monthly-recorded body weight (BW), hip height, body length, and chest girth from birth to 12 mo of age in Holstein animals, followed by stepwise multiple regression of independent or lowly-linked markers from GWAS loci using conditional and joint association analyses (COJO). Subsequently, we conducted a multitrait meta-analysis to detect pleiotropic markers. Based on the single-trait GWAS, we identified 170 significant SNPs, in which 59 of them remained significant after the COJO analyses. The most significant SNP, located at BTA7:3,676,741, explained 2.93% of the total phenotypic variance for BW6 (BW at 6 mo of age). We identified 17 SNPs with potential pleiotropic effects based on the multitrait meta-analyses, which resulted in 3 additional SNPs in comparison to those detected based on the single-trait GWAS. The identified quantitative trait loci regions overlap with genes known to influence human growth-related traits. According to positional and functional analyses, we proposed HMGA2, HNF4G, MED13L, BHLHE40, FRZB, DMP1, TRIB3, and GATAD2A as important candidate genes influencing the studied traits. The combination of single-trait GWAS and meta-analyses of GWAS results improved the efficiency of detecting associated SNPs, and provided new insights into the genetic mechanisms of growth and development in Holstein cattle.

IMF deposition ceRNA network analysis and functional study of HIF1a in yak

The concentration of intramuscular fat (IMF) is a crucial determinant of yak meat quality. However, the molecular mechanisms that regulate IMF in yak remain largely elusive. In our study, we conducted transcriptome sequencing on the longissimus dorsi muscle tissues of yaks with varying IMF contents. We then filtered differentially expressed genes (DEGs), microRNAs (DEMs), and long non-coding RNAs (DELs) to elucidate potential regulatory pathways of adipogenesis in yaks. Overall, our research sheds light on an array of potential mRNAs and noncoding RNAs implicated in IMF deposition and elaborates on the role of HIF1α in yaks. These findings contribute valuable insights that can serve as a guide for further research into the molecular mechanisms governing IMF deposition.

Genome-wide population structure, homozygosity, and heterozygosity patterns of Nero Siciliano pig in the framework of Italian and cosmopolitan breeds

Analysis of genomic data is becoming more and more common for the effective management of livestock breeding programmes, even in the case of local populations. In this work, the genome-wide data of Nero Siciliano pig breed were compared to that of wild boar, Italian local and cosmopolitan breeds to investigate its genetic structure, and runs of homozygosity (ROH) and heterozygosity patterns. The Nero Siciliano has been reported to have the highest rate of genetic diversity among the Italian breeds, and a genetic variability comparable to that of the cosmopolitan breeds. Analyses of genomic structure and relationships underlined its proximity to wild boar, and an internal substructure probably linked to different family lines. The breed showed a low value of inbreeding estimated from ROH, and the highest diversity index among the Italian breeds, even if lower than that of the cosmopolitans. Four ROH islands in three chromosomes (SSC8, SSC11, and SSC14) and one heterozygosity-rich region (SSC1) were identified in Nero Siciliano, highlighting genomic regions related to productive QTL. Across breeds, SSC8 and SSC14 were the chromosomes with most ROH islands, with Mora Romagnola and wild boar showing the highest level of autozygosity. Chromosomes SSC2, SSC6, SSC8 and SSC13 showed the majority of runs of heterozygosity regions, mainly found in the cosmopolitan pig breeds, which reported several genes associated with health-related QTL. The outlined results can help to better identify the genomic profile of this local breed in order to plan matings, maintain adequate internal diversity and exploit the production system.

Candidate genes for domestication and resistance to cold climate according to whole genome sequencing data of Russian cattle and sheep breeds

It is known that different species of animals, when living in the same environmental conditions, can form similar phenotypes. The study of the convergent evolution of several species under the influence of the same environmental factor makes it possible to identify common mechanisms of genetic adaptation. Local cattle and sheep breeds have been formed over thousands of years under the influence of domestication, as well as selection aimed at adaptation to the local environment and meeting human needs. Previously, we identified a number of candidate genes in genome regions potentially selected during domestication and adaptation to the climatic conditions of Russia, in local breeds of cattle and sheep using whole genome genotyping data. However, these data are of low resolution and do not reveal most nucleotide substitutions. The aim of the work was to create, using the whole genome sequencing data, a list of genes associated with domestication, selection and adaptation in Russian cattle and sheep breeds, as well as to identify candidate genes and metabolic pathways for selection for cold adaptation. We used our original data on the search for signatures of selection in the genomes of Russian cattle (Yakut, Kholmogory, Buryat, Wagyu) and sheep (Baikal, Tuva) breeds. We used the HapFLK, DCMS, FST and PBS methods to identify DNA regions with signatures of selection. The number of candidate genes in potentially selective regions was 946 in cattle and 151 in sheep. We showed that the studied Russian cattle and sheep breeds have at least 10 genes in common, apparently involved in the processes of adaptation/selection, including adaptation to a cold climate, including the ASTN2, PM20D1, TMEM176A, and GLIS1 genes. Based on the intersection with the list of selected genes in at least two Arctic/Antarctic mammal species, 20 and 8 genes, have been identified in cattle and sheep, respectively, that are potentially involved in cold adaptation. Among them, the most promising for further research are the ASPH, NCKAP5L, SERPINF1, and SND1 genes. Gene ontology analysis indicated the existence of possible common biochemical pathways for adaptation to cold in domestic and wild mammals associated with cytoskeleton disassembly and apoptosis.

Integrated Comparative Transcriptome and circRNA-lncRNA-miRNA-mRNA ceRNA Regulatory Network Analyses Identify Molecular Mechanisms Associated with Intramuscular Fat Content in Beef Cattle

Intramuscular fat content (IMF), one of the most important carcass traits in beef cattle, is controlled by complex regulatory factors. At present, molecular mechanisms involved in regulating IMF and fat metabolism in beef cattle are not well understood. Our objective was to integrate comparative transcriptomic and competing endogenous RNA (ceRNA) network analyses to identify candidate messenger RNAs (mRNAs) and regulatory RNAs involved in molecular regulation of longissimus dorsi muscle (LDM) tissue for IMF and fat metabolism of 5 beef cattle breeds (Angus, Chinese Simmental, Luxi, Nanyang, and Shandong Black). In total, 34 circRNAs, 57 lncRNAs, 15 miRNAs, and 374 mRNAs were identified by integrating gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses. Furthermore, 7 key subnets with 16 circRNAs, 43 lncRNAs, 7 miRNAs, and 237 mRNAs were detected through clustering analyses, whereas GO enrichment analysis of identified RNAs revealed 48, 13, and 28 significantly enriched GO terms related to IMF in biological process, molecular function, and cellular component categories, respectively. The main metabolic-signaling pathways associated with IMF and fat metabolism that were enriched included metabolic, calcium, cGMP-PKG, thyroid hormone, and oxytocin signaling pathways. Moreover, MCU, CYB5R1, and BAG3 genes were common among the 10 comparative groups defined as important candidate marker genes for fat metabolism in beef cattle. Contributions of transcriptome profiles from various beef breeds and a competing endogenous RNA (ceRNA) regulatory network underlying phenotypic differences in IMF provided novel insights into molecular mechanisms associated with meat quality.

Genome-wide association analysis of milk production, somatic cell score, and body conformation traits in Holstein cows

Milk production and body conformation traits are critical economic traits for dairy cows. To understand the basic genetic structure for those traits, a genome wide association study was performed on milk yield, milk fat yield, milk fat percentage, milk protein yield, milk protein percentage, somatic cell score, body form composite index, daily capacity composite index, feed, and leg conformation traits, based on the Illumina Bovine HD100k BeadChip. A total of 57, 12 and 26 SNPs were found to be related to the milk production, somatic cell score and body conformation traits in the Holstein cattle. Genes with pleiotropic effect were also found in this study. Seven significant SNPs were associated with multi-traits and were located on the PLEC, PLEKHA5, TONSL, PTGER4, and LCORL genes. In addition, some important candidate genes, like GPAT3, CEBPB, AGO2, SLC37A1, and FNDC3B, were found to participate in fat metabolism or mammary gland development. These results can be used as candidate genes for milk production, somatic cell score, and body conformation traits of Holstein cows, and are helpful for further gene function analysis to improve milk production and quality.

Insights into the architecture of human-induced polygenic selection in Duroc pigs

Background

As one of the most utilized commercial composite boar lines, Duroc pigs have been introduced to China and undergone strongly human-induced selection over the past decades. However, the efficiencies and limitations of previous breeding of Chinese Duroc pigs are largely understudied. The objective of this study was to uncover directional polygenic selection in the Duroc pig genome, and investigate points overlooked in the past breeding process.

Results

Here, we utilized the Generation Proxy Selection Mapping (GPSM) on a dataset of 1067 Duroc pigs with 8,766,074 imputed SNPs. GPSM detected a total of 5649 putative SNPs actively under selection in the Chinese Duroc pig population, and the potential functions of the selection regions were mainly related to production, meat and carcass traits. Meanwhile, we observed that the allele frequency of variants related to teat number (NT) relevant traits was also changed, which might be influenced by genes that had pleiotropic effects. First, we identified the direction of selection on NT traits by \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hat{G}$$\end{document}G^, and further pinpointed large-effect genomic regions associated with NT relevant traits by selection signature and GWAS. Combining results of NT relevant traits-specific selection signatures and GWAS, we found three common genome regions, which were overlapped with QTLs related to production, meat and carcass traits besides “teat number” QTLs. This implied that there were some pleiotropic variants underlying NT and economic traits. We further found that rs346331089 has pleiotropic effects on NT and economic traits, e.g., litter size at weaning (LSW), litter weight at weaning (LWW), days to 100 kg (D100), backfat thickness at 100 kg (B100), and loin muscle area at 100 kg (L100) traits.

Conclusions

The selected loci that we identified across methods displayed the past breeding process of Chinese Duroc pigs, and our findings could be used to inform future breeding decision.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40104-022-00751-x.

Use of a graph neural network to the weighted gene co-expression network analysis of Korean native cattle

In the general framework of the weighted gene co-expression network analysis (WGCNA), a hierarchical clustering algorithm is commonly used to module definition. However, hierarchical clustering depends strongly on the topological overlap measure. In other words, this algorithm may assign two genes with low topological overlap to different modules even though their expression patterns are similar. Here, a novel gene module clustering algorithm for WGCNA is proposed. We develop a gene module clustering network (gmcNet), which simultaneously addresses single-level expression and topological overlap measure. The proposed gmcNet includes a “co-expression pattern recognizer” (CEPR) and “module classifier”. The CEPR incorporates expression features of single genes into the topological features of co-expressed ones. Given this CEPR-embedded feature, the module classifier computes module assignment probabilities. We validated gmcNet performance using 4,976 genes from 20 native Korean cattle. We observed that the CEPR generates more robust features than single-level expression or topological overlap measure. Given the CEPR-embedded feature, gmcNet achieved the best performance in terms of modularity (0.261) and the differentially expressed signal (27.739) compared with other clustering methods tested. Furthermore, gmcNet detected some interesting biological functionalities for carcass weight, backfat thickness, intramuscular fat, and beef tenderness of Korean native cattle. Therefore, gmcNet is a useful framework for WGCNA module clustering.

Within- and between-Breed Selection Signatures in the Original and Improved Valachian Sheep

This study explored the genomic diversity and selection signatures in two Slovakian national breeds, the Original Valachian and the Improved Valachian sheep. As they are an important animal genetic resource within the country, but with decreasing population size, our aim is to identify potentially valuable genomic regions. A total of 97 sheep (18 male and 79 female) from the Original Valachian, and 69 sheep (25 male and 44 female) from the Improved Valachian populations were genotyped using the GeneSeek GGP Ovine 50 K chip. The inbreeding levels were assessed with runs of homozygosity (ROH). The selection signatures within breeds were identified based on the top 1% of most homozygous regions within the breed, the so-called ROH islands. The selection signatures between breeds were assessed based on variance in linkage disequilibrium. Overall, we have identified selection signatures with quantitative trait loci (QTL) and genes pointing towards all three production purposes of the Valachian sheep, milk, meat, and wool, including their quality characteristics. Another group with apparent large importance was the various traits related to health and resistance to parasites, which is well in line with the sturdy nature of this breed.

Identifying Key Genes and Functionally Enriched Pathways of Diverse Adipose Tissue Types in Cattle

Background: Fat is a tissue that not just stores energy and plays a protective role; it is also a vital endocrine organ that generates and integrates signals to influence metabolism. Meanwhile, the excessive accumulation of lipids in adipose tissue can lead to metabolic disturbance and diseases. To date, the complicated molecular mechanisms of bovine adipose tissue are still unknown. This study aimed to identify key genes and functionally enriched pathways in various adipose tissue types.

Results: The RNAseq data of 264 samples were downloaded from Gene Expression Omnibus (GEO) and analyzed by weighted gene co-expression network analysis (WGCNA). We identified 19 modules that significantly associated with at least one adipose tissue type. The brown module from GSE39618 was most closely associated with intramuscular fat tissue, which contained 550 genes. These genes were significantly enriched in pathways that related to inflammation and disease, such as TNF signaling pathway, IL-17 signaling pathway, and NF-kappa B signaling pathway. The pink module (GSE39618) that contained 58 genes was most closely associated with omental fat tissue. The turquoise (GSE39618), blue (GSE116775), and yellow (GSE65125) module were most closely associated with subcutaneous fat tissue. Genes in these modules were significantly enriched in pathways related to fat metabolism, such as the PPAR signaling pathway, fatty acid metabolism and PI3K-Akt signaling pathway. At last, key genes for intramuscular fat (PTGS2 and IL6), omental fat (ARHGEF5 and WT1), and subcutaneous fat (KIT, QR6Q1, PKD2L1, etc.) were obtained and verified. In addition, it was found that IL10 and VCAM1 might be potential genes to distinguish adipose and muscle.

Conclusion: The study applied WGCNA to generate a landscape of adipose tissue and provide a basis for identifying potential pathways and hub genes of different adipose tissue types.

Forensic DNA phenotyping: Canis familiaris breed classification and skeletal phenotype prediction using functionally significant skeletal SNPs and indels

This review highlights a novel application of breed identification and prediction of skeletal traits in forensic investigations using canine DNA evidence. Currently, genotyping methods used for canine breed classification involve the application of highly polymorphic short tandem repeats in addition to larger commercially available SNP arrays. Both applications face technical challenges. An additional approach to breed identification could be through genotyping SNPs and indels that characterise the array of skeletal differences displayed across domestic dog populations. Research has shown that a small number of genetic variants of large effect drive differences in skeletal phenotypes among domestic dog breeds. This feature makes functionally significant canine skeletal variants a cost-effective target for forensic investigators to classify individuals according to their breed. Further analysis of these skeletal variants would enable the prediction of external appearance. To date, functionally significant genes with genetic variants associated with differences in size, bulk, skull shape, ear shape, limb length, digit type, and tail morphology have been uncovered. Recommendations of a cost-effective genotyping method that can be readily designed and applied by forensic investigators have been given. Further advances to improve the field of canine skeletal forensic DNA phenotyping include the refinement of phenotyping methods, further biological validation of the skeletal genetic variants and establishing a publicly available database for storage of allele frequencies of the skeletal genetic variants in the wider domestic dog population.

Chromatin accessibility and regulatory vocabulary across indicine cattle tissues

BACKGROUND

Spatiotemporal changes in the chromatin accessibility landscape are essential to cell differentiation, development, health, and disease. The quest of identifying regulatory elements in open chromatin regions across different tissues and developmental stages is led by large international collaborative efforts mostly focusing on model organisms, such as ENCODE. Recently, the Functional Annotation of Animal Genomes (FAANG) has been established to unravel the regulatory elements in non-model organisms, including cattle. Now, we can transition from prediction to validation by experimentally identifying the regulatory elements in tropical indicine cattle. The identification of regulatory elements, their annotation and comparison with the taurine counterpart, holds high promise to link regulatory regions to adaptability traits and improve animal productivity and welfare.

RESULTS

We generate open chromatin profiles for liver, muscle, and hypothalamus of indicine cattle through ATAC-seq. Using robust methods for motif discovery, motif enrichment and transcription factor binding sites, we identify potential master regulators of the epigenomic profile in these three tissues, namely HNF4, MEF2, and SOX factors, respectively. Integration with transcriptomic data allows us to confirm some of their target genes. Finally, by comparing our results with Bos taurus data we identify potential indicine-specific open chromatin regions and overlaps with indicine selective sweeps.

CONCLUSIONS

Our findings provide insights into the identification and analysis of regulatory elements in non-model organisms, the evolution of regulatory elements within two cattle subspecies as well as having an immediate impact on the animal genetics community in particular for a relevant productive species such as tropical cattle.

Nr2e1 deficiency aggravates insulin resistance and chronic inflammation of visceral adipose tissues in a diet-induced obese mice model

AIMS

To investigate the nuclear receptor subfamily 2 group E member 1 (Nr2e1) expression in adipose tissues of obese mice and assess the role of Nr2e1 in insulin resistance and chronic inflammation of the adipose tissues.

MAIN METHODS

An obese model was established in Nr2e1 knockout (KO) mice and their wild type (WT) littermates through a long-term high-fat diet (HFD) feeding regime. The epididymal fat weight, body weight, and daily food intake were recorded. The blood lipid profile, blood inflammatory factors, and the levels of fasting blood glucose (FBG) and fasting insulin were determined. We estimated insulin resistance by the homeostasis model assessment (HOMA). The expression of inflammatory factors and F4/80 was examined by polymerase chain reaction (PCR) and western blotting to assess adipose tissues inflammation. We also determined the molecules of insulin signaling and the nuclear factor kappa B (NF-κB) pathway by western blotting.

KEY FINDINGS

The Nr2e1 expression was upregulated in WT obese mice when compared with that in control mice. Despite a lower body weight and epididymal fat mass in Nr2e1^-/- mice, these rats showed increased inflammatory cytokines secretion, more pronounced hyperlipidemia, and impaired insulin sensitivity after HFD treatment. Further investigation revealed that Nr2e1 deletion affected the expression of insulin signaling and NF-κB pathway-related molecules in visceral adipose tissues.

SIGNIFICANCE

Nr2e1 may act as a potential target to improve insulin sensitivity and inflammation in obesity and related complications.

Identification and co-expression analysis of long noncoding RNAs and mRNAs involved in the deposition of intramuscular fat in Aohan fine-wool sheep

BACKGROUND

Intramuscular fat (IMF) content has become one of the most important indicators for measuring meat quality, and levels of IMF are affected by various genes. Long non-coding RNAs (lncRNAs) are widely expressed non-coding RNAs that play an important regulatory role in a variety of biological processes; however, research on the lncRNAs involved in sheep IMF deposition is still in its infancy. Aohan fine-wool sheep (AFWS), one of China's most important meat-hair, dual-purpose sheep breed, provides a great model for studying the role of lncRNAs in the regulation of IMF deposition. We identified lncRNAs by RNA sequencing in Longissimus thoracis et lumborum (LTL) samples of sheep at two ages: 2 months (Mth-2) and 12 months (Mth-12).

RESULTS

We identified a total of 26,247 genes and 6935 novel lncRNAs in LTL samples of sheep. Among these, 199 mRNAs and 61 lncRNAs were differentially expressed. We then compared the structural characteristics of lncRNAs and mRNAs. We obtained target genes of differentially expressed lncRNAs (DELs) and performed enrichment analyses using Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG). We found that target mRNAs were enriched in metabolic processes and developmental pathways. One pathway was significantly enriched, namely tight junction. Based on the analysis of critical target genes, we obtained seven candidate lncRNAs that potentially regulated lipid deposition and constructed a lncRNA-mRNA co-expression network that included MSTRG.4051.3-FZD4, MSTRG.16157.3-ULK1, MSTRG.21053.3-PAQR3, MSTRG.19941.2-TPI1, MSTRG.12864.1-FHL1, MSTRG.2469.2-EXOC6 and MSTRG.21381.1-NCOA1. We speculated that these candidate lncRNAs might play a role by regulating the expression of target genes. We randomly selected five mRNAs and five lncRNAs to verify the accuracy of the sequencing data by qRT-PCR.

CONCLUSIONS

Our study identified the differentially expressed mRNAs and lncRNAs during intramuscular lipid deposition in Aohan fine-wool sheep. The work may widen the knowledge about the annotation of the sheep genome and provide a working basis for investigating intramuscular fat deposition in sheep.

Genome-Wide Association Analysis of Growth Curve Parameters in Chinese Simmental Beef Cattle

The objective of the present study was to perform a genome-wide association study (GWAS) for growth curve parameters using nonlinear models that fit original weight-age records. In this study, data from 808 Chinese Simmental beef cattle that were weighed at 0, 6, 12, and 18 months of age were used to fit the growth curve. The Gompertz model showed the highest coefficient of determination (R² = 0.954). The parameters' mature body weight (A), time-scale parameter (b), and maturity rate (K) were treated as phenotypes for single-trait GWAS and multi-trait GWAS. In total, 9, 49, and 7 significant SNPs associated with A, b, and K were identified by single-trait GWAS; 22 significant single nucleotide polymorphisms (SNPs) were identified by multi-trait GWAS. Among them, we observed several candidate genes, including PLIN3, KCNS3, TMCO1, PRKAG3, ANGPTL2, IGF-1, SHISA9, and STK3, which were previously reported to associate with growth and development. Further research for these candidate genes may be useful for exploring the full genetic architecture underlying growth and development traits in livestock.

Comparative transcriptome analysis reveals that PCK1 is a potential gene affecting IMF deposition in buffalo

BACKGROUND

In China, although buffaloes are abundant, beef is mainly obtained from cattle, and this preference is mainly attributed to the low intramuscular fat (IMF) content of buffalo. Genetic factors are an important driver that affects IMF deposition.

RESULTS

To reveal the intrinsic factors responsible for the low IMF content of buffalo, mRNA expression patterns in muscle and adipose tissue between buffalo and cattle were characterized by RNA sequencing analysis. The IMF content in Nanyang cattle was higher than that in Xinyang buffalo. A total of 1566 mRNAs expressed in adipose tissue showed differential expression between the longissimus dorsi muscles of buffalo and cattle. Functional annotation suggested a difference in the glycolysis/gluconeogenesis pathway between the two species. The results of RT-qPCR analysis and gain-of-function experiments confirmed the positive association between the IMF content and phosphoenolpyruvate carboxykinase 1 (PCK1) expression in buffalo. In both mouse C2C12 cells and cultured bovine myocytes, the activity of the PCK1 promoter in buffalo is lower than that in cattle. However, in mouse 3T3-L1 adipocytes and cultured bovine adipocytes, the activity of PCK1 in buffalo promoter is higher than that in cattle.

CONCLUSIONS

These results indicate the important role of PCK1 in buffalo IMF deposition and illustrate the differences between buffalo and cattle promoter activity that drive PCK1 expression. This research helps to establish a foundation for further studies investigating IMF deposition in buffalo.

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

Simple Summary

Hanwoo is an indigenous cattle breed in Korea and popular for meat production owing to its rapid growth and high-quality meat. Its yearling weight and carcass traits (backfat thickness, carcass weight, eye muscle area, and marbling score) are economically important for the selection of young and proven bulls. In recent decades, the advent of high throughput genotyping technologies has made it possible to perform genome-wide association studies (GWAS) for the detection of genomic regions associated with traits of economic interest in different species. In this study, we conducted a weighted single-step genome-wide association study which combines all genotypes, phenotypes and pedigree data in one step (ssGBLUP). It allows for the use of all SNPs simultaneously along with all phenotypes from genotyped and ungenotyped animals. Our results revealed 33 relevant genomic regions related to the traits of interest. Gene set enrichment analysis indicated that the identified candidate genes were related to biological processes and functional terms that were involved in growth and lipid metabolism. In conclusion, these results suggest that the incorporation of GWAS results and network analysis can help us to better understand the genetic bases underlying growth and carcass traits.

Abstract

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

Genome-wide association studies reveal novel loci associated with carcass and body measures in beef cattle

Genomic selection has an essential role in the livestock economy by increasing selection productivity. Genomics provides a mechanism to increase the rate of genetic gain using marker-assisted selection. Various quantitative trait loci (QTL) associated with body, carcass and meat quality traits in beef cattle have been found. It is widely accepted that QTL traits in livestock species are regulated by several genes and factors from the environment. Genome-wide association studies (GWAS) are a powerful approach in identifying QTL and to establish genomic regions harboring the genes and polymorphisms associated with specific characteristics in beef cattle. Due to their impact on economic returns, growth, carcass and meat quality traits of cattle are frequently used as essential criteria in selection in breeding programs., GWAS has been used in beef cattle breeding and genetic program and some progress has been made. Furthermore, numerous genes and markers related to productivity traits in beef cattle have been found. This review summarizes the advances in the use of GWAS in beef cattle production and outlines the associations with growth, carcass, and meat quality.

Molecular characterization, tissue expression and polymorphisms of buffalo PPARGC1A gene

PPARGC1A exerts important functions in activating many nuclear receptors and transcription factors that are related to energy balance. Previous studies have shown that PPARGC1A gene is associated with lactation traits of dairy cattle. However, the functional role of the buffalo PPARGC1A gene is still unknown. In this work, the complete coding sequence (CDS) of buffalo PPARGC1A was isolated and characterized for swamp and river buffalo. The CDS length of PPARGC1A for both types of buffalo was the same, which was composed of 2394 nucleotides and encoded a peptide composed of 797 amino acid residues. This protein belonged to a hydrophilic protein and contained one RRM_PPARGC1A domain (AA 674–764) without a signal peptide or a transmembrane domain. The differential expressions of this gene in 10 buffalo tissues in lactation and non-lactation displayed that the PPARGC1A was highly expressed in the muscle, heart, liver, brain and kidney of both non-lactating and lactating periods, but its expression was significantly different in the muscle, heart, liver, small intestine, mammary gland, rumen, spleen and lung between the two periods. Eight single nucleotide polymorphisms (SNPs) were found in buffalo, in which the c.778C>T, c.1257G>A and c.1311G>A were shared by two types of buffalo with similar allele frequencies, while the c.419C>T, c.759A>G, c.920C>A, c.926G>A and c.1509A>T were only observed in river buffalo. The SNP419, SNP920 and SNP926 were non-synonymous, which led to the amino acid changes of p.Ser140Phe, p.Pro307His and p.Arg309Lys. Seven nucleotide differential sites were identified in the PPARGC1A gene between buffalo and other Bovidae species. Phylogenetic analysis indicated that buffaloes were independently clustered into one branch, but they were closely related to the species of the Bos genus. The results indicate that buffalo PPARGC1A is an inducible transcriptional coactivator involved in regulating carbohydrate and fat metabolism. It can exert a functional role in a variety of buffalo tissues and may participate in milk fat synthesis and development in the mammary gland.

Polymorphisms in 5' proximal regulating region of THRSP gene are associated with fat production in pigs

Chinese and imported pig breeds differ in fat production potential, which is associated with the polymorphisms in the 5' proximal regulating region (5'PRR) of thyroid hormone responsive gene (THRSP). In three Chinese breeds (Dingyuan, CDY; Wannanhua, CWH; and Jixi, CJX) and one introduced breed (Yorkshire, YKS), three variant sites were located at T/C-400, A/G-376, and G/A-98 in the 5'PRR. Chinese pig breeds had higher C-400 allele frequencies than YKS. The frequencies of A-376 in CDY and G-376 in CWH were about 0.8. G-98 allele frequencies in CWH and YKS were 0.8617 and 0.8149, respectively. TGG was the dominant haplotype in YKS, CGG in CWH and CJX, and CAA in CDY. According to haplotype frequency, four breeds were clustered into three types, which was consistent with the geographical distribution of the breeds. In CDY, the average backfat thickness (BFT) was the highest with the CC-400 genotype, followed by CT-400 and TT-400 genotypes. In YKS, the pigs with CC-400 or CT-400 genotypes had higher BFT and average daily weight gain, whereas those with CC-400 or TT-400 genotypes had larger lion-eye area. No significant difference was observed in carcass traits among different genotypes at the A/G-376 and G/A-98 loci. The mRNA abundance of THRSP expression for the CCAGAG genotype was significantly higher than that for CTAGAG or TTAGAG genotype. These results indicated that the polymorphisms and genotype distribution of THRSP were closely related to the potential for fat production in pig breeds, which were the result of adaptation to artificial selection and natural selection.

Use of gene expression profile to identify potentially relevant transcripts to myofibrillar fragmentation index trait

The myofibrillar fragmentation index (MFI) is an indicative trait for meat tenderness. Longissimus thoracis muscle samples from the 20 most extreme bulls (out of 80 bulls set) for MFI (high (n = 10) and low (n = 10) groups) trait were used to perform transcriptomic analysis, using RNA Sequencing (RNA-Seq). An average of 24.616 genes was expressed in the Nellore muscle transcriptome analysis. A total of 96 genes were differentially expressed (p value ≤ 0.001) between the two groups of divergent bulls for MFI. The HEBP2 and BDH1 genes were overexpressed in animals with high MFI. The MYBPH and MYL6, myosin encoders, were identified. The differentially expressed genes were related to increase mitochondria efficiency, especially in cells under oxidative stress conditions, and these also were related to zinc and calcium binding, membrane transport, and muscle constituent proteins, such as actin and myosin. Most of those genes were involved in metabolic pathways of oxidation-reduction, transport of lactate in the plasma membrane, and muscle contraction. This is the first study applying MFI phenotypes in transcriptomic studies to identify and understand differentially expressed genes for beef tenderness. These results suggest that differences detected in gene expression between high and low MFI animals are related to reactive mechanisms and structural components of oxidative fibers under the condition of cellular stress. Some genes may be selected as positional candidate genes to beef tenderness, MYL6, MYBPH, TRIM63, TRIM55, TRIOBP, and CHRNG genes. The use of MFI phenotypes could enhance results of meat tenderness studies.

Genetic Architecture of Carcass and Meat Quality Traits in Montana Tropical® Composite Beef Cattle

The Montana Tropical^® Composite is a recently developed beef cattle population that is rapidly expanding in Brazil and other tropical countries. This is mainly due to its improved meat quality and adaptation to tropical climate conditions compared to Zebu and Taurine cattle breeds, respectively. This study aimed to investigate the genetic architecture of ultrasound-based carcass and meat quality traits in Montana Tropical^® Composite beef cattle. Therefore, we estimated variance components and genetic parameters and performed genome-wide association studies using the weighted single-step Genomic Best Linear Unbiased Prediction (GBLUP) approach. A pedigree dataset containing 28,480 animals was used, in which 1,436 were genotyped using a moderate-density Single Nucleotide Polymorphism panel (30K; 30,105 SNPs). A total of 9,358, 5,768, 7,996, and 1,972 phenotypic records for the traits Longissimus muscle area (LMA), backfat thickness (BFT), rump fat thickness (RFT), and for marbling score (MARB), respectively, were used for the analyses. Moderate to high heritability estimates were obtained and ranged from 0.16 ± 0.03 (RFT) to 0.33 ± 0.05 (MARB). A high genetic correlation was observed between BFT and RFT (0.97 ± 0.02), suggesting that a similar set of genes affects both traits. The most relevant genomic regions associated with LMA, BFT, RFT, and MARB were found on BTA10 (5.4–5.8 Mb), BTA27 (25.2–25.5 Mb), BTA18 (60.6–61.0 Mb), and BTA21 (14.8–15.4 Mb). Two overlapping genomic regions were identified for RFT and MARB (BTA13:47.9–48.1 Mb) and for BFT and RFT (BTA13:61.5–62.3 Mb). Candidate genes identified in this study, including PLAG1, LYN, WWOX, and PLAGL2, were previously reported to be associated with growth, stature, skeletal muscle growth, fat thickness, and fatty acid composition. Our results indicate that ultrasound-based carcass and meat quality traits in the Montana Tropical^® Composite beef cattle are heritable, and therefore, can be improved through selective breeding. In addition, various novel and already known genomic regions related to these traits were identified, which contribute to a better understanding of the underlying genetic background of LMA, BFT, RFT, and MARB in the Montana Tropical Composite population.

Combining information from genome-wide association and multi-tissue gene expression studies to elucidate factors underlying genetic variation for residual feed intake in Australian Angus cattle

BACKGROUND

Genome-wide association studies (GWAS) are extensively used to identify single nucleotide polymorphisms (SNP) underlying the genetic variation of complex traits. However, much uncertainly often still exists about the causal variants and genes at quantitative trait loci (QTL). The aim of this study was to identify QTL associated with residual feed intake (RFI) and genes in these regions whose expression is also associated with this trait. Angus cattle (2190 steers) with RFI records were genotyped and imputed to high density arrays (770 K) and used for a GWAS approach to identify QTL associated with RFI. RNA sequences from 126 Angus divergently selected for RFI were analyzed to identify the genes whose expression was significantly associated this trait with special attention to those genes residing in the QTL regions.

RESULTS

The heritability for RFI estimated for this Angus population was 0.3. In a GWAS, we identified 78 SNPs associated with RFI on six QTL (on BTA1, BTA6, BTA14, BTA17, BTA20 and BTA26). The most significant SNP was found on chromosome BTA20 (rs42662073) and explained 4% of the genetic variance. The minor allele frequencies of significant SNPs ranged from 0.05 to 0.49. All regions, except on BTA17, showed a significant dominance effect. In 1 Mb windows surrounding the six significant QTL, we found 149 genes from which OAS2, STC2, SHOX, XKR4, and SGMS1 were the closest to the most significant QTL on BTA17, BTA20, BTA1, BTA14, and BTA26, respectively. In a 2 Mb windows around the six significant QTL, we identified 15 genes whose expression was significantly associated with RFI: BTA20) NEURL1B and CPEB4; BTA17) RITA1, CCDC42B, OAS2, RPL6, and ERP29; BTA26) A1CF, SGMS1, PAPSS2, and PTEN; BTA1) MFSD1 and RARRES1; BTA14) ATP6V1H and MRPL15.

CONCLUSIONS

Our results showed six QTL regions associated with RFI in a beef Angus population where five of these QTL contained genes that have expression associated with this trait. Therefore, here we show that integrating information from gene expression and GWAS studies can help to better understand the genetic mechanisms that determine variation in complex traits.

Association of Twelve Candidate Gene Polymorphisms with the Intramuscular Fat Content and Average Backfat Thickness of Chinese Suhuai Pigs

Simple Summary

Appropriate intramuscular fat content (IFC) is the goal of consumers and the direction that breeders must pursue. However, it is difficult to improve the IFC but not average backfat thickness (ABT) by traditional breeding methods, and pigs must be slaughtered to accurately measure IFC. Marker-assisted selection (MAS) provides an economic and efficient method to improve the IFC in pigs. Our research indicated that the FABP3 (rs1110770079) single nucleotide polymorphism (SNP) could be a candidate gene associated with IFC (but not ABT), and IFC could be improved by selecting the individuals with a favorable genotype (GG) of FABP3 (rs1110770079) SNP for pig breeding.

Abstract

The present study aimed to identify the molecular markers for genes that influence intramuscular fat content (IFC), but not average backfat thickness (ABT). A total of 330 Suhuai pigs were slaughtered, and measurements of IFC and ABT were obtained. Phenotypic and genetic correlations between IFC and ABT were calculated. Thirteen single nucleotide polymorphisms (SNPs) among 12 candidate genes for IFC were analyzed, including FABP3, LIPE, IGF1, IGF2, LEP, LEPR, MC4R, PHKG1, RETN, RYR1, SCD, and UBE3C. Associations of the evaluated SNPs with IFCIFC and ABT were performed. Our results showed that the means of IFC and ABT were 1.99 ± 0.03 % and 26.68 ± 0.28 mm, respectively. The coefficients of variation (CVs) of IFC and ABT were 31.21% and 19.36%, respectively. The phenotypic and genetic correlations between IFC and ABT were moderate. Only the FABP3 (rs1110770079) was associated with IFC (p < 0.05) but not with ABT. Besides, there was a tendency for associations of RYR1 (rs344435545) and SCD (rs80912566) with IFC (p < 0.1). Our results indicated that the FABP3 (rs1110770079) SNP could be used as a marker to improve IFC without changing ABT in the Suhuai pig breeding system.

Comparative Transcriptomic and Proteomic Analyses Identify Key Genes Associated With Milk Fat Traits in Chinese Holstein Cows

Milk fat is the most important energy substance in milk and contributes to its quality and health benefits. However, the genetic mechanisms underlying milk fat synthesis are not fully understood. The development of RNA sequencing and tandem mass tag technologies has facilitated the identification of eukaryotic genes associated with complex traits. In this study, we used these methods to obtain liver transcriptomic and proteomic profiles of Chinese Holstein cows (n = 6). Comparative analyses of cows with extremely high vs. low milk fat percentage phenotypes yielded 321 differentially expressed genes (DEGs) and 76 differentially expressed proteins (DEPs). Functional annotation of these DEGs and DEPs revealed 26 genes that were predicted to influence lipid metabolism through insulin, phosphatidylinositol 3-kinase/Akt, mitogen-activated protein kinase, 5′ AMP-activated protein kinase, mammalian target of rapamycin, and peroxisome proliferator-activated receptor signaling pathways; these genes are considered as the most promising candidate regulators of milk fat synthesis. The findings of this study enhance the understanding of the genetic basis and molecular mechanisms of milk fat synthesis, which could lead to the development of cow breeds that produce milk with higher nutritional value.

Bayes Factor-Based Regulatory Gene Network Analysis of Genome-Wide Association Study of Economic Traits in a Purebred Swine Population

Early stage prediction of economic trait performance is important and directly linked to profitability of farm pig production. Genome-wide association study (GWAS) has been applied to find causative genomic regions of traits. This study established a regulatory gene network using GWAS for critical economic pig characteristics, centered on easily measurable body fat thickness in live animals. We genotyped 2,681 pigs using Illumina Porcine SNP60, followed by GWAS to calculate Bayes factors for 47,697 single nucleotide polymorphisms (SNPs) of seven traits. Using this information, SNPs were annotated with specific genes near genome locations to establish the association weight matrix. The entire network consisted of 226 nodes and 6,921 significant edges. For in silico validation of their interactions, we conducted regulatory sequence analysis of predicted target genes of transcription factors (TFs). Three key regulatory TFs were identified to guarantee maximum coverage: AT-rich interaction domain 3B (ARID3B), glial cell missing homolog 1 (GCM1), and GLI family zinc finger 2 (GLI2). We identified numerous genes targeted by ARID3B, associated with cellular processes. GCM1 and GLI2 were involved in developmental processes, and their shared target genes regulated multicellular organismal process. This system biology-based function analysis might contribute to enhancing understanding of economic pig traits.

Genetic variants in the SIRT6 transcriptional regulatory region affect gene activity and carcass quality traits in indigenous Chinese beef cattle (Bos taurus)

Background

The aim of this study was to analyze potential influences of polymorphisms within the regulatory region of the bovine SIRT6 gene on carcass quality traits. Expression analyses suggested that SIRT6 gene is predominately expressed in kidney, compared with other tissues. In 535 indigenous Chinese beef cattle, two novel single nucleotide polymorphisms (SNPs) were identified within the promoter region of the SIRT6 gene.

Results

Association analysis indicated that G allele of the c.-1100 A > G had a positive effect on fat deposition, and the Hap4/4 diplotype had more favourable results than other dipoltypes with respect to the evaluation of carcass quality traits. Furthermore, promoter activity associated with the Hap3 haplotype was measured at higher levels than the Hap1 haplotype, which would be in agreement with the previously described association analysis.

Conclusion

The SIRT6 promoter variants significantly affect transcriptional levels and subsequently significantly influence bovine intramscular fat content.

Genomic regions underlying uniformity of yearling weight in Nellore cattle evaluated under different response variables

BACKGROUND

In livestock, residual variance has been studied because of the interest to improve uniformity of production. Several studies have provided evidence that residual variance is partially under genetic control; however, few investigations have elucidated genes that control it. The aim of this study was to identify genomic regions associated with within-family residual variance of yearling weight (YW; N = 423) in Nellore bulls with high density SNP data, using different response variables. For this, solutions from double hierarchical generalized linear models (DHGLM) were used to provide the response variables, as follows: a DGHLM assuming non-null genetic correlation between mean and residual variance (rmv ≠ 0) to obtain deregressed EBV for mean (dEBVm) and residual variance (dEBVv); and a DHGLM assuming rmv = 0 to obtain two alternative response variables for residual variance, dEBVv_r0 and log-transformed variance of estimated residuals (ln_[Formula: see text]).

RESULTS

The dEBVm and dEBVv were highly correlated, resulting in common regions associated with mean and residual variance of YW. However, higher effects on variance than the mean showed that these regions had effects on the variance beyond scale effects. More independent association results between mean and residual variance were obtained when null rmv was assumed. While 13 and 4 single nucleotide polymorphisms (SNPs) showed a strong association (Bayes Factor > 20) with dEBVv and ln_[Formula: see text], respectively, only suggestive signals were found for dEBVv_r0. All overlapping 1-Mb windows among top 20 between dEBVm and dEBVv were previously associated with growth traits. The potential candidate genes for uniformity are involved in metabolism, stress, inflammatory and immune responses, mineralization, neuronal activity and bone formation.

CONCLUSIONS

It is necessary to use a strategy like assuming null rmv to obtain genomic regions associated with uniformity that are not associated with the mean. Genes involved not only in metabolism, but also stress, inflammatory and immune responses, mineralization, neuronal activity and bone formation were the most promising biological candidates for uniformity of YW. Although no clear evidence of using a specific response variable was found, we recommend consider different response variables to study uniformity to increase evidence on candidate regions and biological mechanisms behind it.

Genetic diversity and selection signatures of the beef 'Charolais de Cuba' breed

In this study, we used BovineSNP50 Genotyping BeadChip data to estimate the structure, putative ancestral origin as well as to identify regions with selective sweeps that may have had an important role in the adaptation to tropical conditions of the 'Charolais de Cuba' (CHCU) breed. According to a principal component analysis, CHCU samples cluster together with taurine breeds with an estimated 93% of taurus ancestral alleles. Despite the short period since importation, we detected differentiation (Fst = 0.049) between the French Charolaise (CHA) and CHCU. However, CHA breed was the closest breed to CHCU followed by other hybrids breed with a clear CHA origin. Linkage disequilibrium (r²) decay tends to be lower in CHCU compared to CHA probably due to a less intense artificial selection programs of CHCU. Signals of recent adaptation to tropical conditions between CHCU and CHA were identified. Genes mapping within those regions reflect different functions related to immunity, metabolic changes and heat tolerance (CHCU) and muscle development and meat quality (CHA) that may have had an important role in the phenotypic differentiation of these breeds. Further studies will expand our knowledge on the molecular basis of adaptation of cattle to tropical conditions and molecular process associated with meat quality traits.

Genome-wide scan reveals population stratification and footprints of recent selection in Nelore cattle

Background

This study aimed at (1) assessing the genomic stratification of experimental lines of Nelore cattle that have experienced different selection regimes for growth traits, and (2) identifying genomic regions that have undergone recent selection. We used a sample of 763 animals genotyped with the Illumina BovineHD BeadChip, among which 674 animals originated from two lines that are maintained under directional selection for increased yearling body weight and 89 animals from a control line that is maintained under stabilizing selection.

Results

Multidimensional analysis of the genomic dissimilarity matrix and admixture analysis revealed a substantial level of population stratification between the directional selection lines and the stabilizing selection control line. Two of the three tests used to detect selection signatures (F_ST, XP-EHH and iHS) revealed six candidate regions with indications of selection, which strongly indicates truly positive signals. The set of identified candidate genes included several genes with roles that are functionally related to growth metabolism, such as COL14A1, CPT1C, CRH, TBC1D1, and XKR4.

Conclusions

The current study identified genetic stratification that resulted from almost four decades of divergent selection in an experimental Nelore population, and highlighted autosomal genomic regions that present patterns of recent selection. Our findings provide a basis for a better understanding of the metabolic mechanism that underlies the growth traits, which are modified by selection for yearling body weight.

Electronic supplementary material

The online version of this article (10.1186/s12711-018-0381-2) contains supplementary material, which is available to authorized users.

A Key Gene, PLIN1, Can Affect Porcine Intramuscular Fat Content Based on Transcriptome Analysis

Intramuscular fat (IMF) content is an important indicator for meat quality evaluation. However, the key genes and molecular regulatory mechanisms affecting IMF deposition remain unclear. In the present study, we identified 75 differentially expressed genes (DEGs) between the higher (H) and lower (L) IMF content of pigs using transcriptome analysis, of which 27 were upregulated and 48 were downregulated. Notably, Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis indicated that the DEG perilipin-1 (PLIN1) was significantly enriched in the fat metabolism-related peroxisome proliferator-activated receptor (PPAR) signaling pathway. Furthermore, we determined the expression patterns and functional role of porcine PLIN1. Our results indicate that PLIN1 was highly expressed in porcine adipose tissue, and its expression level was significantly higher in the H IMF content group when compared with the L IMF content group, and expression was increased during adipocyte differentiation. Additionally, our results confirm that PLIN1 knockdown decreases the triglyceride (TG) level and lipid droplet (LD) size in porcine adipocytes. Overall, our data identify novel candidate genes affecting IMF content and provide new insight into PLIN1 in porcine IMF deposition and adipocyte differentiation.

TT genotype of rs2941484 in the human HNF4G gene is associated with hyperuricemia in Chinese Han men

The aim of the study is to investigate the association between the human hepatocyte nuclear factor 4 gamma (HNF4G) gene and hyperuricemia in Chinese Han population. A total of 414 hyperuricemia patients and 406 gender and age-matched normouricemic controls were enrolled. Four single nucleotide polymorphisms were genotyped as genetic markers for the human HNF4G gene (rs2977939, rs1805098, rs2941484, rs4735692). Data were analyzed for two separate groups: men and women. For rs2941484, the genotype distribution frequency in hyperuricemic subjects and was significantly different from that in normouricemic controls in men (P = 0.038). Meanwhile, in recessive model of rs2941484, the distribution frequency of TT genotype and CC+CT genotypes also differed significantly between the hyperuricemia men and normouricemic men (P = 0.011). For the other 3 SNPs in both men and women, there was no difference in the genotype and allele and distribution frequency between the hyperuricemia patients and normouricemic controls. In men, after adjustments for BMI, SBP, DBP, fasting glucose, total cholesterol, triglycerides, low density lipoprotein cholesterol and creatinine, the men with the TT genotype of rs2941484 were found to have significantly higher probability of suffering from hyperuricemia than the ones with CT and CC genotypes (OR = 2.170, P < 0.001). Therefore, TT genotype of rs2941484 in the human HNF4G gene might be a gender-specific genetic marker for hyperuricemia in Chinese Han men.

Multiple-trait genomewide mapping and gene network analysis for scrotal circumference growth curves in Brahman cattle

Fertility traits are economically important in cattle breeding programs. Scrotal circumference (SC) measures are repeatable, easily obtained, highly heritable, and positively correlated with female fertility traits and sperm quality traits in males. A useful approach to summarize SC measures over time is using nonlinear models, which summarize specific measures of SC in a few parameters with biological interpretation. This approach facilitates the selection of bulls with larger SC and maturity index (K), that is, early maturing animals. Because SC is a sex-limited trait, identifying the underlying genomics of growth curve parameters will allow selection across both males and females. We reported the first multitrait genomewide association study (GWAS) of estimated growth curve parameters for SC data in Brahman cattle. Five widely used nonlinear models were tested to fit a total of 3,612 SC records, measured at 6, 12, 18, and 24 mo of age. The von Bertalanffy model, individually fitted for each animal, best fit this SC data. Parameter estimates SC at maturity (A) and K as well as SC at all ages were jointly analyzed in a GWAS to identify 1-Mb regions most strongly associated with each trait. Heritabilities were 0.25 for K and 0.32 for A and ranged from 0.51 to 0.72 for SC at 6 (SC6), 12 (SC12), 18 (SC18), and 24 mo of age (SC24). An overlapping window on chromosome 14 explaining around 0.8% of genetic variance for K, SC12, SC18, and SC24 was observed. The major positional candidate genes within 1 Mb upstream and downstream of this overlapping window were , , , and . Windows of 1 Mb explaining more than 0.4% of each trait on chromosomes 1, 3, 6, 7, 14, 17, 18, 24, 25, and 26 were identified. Pathways and net-work analyses were indicated through transcription factors playing a role on fertility traits: , , , , , , and . Further validation studies on larger populations or other breeds are required to validate these findings and to improve our understanding of the biology and complex genetic architecture of traits associated with scrotal growth and male fertility in cattle.

TRIENNIAL GROWTH AND DEVELOPMENT SYMPOSIUM: Factors influencing bovine intramuscular adipose tissue development and cellularity

Appearance, distribution, and amount of intramuscular fat (IMF), often referred to as marbling, are highly variable and depend on environmental and genetic factors. On the molecular level, the concerted action of several drivers, including hormones, receptors, transcription factors, etc., determines where clusters of adipocytes arise. Therefore, the aim of future studies remains to identify such factors as biological markers of IMF to increase the ability to identify animals that deposit IMF early in age to increase efficiency of high-quality meat production. In an attempt to unravel the cellular development of marbling, we investigated the abundance of markers for adipogenic differentiation during fattening of cattle and the transcriptome of muscle and dissected IMF. Markers of different stages of adipogenic differentiation are well known from cell culture experiments. They are usually transiently expressed, such as delta-like homolog 1 (DLK1) that is abundant in preadipocytes and absent during differentiation to mature adipocytes. It is even a greater challenge to detect those markers in live animals. Within skeletal muscles, hyperplasia and hypertrophy of adipocytes can be observed throughout life. Therefore, development of marbling requires, on the cellular level, recruitment, proliferation, and differentiation of adipogenic cells to store excess energy in the form of lipids in new cells. In a recent study, we investigated the localization and abundance of early markers of adipogenic differentiation, such as DLK1, in bovine muscle tissue. An inverse relationship between IMF content and number of DLK1-positive cells in bovine muscle was demonstrated. Considering the cellular environment of differentiating adipocytes in muscle and the secretory action of adipocytes and myocytes, it becomes obvious that cross talk between cells via adipokines and myokines may be important for IMF development. Secreted proteins can act on other cells, inhibiting or stimulating their function via autocrine and paracrine actions. Such factors with potential influence on IMF, among them, agouti signaling protein and thrombospondin 4, were identified in transcriptome analyses and further investigated. Furthermore, results from transcriptome analysis indicate involvement of genes that are not directly related to adipogenesis and lipid metabolism, providing new candidates for future research.

Candidate genes for male and female reproductive traits in Canchim beef cattle

Background

Beef cattle breeding programs in Brazil have placed greater emphasis on the genomic study of reproductive traits of males and females due to their economic importance. In this study, genome-wide associations were assessed for scrotal circumference at 210 d of age, scrotal circumference at 420 d of age, age at first calving, and age at second calving, in Canchim beef cattle. Data quality control was conducted resulting in 672,778 SNPs and 392 animals.

Results

Associated SNPs were observed for scrotal circumference at 420 d of age (435 SNPs), followed by scrotal circumference at 210 d of age (12 SNPs), age at first calving (six SNPs), and age at second calving (four SNPs). We investigated whether significant SNPs were within genic or surrounding regions. Biological processes of genes were associated with immune system, multicellular organismal process, response to stimulus, apoptotic process, cellular component organization or biogenesis, biological adhesion, and reproduction.

Conclusions

Few associations were observed for scrotal circumference at 210 d of age, age at first calving, and age at second calving, reinforcing their polygenic inheritance and the complexity of understanding the genetic architecture of reproductive traits. Finding many associations for scrotal circumference at 420 d of age in various regions of the Canchim genome also reveals the difficulty of targeting specific candidate genes that could act on fertility; nonetheless, the high linkage disequilibrium between loci herein estimated could aid to overcome this issue. Therefore, all relevant information about genomic regions influencing reproductive traits may contribute to target candidate genes for further investigation of causal mutations and aid in future genomic studies in Canchim cattle to improve the breeding program.

Electronic supplementary material

The online version of this article (doi:10.1186/s40104-017-0199-8) contains supplementary material, which is available to authorized users.

The effect of Bu Zhong Yi Qi decoction on simulated weightlessness‑induced muscle atrophy and its mechanisms

Microgravity has been previously demonstrated to induce skeletal muscle atrophy, loss of muscle force and disorders in myogenesis and metabolism. Current pharmacological strategies exhibit poor efficacy. Bu Zhong Yi Qi decoction (BZ) is a well-known traditional Chinese medicine decoction used for myasthenia gravis. In the present study, its effect on unloading induced muscle atrophy was investigated. The mousetail suspension model was used to simulate weightlessness induced muscle atrophy. The results indicated that BZ could significantly protect muscles from simulated weightlessness-induced atrophy. To elucidate the underlying mechanisms, drugCIPHER-CS methods were introduced to predict its potential targets, significantly enriched pathways and biological processes. The results demonstrated that the calcium signaling pathway, citrate cycle, biosynthetic and lipid metabolic process are affected by BZ. Among the targets, nuclear receptor corepressor 1 (NCoR1) is one of the most important proteins involved in myogenesis and metabolism. The results indicated that BZ significantly downregulated NCoR 1 expression, and further induced muscle differentiation and metabolism by regulating NCoR1-associated gene expression in vivo and in vitro. In summary, the present study indicated that may be effective in combating weightlessness-induced muscle atrophy. Combined with bioinformatics, the underlying mechanism for this decoction was investigated, which provided an improved understanding of this decoction.

Locus-specific ancestry to detect recent response to selection in admixed Swiss Fleckvieh cattle

Identification of selection signatures is one of the current endeavors of evolutionary genetics. Admixed populations may be used to infer post-admixture selection. We calculated local ancestry for Swiss Fleckvieh, a composite of Simmental (SI) and Red Holstein Friesian (RHF), to infer such signals. Illumina Bovine SNP50 BeadChip data for 300 admixed, 88 SI and 97 RHF bulls were used. The average RHF ancestry across the whole genome was 0.70. To identify regions with high deviation from average, we considered two significance thresholds, based on a permutation test and extreme deviation from normal distribution. Regions on chromosomes 13 (46.3-47.3 Mb) and 18 (18.7-25.9 Mb) passed both thresholds in the direction of increased SI. Extended haplotype homozygosity within (iHS) and between (Rsb) populations was calculated to explore additional patterns of pre- and post-admixture selection signals. The Rsb score of admixed and SI was significant in a wide region of chromosome 18 (6.6-24.6 Mb) overlapped with one area of strong local ancestry deviation. FTO, with pleiotropic effect on milk and fertility, NOD2 on dairy and NKD1 and SALL1 on fertility traits are located there. Genetic differentiation of RHF and SI (F_st ), an alternative indicator of pre-admixture selection in pure populations, was calculated. No considerable overlap of peaks of local ancestry deviations and F_st was observed. We found two regions with significant signatures of post-admixture selection in this very young composite, applying comparatively stringent significance thresholds. The signals cover relatively large genomic areas and did not allow pinpointing of the gene(s) responsible for the apparent shift in ancestry proportions.

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.

Multi-breed and multi-trait co-association analysis of meat tenderness and other meat quality traits in three French beef cattle breeds

Background

Studies to identify markers associated with beef tenderness have focused on Warner–Bratzler shear force (WBSF) but the interplay between the genes associated with WBSF has not been explored. We used the association weight matrix (AWM), a systems biology approach, to identify a set of interacting genes that are co-associated with tenderness and other meat quality traits, and shared across the Charolaise, Limousine and Blonde d’Aquitaine beef cattle breeds.

Results

Genome-wide association studies were performed using ~500K single nucleotide polymorphisms (SNPs) and 17 phenotypes measured on more than 1000 animals for each breed. First, this multi-trait approach was applied separately for each breed across 17 phenotypes and second, between- and across-breed comparisons at the AWM and functional levels were performed. Genetic heterogeneity was observed, and most of the variants that were associated with WBSF segregated within rather than across breeds. We identified 206 common candidate genes associated with WBSF across the three breeds. SNPs in these common genes explained between 28 and 30 % of the phenotypic variance for WBSF. A reduced number of common SNPs mapping to the 206 common genes were identified, suggesting that different mutations may target the same genes in a breed-specific manner. Therefore, it is likely that, depending on allele frequencies and linkage disequilibrium patterns, a SNP that is identified for one breed may not be informative for another unrelated breed. Well-known candidate genes affecting beef tenderness were identified. In addition, some of the 206 common genes are located within previously reported quantitative trait loci for WBSF in several cattle breeds. Moreover, the multi-breed co-association analysis detected new candidate genes, regulators and metabolic pathways that are likely involved in the determination of meat tenderness and other meat quality traits in beef cattle.

Conclusions

Our results suggest that systems biology approaches that explore associations of correlated traits increase statistical power to identify candidate genes beyond the one-dimensional approach. Further studies on the 206 common genes, their pathways, regulators and interactions will expand our knowledge on the molecular basis of meat tenderness and could lead to the discovery of functional mutations useful for genomic selection in a multi-breed beef cattle context.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-016-0216-y) contains supplementary material, which is available to authorized users.