Transcriptional response to an alternative diet on liver, muscle, and rumen of beef cattle.

Feed cost represents a major economic determinant within cattle production, amounting to an estimated 75% of the total variable costs. Consequently, comprehensive approaches such as optimizing feed utilization through alternative feed sources, alongside the selection of feed-efficient animals, are of great significance. Here, we investigate the effect of two diets, traditional corn-grain fed and alternative by-product based, on 14 phenotypes related to feed, methane emission and production efficiency and on multi-tissue transcriptomics data from liver, muscle, and rumen wall, derived from 52 Nellore bulls, 26 on each diet. To this end, diets were contrasted at the level of phenotype, gene expression, and gene-phenotype network connectivity. As regards the phenotypic level, at a P value < 0.05, significant differences were found in favour of the alternative diet for average daily weight gain at finishing, dry matter intake at finishing, methane emission, carcass yield and subcutaneous fat thickness at the rib-eye muscle area. In terms of the transcriptional level of the 14,776 genes expressed across the examined tissues, we found 487, 484, and 499 genes differentially expressed due to diet in liver, muscle, and rumen, respectively (P value < 0.01). To explore differentially connected phenotypes across both diet-based networks, we focused on the phenotypes with the largest change in average number of connections within diets and tissues, namely methane emission and carcass yield, highlighting, in particular, gene expression changes involving SREBF2, and revealing the largest differential connectivity in rumen and muscle, respectively. Similarly, from examination of differentially connected genes across diets, the top-ranked most differentially connected regulators within each tissue were MEOX1, PTTG1, and BASP1 in liver, muscle, and rumen, respectively. Changes in gene co-expression patterns suggest activation or suppression of specific biological processes and pathways in response to dietary interventions, consequently impacting the phenotype. The identification of genes that respond differently to diets and their associated phenotypic effects serves as a crucial stepping stone for further investigations, aiming to build upon our discoveries. Ultimately, such advancements hold the promise of improving animal welfare, productivity, and sustainability in livestock farming.

Multi-omics data elucidate parasite-host-microbiota interactions and resistance to Haemonchus contortus in sheep.

BACKGROUND: The integration of molecular data from hosts, parasites, and microbiota can enhance our understanding of the complex biological interactions underlying the resistance of hosts to parasites. Haemonchus contortus, the predominant sheep gastrointestinal parasite species in the tropics, causes significant production and economic losses, which are further compounded by the diminishing efficiency of chemical control owing to anthelmintic resistance. Knowledge of how the host responds to infection and how the parasite, in combination with microbiota, modulates host immunity can guide selection decisions to breed animals with improved parasite resistance. This understanding will help refine management practices and advance the development of new therapeutics for long-term helminth control. METHODS: Eggs per gram (EPG) of feces were obtained from Morada Nova sheep subjected to two artificial infections with H. contortus and used as a proxy to select animals with high resistance or susceptibility for transcriptome sequencing (RNA-seq) of the abomasum and 50 K single-nucleotide genotyping. Additionally, RNA-seq data for H. contortus were generated, and amplicon sequence variants (ASV) were obtained using polymerase chain reaction amplification and sequencing of bacterial and archaeal 16S ribosomal RNA genes from sheep feces and rumen content. RESULTS: The heritability estimate for EPG was 0.12. GAST, GNLY, IL13, MGRN1, FGF14, and RORC genes and transcripts were differentially expressed between resistant and susceptible animals. A genome-wide association study identified regions on chromosomes 2 and 11 that harbor candidate genes for resistance, immune response, body weight, and adaptation. Trans-expression quantitative trait loci were found between significant variants and differentially expressed transcripts. Functional co-expression modules based on sheep genes and ASVs correlated with resistance to H. contortus, showing enrichment in pathways of response to bacteria, immune and inflammatory responses, and hub features of the Christensenellaceae, Bacteroides, and Methanobrevibacter genera; Prevotellaceae family; and Verrucomicrobiota phylum. In H. contortus, some mitochondrial, collagen-, and cuticle-related genes were expressed only in parasites isolated from susceptible sheep. CONCLUSIONS: The present study identified chromosome regions, genes, transcripts, and pathways involved in the elaborate interactions between the sheep host, its gastrointestinal microbiota, and the H. contortus parasite. These findings will assist in the development of animal selection strategies for parasite resistance and interdisciplinary approaches to control H. contortus infection in sheep.

Ruminal and feces metabolites associated with feed efficiency, water intake and methane emission in Nelore bulls.

The objectives of this study were twofold: (1) to identify potential differences in the ruminal and fecal metabolite profiles of Nelore bulls under different nutritional interventions; and (2) to identify metabolites associated with cattle sustainability related-traits. We used different nutritional interventions in the feedlot: conventional (Conv; n = 26), and by-product (ByPr, n = 26). Thirty-eight ruminal fluid and 27 fecal metabolites were significantly different (P < 0.05) between the ByPr and Conv groups. Individual dry matter intake (DMI), residual feed intake (RFI), observed water intake (OWI), predicted water intake (WI), and residual water intake (RWI) phenotypes were lower (P < 0.05) in the Conv group, while the ByPr group exhibited lower methane emission (ME) (P < 0.05). Ruminal fluid dimethylamine was significantly associated (P < 0.05) with DMI, RFI, FE (feed efficiency), OWI and WI. Aspartate was associated (P < 0.05) with DMI, RFI, FE and WI. Fecal C22:1n9 was significantly associated with OWI and RWI (P < 0.05). Fatty acid C14:0 and hypoxanthine were significantly associated with DMI and RFI (P < 0.05). The results demonstrated that different nutritional interventions alter ruminal and fecal metabolites and provided new insights into the relationship of these metabolites with feed efficiency and water intake traits in Nelore bulls.

Repressive epigenetic mechanisms, such as the H3K27me3 histone modification, were predicted to affect muscle gene expression and its mineral content in Nelore cattle.

Epigenetic repression has been linked to the regulation of different cell states. In this study, we focus on the influence of this repression, mainly by H3K27me3, over gene expression in muscle cells, which may affect mineral content, a phenotype that is relevant to muscle function and beef quality. Based on the inverse relationship between H3K27me3 and gene expression (i.e., epigenetic repression) and on contrasting sample groups, we computationally predicted regulatory genes that affect muscle mineral content. To this end, we applied the TRIAGE predictive method followed by a rank product analysis. This methodology can predict regulatory genes that might be affected by repressive epigenetic regulation related to mineral concentration. Annotation of orthologous genes, between human and bovine, enabled our investigation of gene expression in the Longissimus thoracis muscle of Bos indicus cattle. The animals under study had a contrasting mineral content in their muscle cells. We identified candidate regulatory genes influenced by repressive epigenetic mechanisms, linking histone modification to mineral content in beef samples. The discovered candidate genes take part in multiple biological pathways, i.e., impulse transmission, cell signalling, immunological, and developmental pathways. Some of these genes were previously associated with mineral content or regulatory mechanisms. Our findings indicate that epigenetic repression can partially explain the gene expression profiles observed in muscle samples with contrasting mineral content through the candidate regulators here identified.

Bivariate GWAS reveals pleiotropic regions among feed efficiency and beef quality-related traits in Nelore cattle.

Feed-efficient cattle selection is among the most leading solutions to reduce cost for beef cattle production. However, technical difficulties in measuring feed efficiency traits had limited the application in livestock. Here, we performed a Bivariate Genome-Wide Association Study (Bi-GWAS) and presented candidate biological mechanisms underlying the association between feed efficiency and meat quality traits in a half-sibling design with 353 Nelore steers derived from 34 unrelated sires. A total of 13 Quantitative Trait Loci (QTL) were found explaining part of the phenotypic variations. An important transcription factor of adipogenesis in cattle, the TAL1 (rs133408775) gene located on BTA3 was associated with intramuscular fat and average daily gain (IMF-ADG), and a region located on BTA20, close to CD180 and MAST4 genes, both related to fat accumulation. We observed a low positive genetic correlation between IMF-ADG (r = 0.30 ± 0.0686), indicating that it may respond to selection in the same direction. Our findings contributed to clarifying the pleiotropic modulation of the complex traits, indicating new QTLs for bovine genetic improvement.

Rumen and fecal microbiomes are related to diet and production traits in Bos indicus beef cattle.

Background: Ruminants harbor a complex microbial community within their gastrointestinal tract, which plays major roles in their health and physiology. Brazil is one of the largest producers of beef in the world and more than 90% of the beef cattle herds are composed of pure and crossbred Nelore (Bos indicus). Despite its importance to the Brazilian economy and human feeding, few studies have characterized the Nelore microbiome. Therefore, using shotgun metagenomics, we investigated the impact of diet on the composition and functionality of the Nelore microbiome, and explored the associations between specific microbial taxa and their functionality with feed efficiency and methane emission. Results: The ruminal microbiome exhibited significantly higher microbial diversity, distinctive taxonomic profile and variations in microbial functionality compared to the fecal microbiome, highlighting the distinct contributions of the microbiomes of these environments. Animals subjected to different dietary treatments exhibited significant differences in their microbiomes' archaeal diversity and in the abundance of 89 genera, as well as in the functions associated with the metabolism of components of each diet. Moreover, depending on the diet, feed-efficient animals and low methane emitters displayed higher microbial diversity in their fecal microbiome. Multiple genera were associated with an increase or decrease of the phenotypes. Upon analyzing the functions attributed to these taxa, we observed significant differences on the ruminal taxa associated with feed efficient and inefficient cattle. The ruminal taxa that characterized feed efficient cattle stood out for having significantly more functions related to carbohydrate metabolism, such as monosaccharides, di-/oligosaccharides and amino acids. The taxa associated with methane emission had functions associated with methanogenesis and the production of substrates that may influence methane production, such as hydrogen and formate. Conclusion: Our findings highlight the significant role of diet in shaping Nelore microbiomes and how its composition and functionality may affect production traits such as feed efficiency and methane emission. These insights provide valuable support for the implementation of novel feeding and biotechnological strategies.

Differential Allele-Specific Expression Revealed Functional Variants and Candidate Genes Related to Meat Quality Traits in B. indicus Muscle.

Traditional transcriptomics approaches have been used to identify candidate genes affecting economically important livestock traits. Regulatory variants affecting these traits, however, remain under covered. Genomic regions showing allele-specific expression (ASE) are under the effect of cis-regulatory variants, being useful for improving the accuracy of genomic selection models. Taking advantage of the better of these two methods, we investigated single nucleotide polymorphisms (SNPs) in regions showing differential ASE (DASE SNPs) between contrasting groups for beef quality traits. For these analyses, we used RNA sequencing data, imputed genotypes and genomic estimated breeding values of muscle-related traits from 190 Nelore (Bos indicus) steers. We selected 40 contrasting unrelated samples for the analysis (N = 20 animals per contrasting group) and used a beta-binomial model to identify ASE SNPs in only one group (i.e., DASE SNPs). We found 1479 DASE SNPs (FDR ≤ 0.05) associated with 55 beef-quality traits. Most DASE genes were involved with tenderness and muscle homeostasis, presenting a co-expression module enriched for the protein ubiquitination process. The results overlapped with epigenetics and phenotype-associated data, suggesting that DASE SNPs are potentially linked to cis-regulatory variants affecting simultaneously the transcription and phenotype through chromatin state modulation.

Allele-specific expression reveals functional SNPs affecting muscle-related genes in bovine.

Single nucleotide polymorphisms showing allele-specific expression (ASE SNPs) are useful for cis-regulatory variants discovery. Despite this potential, there are expensive costs involved in genome-level ASE analysis for large sample sizes. If different data resolutions are available, genotype imputation can be used to mitigate this limitation. Aiming to increase the power to detect regulatory variants, we used a large dataset (>4 million) of imputed SNP genotypes and RNA-Seq data from 190 Nelore steers. Differences between major and minor allele expressions in muscle were tested with a Binomial Test. We identified 38,177 ASE SNPs (FDR ≤ 0.05) within 7304 linkage disequilibrium blocks. After that, we searched for aseQTLs (i.e., neighboring SNPs potentially regulating the ASE SNPs' allelic expression) by comparing the ASE of heterozygous to homozygous sample groups under a Wilcoxon Rank Sum test. We identified 21,543 aseQTLs potentially regulating 430 ASE SNPs (FDR ≤ 0.05). A total of 3333 cis-eQTLs (being 2098 ASE SNPs and 1075 aseQTLs) were associated with the expression of 758 transcripts (FDR ≤ 0.05), demonstrating the cis-regulatory effect of these ASE SNPs and aseQTLs. Data integration showed reproducibility with previous studies in bovine ASE and genomic imprinting. Furthermore, we identified 36,756 novel ASE regions due to the imputation approach. Comparisons with epigenetics data from Functional Annotation of Animal Genomes (FAANG) suggest a regulatory potential of the ASE-related SNPs. The affected genes were enriched in metabolic pathways essential for muscle homeostasis. These findings reinforce the potential of using ASE for discovering cis-regulatory SNPs that may affect muscle-related traits.

Modeling microRNA-driven post-transcriptional regulation using exon-intron split analysis in pigs.

The contribution of microRNAs (miRNAs) to mRNA post-transcriptional regulation has often been explored by the post hoc selection of downregulated genes and determining whether they harbor binding sites for miRNAs of interest. This approach, however, does not discriminate whether these mRNAs are also downregulated at the transcriptional level. Here, we have characterized the transcriptional and post-transcriptional changes in mRNA expression in two porcine tissues: gluteus medius muscle of fasted and fed Duroc gilts and adipose tissue of lean and obese Duroc-Göttingen minipigs. Exon-intron split analysis of RNA-seq data allowed us to identify downregulated mRNAs with high post-transcriptional signals in fed or obese states, and we assessed whether they harbor binding sites for upregulated miRNAs in any of these two physiological states. We found 26 downregulated mRNAs with high post-transcriptional signals in the muscle of fed gilts and 21 of these were predicted targets of miRNAs upregulated in fed pigs. For adipose tissue, 44 downregulated mRNAs in obese minipigs displayed high post-transcriptional signals, and 25 of these were predicted targets of miRNAs upregulated in the obese state. These results suggest that the contribution of miRNAs to mRNA repression is more prominent in the skeletal muscle system. Finally, we identified several genes that may play relevant roles in the energy homeostasis of the pig skeletal muscle (DKK2 and PDK4) and adipose (SESN3 and ESRRG) tissues. By differentiating transcriptional from post-transcriptional changes in mRNA expression, exon-intron split analysis provides a valuable view of the regulation of gene expression, complementary to canonical differential expression analyses.

Assessing the levels of intraspecific admixture and interspecific hybridization in Iberian wild goats (Capra pyrenaica).

Iberian wild goats (Capra pyrenaica, also known as Iberian ibex, Spanish ibex, and Spanish wild goat) underwent strong genetic bottlenecks during the 19th and 20th centuries due to overhunting and habitat destruction. From the 1970s to 1990s, augmentation translocations were frequently carried out to restock Iberian wild goat populations (very often with hunting purposes), but they were not systematically planned or recorded. On the other hand, recent data suggest the occurrence of hybridization events between Iberian wild goats and domestic goats (Capra hircus). Augmentation translocations and interspecific hybridization might have contributed to increase the diversity of Iberian wild goats. With the aim of investigating this issue, we have genotyped 118 Iberian wild goats from Tortosa-Beceite, Sierra Nevada, Muela de Cortes, Gredos, Batuecas, and Ordesa and Monte Perdido by using the Goat SNP50 BeadChip (Illumina). The analysis of genotypic data indicated that Iberian wild goat populations are strongly differentiated and display low diversity. Only three Iberian wild goats out from 118 show genomic signatures of mixed ancestry, a result consistent with a scenario in which past augmentation translocations have had a limited impact on the diversity of Iberian wild goats. Besides, we have detected eight Iberian wild goats from Tortosa-Beceite with signs of domestic goat introgression. Although rare, hybridization with domestic goats could become a potential threat to the genetic integrity of Iberian wild goats; hence, measures should be taken to avoid the presence of uncontrolled herds of domestic or feral goats in mountainous areas inhabited by this iconic wild ungulate.

Muscle allele-specific expression QTLs may affect meat quality traits in Bos indicus.

Single nucleotide polymorphisms (SNPs) located in transcript sequences showing allele-specific expression (ASE SNPs) were previously identified in the Longissimus thoracis muscle of a Nelore (Bos indicus) population consisting of 190 steers. Given that the allele-specific expression pattern may result from cis-regulatory SNPs, called allele-specific expression quantitative trait loci (aseQTLs), in this study, we searched for aseQTLs in a window of 1 Mb upstream and downstream from each ASE SNP. After this initial analysis, aiming to investigate variants with a potential regulatory role, we further screened our aseQTL data for sequence similarity with transcription factor binding sites and microRNA (miRNA) binding sites. These aseQTLs were overlapped with methylation data from reduced representation bisulfite sequencing (RRBS) obtained from 12 animals of the same population. We identified 1134 aseQTLs associated with 126 different ASE SNPs. For 215 aseQTLs, one allele potentially affected the affinity of a muscle-expressed transcription factor to its binding site. 162 aseQTLs were predicted to affect 149 miRNA binding sites, from which 114 miRNAs were expressed in muscle. Also, 16 aseQTLs were methylated in our population. Integration of aseQTL with GWAS data revealed enrichment for traits such as meat tenderness, ribeye area, and intramuscular fat . To our knowledge, this is the first report of aseQTLs identification in bovine muscle. Our findings indicate that various cis-regulatory and epigenetic mechanisms can affect multiple variants to modulate the allelic expression. Some of the potential regulatory variants described here were associated with the expression pattern of genes related to interesting phenotypes for livestock. Thus, these variants might be useful for the comprehension of the genetic control of these phenotypes.

Multi-Omics Approach Reveals miR-SNPs Affecting Muscle Fatty Acids Profile in Nelore Cattle.

MicroRNAs (miRNAs) are key regulators of gene expression, potentially affecting several biological processes, whose function can be altered by sequence variation. Hence, the integration of single nucleotide polymorphisms (SNP) and miRNAs can explain individual differences in economic traits. To provide new insights into the effects of SNPs on miRNAs and their related target genes, we carried out a multi-omic analysis to identify SNPs in miRNA mature sequences (miR-SNPs) associated with fatty acid (FA) composition in the Nelore cattle. As a result, we identified 3 miR-SNPs in different miRNAs (bta-miR-2419-3p, bta-miR-193a-2, and bta-miR-1291) significantly associated with FA traits (p-value < 0.02, Bonferroni corrected). Among these, the rs110817643C>T, located in the seed sequence of the bta-miR-1291, was associated with different ω6 FAs, polyunsaturated FA, and polyunsaturated:saturated FA ratios. Concerning the other two miR-SNPs, the rs43400521T>C (located in the bta-miR-2419-3p) was associated with C12:0 and C18:1 cis-11 FA, whereas the rs516857374A>G (located in the bta-miR-193a-2) was associated with C18:3 ω6 and ratio of ω6/ω3 traits. Additionally, to identify potential biomarkers for FA composition, we described target genes affected by these miR-SNPs at the mRNA or protein level. Our multi-omics analysis outlines the effects of genetic polymorphism on miRNA, and it highlights miR-SNPs and target candidate genes that control beef fatty acid composition.

Allele-specific expression is widespread in Bos indicus muscle and affects meat quality candidate genes.

Differences between the expression of the two alleles of a gene are known as allele-specific expression (ASE), a common event in the transcriptome of mammals. Despite ASE being a source of phenotypic variation, its occurrence and effects on genetic prediction of economically relevant traits are still unexplored in bovines. Furthermore, as ASE events are likely driven by cis-regulatory mutations, scanning them throughout the bovine genome represents a significant step to elucidate the mechanisms underlying gene expression regulation. To address this question in a Bos indicus population, we built the ASE profile of the skeletal muscle tissue of 190 Nelore steers, using RNA sequencing data and SNPs genotypes from the Illumina BovineHD BeadChip (770 K bp). After quality control, 820 SNPs showed at least one sample with ASE. These SNPs were widespread among all autosomal chromosomes, being 32.01% found in 3'UTR and 31.41% in coding regions. We observed a considerable variation of ASE profile among individuals, which highlighted the need for biological replicates in ASE studies. Functional analysis revealed that ASE genes play critical biological functions in the development and maintenance of muscle tissue. Additionally, some of these genes were previously reported as associated with beef production and quality traits in livestock, thus indicating a possible source of bias on genomic predictions for these traits.

Co-expression network analysis predicts a key role of microRNAs in the adaptation of the porcine skeletal muscle to nutrient supply.

BACKGROUND: The role of non-coding RNAs in the porcine muscle metabolism is poorly understood, with few studies investigating their expression patterns in response to nutrient supply. Therefore, we aimed to investigate the changes in microRNAs (miRNAs), long intergenic non-coding RNAs (lincRNAs) and mRNAs muscle expression before and after food intake. RESULTS: We measured the miRNA, lincRNA and mRNA expression levels in the gluteus medius muscle of 12 gilts in a fasting condition (AL-T0) and 24 gilts fed ad libitum during either 5 h. (AL-T1, N = 12) or 7 h. (AL-T2, N = 12) prior to slaughter. The small RNA fraction was extracted from muscle samples retrieved from the 36 gilts and sequenced, whereas lincRNA and mRNA expression data were already available. In terms of mean and variance, the expression profiles of miRNAs and lincRNAs in the porcine muscle were quite different than those of mRNAs. Food intake induced the differential expression of 149 (AL-T0/AL-T1) and 435 (AL-T0/AL-T2) mRNAs, 6 (AL-T0/AL-T1) and 28 (AL-T0/AL-T2) miRNAs and none lincRNAs, while the number of differentially dispersed genes was much lower. Among the set of differentially expressed miRNAs, we identified ssc-miR-148a-3p, ssc-miR-22-3p and ssc-miR-1, which play key roles in the regulation of glucose and lipid metabolism. Besides, co-expression network analyses revealed several miRNAs that putatively interact with mRNAs playing key metabolic roles and that also showed differential expression before and after feeding. One case example was represented by seven miRNAs (ssc-miR-148a-3p, ssc-miR-151-3p, ssc-miR-30a-3p, ssc-miR-30e-3p, ssc-miR-421-5p, ssc-miR-493-5p and ssc-miR-503) which putatively interact with the PDK4 mRNA, one of the master regulators of glucose utilization and fatty acid oxidation. CONCLUSIONS: As a whole, our results evidence that microRNAs are likely to play an important role in the porcine skeletal muscle metabolic adaptation to nutrient availability.

About the existence of common determinants of gene expression in the porcine liver and skeletal muscle.

BACKGROUND: The comparison of expression QTL (eQTL) maps obtained in different tissues is an essential step to understand how gene expression is genetically regulated in a context-dependent manner. In the current work, we have compared the transcriptomic and eQTL profiles of two porcine tissues (skeletal muscle and liver) which typically show highly divergent expression profiles, in 103 Duroc pigs genotyped with the Porcine SNP60 BeadChip (Illumina) and with available microarray-based measurements of hepatic and muscle mRNA levels. Since structural variation could have effects on gene expression, we have also investigated the co-localization of cis-eQTLs with copy number variant regions (CNVR) segregating in this Duroc population. RESULTS: The analysis of differential expresssion revealed the existence of 1204 and 1490 probes that were overexpressed and underexpressed in the gluteus medius muscle when compared to liver, respectively (|fold-change| > 1.5, q-value < 0.05). By performing genome scans in 103 Duroc pigs with available expression and genotypic data, we identified 76 and 28 genome-wide significant cis-eQTLs regulating gene expression in the gluteus medius muscle and liver, respectively. Twelve of these cis-eQTLs were shared by both tissues (i.e. 42.8% of the cis-eQTLs identified in the liver were replicated in the gluteus medius muscle). These results are consistent with previous studies performed in humans, where 50% of eQTLs were shared across tissues. Moreover, we have identified 41 CNVRs in a set of 350 pigs from the same Duroc population, which had been genotyped with the Porcine SNP60 BeadChip by using the PennCNV and GADA softwares, but only a small proportion of these CNVRs co-localized with the cis-eQTL signals. CONCLUSION: Despite the fact that there are considerable differences in the gene expression patterns of the porcine liver and skeletal muscle, we have identified a substantial proportion of common cis-eQTLs regulating gene expression in both tissues. Several of these cis-eQTLs influence the mRNA levels of genes with important roles in meat (CTSF) and carcass quality (TAPT1), lipid metabolism (TMEM97) and obesity (MARC2), thus evidencing the practical importance of dissecting the genetic mechanisms involved in their expression.

Comparing the mRNA expression profile and the genetic determinism of intramuscular fat traits in the porcine gluteus medius and longissimus dorsi muscles.

BACKGROUND: Intramuscular fat (IMF) content and composition have a strong impact on the nutritional and organoleptic properties of porcine meat. The goal of the current work was to compare the patterns of gene expression and the genetic determinism of IMF traits in the porcine gluteus medius (GM) and longissimus dorsi (LD) muscles. RESULTS: A comparative analysis of the mRNA expression profiles of the pig GM and LD muscles in 16 Duroc pigs with available microarray mRNA expression measurements revealed the existence of 106 differentially expressed probes (fold-change > 1.5 and q-value < 0.05). Amongst the genes displaying the most significant differential expression, several loci belonging to the Hox transcription factor family were either upregulated (HOXA9, HOXA10, HOXB6, HOXB7 and TBX1) or downregulated (ARX) in the GM muscle. Differences in the expression of genes with key roles in carbohydrate and lipid metabolism (e.g. FABP3, ORMDL1 and SLC37A1) were also detected. By performing a GWAS for IMF content and composition traits recorded in the LD and GM muscles of 350 Duroc pigs, we identified the existence of one region on SSC14 (110-114 Mb) displaying significant associations with C18:0, C18:1(n-7), saturated and unsaturated fatty acid contents in both GM and LD muscles. Moreover, we detected several genome-wide significant associations that were not consistently found in both muscles. Further studies should be performed to confirm whether these associations are muscle-specific. Finally, the performance of an eQTL scan for 74 genes, located within GM QTL regions and with available microarray measurements of gene expression, made possible to identify 14 cis-eQTL regulating the expression of 14 loci, and six of them were confirmed by RNA-Seq. CONCLUSIONS: We have detected significant differences in the mRNA expression patterns of the porcine LD and GM muscles, evidencing that the transcriptomic profile of the skeletal muscle tissue is affected by anatomical, metabolic and functional factors. A highly significant association with IMF composition on SSC14 was replicated in both muscles, highlighting the existence of a common genetic determinism, but we also observed the existence of a few associations whose magnitude and significance varied between LD and GM muscles.

Analysing the Expression of Eight Clock Genes in Five Tissues From Fasting and Fed Sows.

In a previous study, we observed that circadian clock genes are differentially expressed in the skeletal muscle of fasting and fed sows. The goal of the current work was to investigate if these genes are also differentially expressed in tissues containing the central (hypothalamus) and peripheral (duodenum, dorsal fat, muscle, and liver) clocks. As animal material, we used 12 sows that fasted 12 h before slaughtering (T0) and 12 sows that were fed ad libitum 7 h prior slaughtering (T2). Tissue samples were collected immediately after slaughter and total RNA was subsequently extracted. The expression of the ARNTL, BHLHE40, CRY2, NPAS2, NR1D1, PER1, PER2, and SIK1 genes was measured by quantitative reverse transcription PCR. The numbers of clock genes showing differential expression before and after feeding varied depending on the tissue i.e., four in dorsal fat and duodenum, six in skeletal muscle, and seven in the liver. In contrast, none of the eight analysed genes displayed a significant differential expression in hypothalamus, the tissue where the central clock resides. This result supports that the differential expression of clock genes in the four tissues mentioned above is probably induced by nutrition and not by the central clock entrained by light. Moreover, we have observed that the NPAS2 and ARNTL genes display positive log2(FC) values in the five tissues under analysis, whilst the CRY2, PER1 (except dorsal fat) and PER2 (except hypothalamus) genes generally show negative log2(FC) values. Such result might be explained by the existence of a negative feedback loop between the ARNTL/NPAS2 and CRY/PER genes. Collectively, these results support that nutrition plays an important role in modulating the timing of porcine peripheral circadian clocks. Such regulation could be essential for coordinating the subsequent metabolic response to nutrient supply.

Genome-wide patterns of homozygosity provide clues about the population history and adaptation of goats.

BACKGROUND: Patterns of homozygosity can be influenced by several factors, such as demography, recombination, and selection. Using the goat SNP50 BeadChip, we genotyped 3171 goats belonging to 117 populations with a worldwide distribution. Our objectives were to characterize the number and length of runs of homozygosity (ROH) and to detect ROH hotspots in order to gain new insights into the consequences of neutral and selection processes on the genome-wide homozygosity patterns of goats. RESULTS: The proportion of the goat genome covered by ROH is, in general, less than 15% with an inverse relationship between ROH length and frequency i.e. short ROH (< 3 Mb) are the most frequent ones. Our data also indicate that ~ 60% of the breeds display low FROH coefficients (< 0.10), while ~ 30 and ~ 10% of the goat populations show moderate (0.10 < FROH < 0.20) or high (> 0.20) FROH values. For populations from Asia, the average number of ROH is smaller and their coverage is lower in goats from the Near East than in goats from Central Asia, which is consistent with the role of the Fertile Crescent as the primary centre of goat domestication. We also observed that local breeds with small population sizes tend to have a larger fraction of the genome covered by ROH compared to breeds with tens or hundreds of thousands of individuals. Five regions on three goat chromosomes i.e. 11, 12 and 18, contain ROH hotspots that overlap with signatures of selection. CONCLUSIONS: Patterns of homozygosity (average number of ROH of 77 and genome coverage of 248 Mb; FROH < 0.15) are similar in goats from different geographic areas. The increased homozygosity in local breeds is the consequence of their small population size and geographic isolation as well as of founder effects and recent inbreeding. The existence of three ROH hotspots that co-localize with signatures of selection demonstrates that selection has also played an important role in increasing the homozygosity of specific regions in the goat genome. Finally, most of the goat breeds analysed in this work display low levels of homozygosity, which is favourable for their genetic management and viability.

Differential expression of mRNA isoforms in the skeletal muscle of pigs with distinct growth and fatness profiles.

BACKGROUND: The identification of genes differentially expressed in the skeletal muscle of pigs displaying distinct growth and fatness profiles might contribute to identify the genetic factors that influence the phenotypic variation of such traits. So far, the majority of porcine transcriptomic studies have investigated differences in gene expression at a global scale rather than at the mRNA isoform level. In the current work, we have investigated the differential expression of mRNA isoforms in the gluteus medius (GM) muscle of 52 Duroc HIGH (increased backfat thickness, intramuscular fat and saturated and monounsaturated fatty acids contents) and LOW pigs (opposite phenotype, with an increased polyunsaturated fatty acids content). RESULTS: Our analysis revealed that 10.9% of genes expressed in the GM muscle generate alternative mRNA isoforms, with an average of 2.9 transcripts per gene. By using two different pipelines, one based on the CLC Genomics Workbench and another one on the STAR, RSEM and DESeq2 softwares, we have identified 10 mRNA isoforms that both pipelines categorize as differentially expressed in HIGH vs LOW pigs (P-value < 0.01 and ±0.6 log2fold-change). Only five mRNA isoforms, produced by the ITGA5, SEMA4D, LITAF, TIMP1 and ANXA2 genes, remain significant after correction for multiple testing (q-value < 0.05 and ±0.6 log2fold-change), being upregulated in HIGH pigs. CONCLUSIONS: The increased levels of specific ITGA5, LITAF, TIMP1 and ANXA2 mRNA isoforms in HIGH pigs is consistent with reports indicating that the overexpression of these four genes is associated with obesity and metabolic disorders in humans. A broader knowledge about the functional attributes of these mRNA variants would be fundamental to elucidate the consequences of transcript diversity on the determinism of porcine phenotypes of economic interest.

Deregulation of obesity-relevant genes is associated with progression in BMI and the amount of adipose tissue in pigs.

The aim of this study was to elucidate the relative impact of three phenotypes often used to characterize obesity on perturbation of molecular pathways involved in obesity. The three obesity-related phenotypes are (1) body mass index (BMI), (2) amount of subcutaneous adipose tissue (SATa), and (3) amount of retroperitoneal adipose tissue (RPATa). Although it is generally accepted that increasing amount of RPATa is 'unhealthy', a direct comparison of the relative impact of the three obesity-related phenotypes on gene expression has, to our knowledge, not been performed previously. We have used multiple linear models to analyze altered gene expression of selected obesity-related genes in tissues collected from 19 female pigs phenotypically characterized with respect to the obesity-related phenotypes. Gene expression was assessed by high-throughput qPCR in RNA from liver, skeletal muscle and abdominal adipose tissue. The stringent statistical approach used in the study has increased the power of the analysis compared to the classical approach of analysis in divergent groups of individuals. Our approach led to the identification of key components of cellular pathways that are modulated in the three tissues in association with changes in the three obesity-relevant phenotypes (BMI, SATa and RPATa). The deregulated pathways are involved in biosynthesis and transcript regulation in adipocytes, in lipid transport, lipolysis and metabolism, and in inflammatory responses. Deregulation seemed more comprehensive in liver (23 genes) compared to abdominal adipose tissue (10 genes) and muscle (3 genes). Notably, the study supports the notion that excess amount of intra-abdominal adipose tissue is associated with a greater metabolic disease risk. Our results provide molecular support for this notion by demonstrating that increasing amount of RPATa has a higher impact on perturbation of cellular pathways influencing obesity and obesity-related metabolic traits compared to increase in BMI and amount of SATa.