Immune status, well-being and gut microbiota in military supplemented with synbiotic ice cream and submitted to field training: a randomised clinical trial.

Strenuous physical activity, sleep deprivation and psychological stress are common features of military field training. The present study aimed to evaluate the effects of supplementation with a synbiotic ice cream on salivary IgA, gastrointestinal symptoms, well-being indicators and gut microbiota in young military participants undergoing field training. Sixty-five military completed the study: one group was supplemented for 30 d with synbiotic ice cream containing: 2·1 × 108 CFU/g for Lactobacillus acidophilus LA-5 and 2·7 × 109 CFU/g for Bifidobacterium animalis BB-12 and 2·3 g of inulin in the 60 g of ice cream at manufacture, and the other with a placebo ice cream. Volunteers were evaluated at pre-supplementation (baseline), post-supplementation and after a 5-d military training. Bifidobacterium and Lactobacillus genera were measured in stool samples and both showed a higher differential abundance post-supplementation and training. Salivary IgA and gastrointestinal symptoms decreased at post-training in both groups (P < 0·05; main effect of time); however, supplementation with synbiotic did not mitigate this effect. Tenseness and sleepiness were decreased in the synbiotic-treated group, but not in the placebo group at post-military training (P = 0·01 and 0·009, respectively; group × time effect). The other well-being indicators were not affected by the synbiotic supplementation. In conclusion, 30 d of synbiotic ice cream supplementation containing inulin, L. acidophilus LA-5 and B. animalis BB-12 favourably modulated gut microbiota and improved tenseness and sleepiness in healthy young military undergoing a 5-d field training. These improvements may be relevant to this population as they may influence the decision-making process in an environment of high physical and psychological stress.

Identification of putative regulatory regions and transcription factors associated with intramuscular fat content traits.

BACKGROUND: Integration of high throughput DNA genotyping and RNA-sequencing data allows for the identification of genomic regions that control gene expression, known as expression quantitative trait loci (eQTL), on a whole genome scale. Intramuscular fat (IMF) content and carcass composition play important roles in metabolic and physiological processes in mammals because they influence insulin sensitivity and consequently prevalence of metabolic diseases such as obesity and type 2 diabetes. However, limited information is available on the genetic variants and mechanisms associated with IMF deposition in mammals. Thus, our hypothesis was that eQTL analyses could identify putative regulatory regions and transcription factors (TFs) associated with intramuscular fat (IMF) content traits. RESULTS: We performed an integrative eQTL study in skeletal muscle to identify putative regulatory regions and factors associated with intramuscular fat content traits. Data obtained from skeletal muscle samples of 192 animals was used for association analysis between 461,466 SNPs and the transcription level of 11,808 genes. This yielded 1268 cis- and 10,334 trans-eQTLs, among which we identified nine hotspot regions that each affected the expression of > 119 genes. These putative regulatory regions overlapped with previously identified QTLs for IMF content. Three of the hotspots respectively harbored the transcription factors USF1, EGR4 and RUNX1T1, which are known to play important roles in lipid metabolism. From co-expression network analysis, we further identified modules significantly correlated with IMF content and associated with relevant processes such as fatty acid metabolism, carbohydrate metabolism and lipid metabolism. CONCLUSION: This study provides novel insights into the link between genotype and IMF content as evident from the expression level. It thereby identifies genomic regions of particular importance and associated regulatory factors. These new findings provide new knowledge about the biological processes associated with genetic variants and mechanisms associated with IMF deposition in mammals.

Differences in the skeletal muscle transcriptome profile associated with extreme values of fatty acids content.

BACKGROUND: Lipids are a class of molecules that play an important role in cellular structure and metabolism in all cell types. In the last few decades, it has been reported that long-chain fatty acids (FAs) are involved in several biological functions from transcriptional regulation to physiological processes. Several fatty acids have been both positively and negatively implicated in different biological processes in skeletal muscle and other tissues. To gain insight into biological processes associated with fatty acid content in skeletal muscle, the aim of the present study was to identify differentially expressed genes (DEGs) and functional pathways related to gene expression regulation associated with FA content in cattle. RESULTS: Skeletal muscle transcriptome analysis of 164 Nellore steers revealed no differentially expressed genes (DEGs, FDR 10%) for samples with extreme values for linoleic acid (LA) or stearic acid (SA), and only a few DEGs for eicosapentaenoic acid (EPA, 5 DEGs), docosahexaenoic acid (DHA, 4 DEGs) and palmitic acid (PA, 123 DEGs), while large numbers of DEGs were associated with oleic acid (OA, 1134 DEGs) and conjugated linoleic acid cis9 trans11 (CLA-c9t11, 872 DEGs). Functional annotation and functional enrichment from OA DEGs identified important genes, canonical pathways and upstream regulators such as SCD, PLIN5, UCP3, CPT1, CPT1B, oxidative phosphorylation mitochondrial dysfunction, PPARGC1A, and FOXO1. Two important genes associated with lipid metabolism, gene expression and cancer were identified as DEGs between animals with high and low CLA-c9t11, specifically, epidermal growth factor receptor (EGFR) and RNPS. CONCLUSION: Only two out of seven classes of molecules of FA studied were associated with large changes in the expression profile of skeletal muscle. OA and CLA-c9t11 content had significant effects on the expression level of genes related to important biological processes associated with oxidative phosphorylation, and cell growth, survival, and migration. These results contribute to our understanding of how some FAs modulate metabolism and may have protective health function.

Genomic structure and marker-derived gene networks for growth and meat quality traits of Brazilian Nelore beef cattle.

BACKGROUND: Nelore is the major beef cattle breed in Brazil with more than 130 million heads. Genome-wide association studies (GWAS) are often used to associate markers and genomic regions to growth and meat quality traits that can be used to assist selection programs. An alternative methodology to traditional GWAS that involves the construction of gene network interactions, derived from results of several GWAS is the AWM (Association Weight Matrices)/PCIT (Partial Correlation and Information Theory). With the aim of evaluating the genetic architecture of Brazilian Nelore cattle, we used high-density SNP genotyping data (~770,000 SNP) from 780 Nelore animals comprising 34 half-sibling families derived from highly disseminated and unrelated sires from across Brazil. The AWM/PCIT methodology was employed to evaluate the genes that participate in a series of eight phenotypes related to growth and meat quality obtained from this Nelore sample. RESULTS: Our results indicate a lack of structuring between the individuals studied since principal component analyses were not able to differentiate families by its sires or by its ancestral lineages. The application of the AWM/PCIT methodology revealed a trio of transcription factors (comprising VDR, LHX9 and ZEB1) which in combination connected 66 genes through 359 edges and whose biological functions were inspected, some revealing to participate in biological growth processes in literature searches. CONCLUSIONS: The diversity of the Nelore sample studied is not high enough to differentiate among families neither by sires nor by using the available ancestral lineage information. The gene networks constructed from the AWM/PCIT methodology were a useful alternative in characterizing genes and gene networks that were allegedly influential in growth and meat quality traits in Nelore cattle.

Lessons for livestock genomics from genome and transcriptome sequencing in cattle and other mammals.

BACKGROUND: Decreasing sequencing costs and development of new protocols for characterizing global methylation, gene expression patterns and regulatory regions have stimulated the generation of large livestock datasets. Here, we discuss experiences in the analysis of whole-genome and transcriptome sequence data. METHODS: We analyzed whole-genome sequence (WGS) data from 132 individuals from five canid species (Canis familiaris, C. latrans, C. dingo, C. aureus and C. lupus) and 61 breeds, three bison (Bison bison), 64 water buffalo (Bubalus bubalis) and 297 bovines from 17 breeds. By individual, data vary in extent of reference genome depth of coverage from 4.9X to 64.0X. We have also analyzed RNA-seq data for 580 samples representing 159 Bos taurus and Rattus norvegicus animals and 98 tissues. By aligning reads to a reference assembly and calling variants, we assessed effects of average depth of coverage on the actual coverage and on the number of called variants. We examined the identity of unmapped reads by assembling them and querying produced contigs against the non-redundant nucleic acids database. By imputing high-density single nucleotide polymorphism data on 4010 US registered Angus animals to WGS using Run4 of the 1000 Bull Genomes Project and assessing the accuracy of imputation, we identified misassembled reference sequence regions. RESULTS: We estimate that a 24X depth of coverage is required to achieve 99.5 % coverage of the reference assembly and identify 95 % of the variants within an individual's genome. Genomes sequenced to low average coverage (e.g., <10X) may fail to cover 10 % of the reference genome and identify <75 % of variants. About 10 % of genomic DNA or transcriptome sequence reads fail to align to the reference assembly. These reads include loci missing from the reference assembly and misassembled genes and interesting symbionts, commensal and pathogenic organisms. CONCLUSIONS: Assembly errors and a lack of annotation of functional elements significantly limit the utility of the current draft livestock reference assemblies. The Functional Annotation of Animal Genomes initiative seeks to annotate functional elements, while a 70X Pac-Bio assembly for cow is underway and may result in a significantly improved reference assembly.

What's in your next-generation sequence data? An exploration of unmapped DNA and RNA sequence reads from the bovine reference individual.

BACKGROUND: Next-generation sequencing projects commonly commence by aligning reads to a reference genome assembly. While improvements in alignment algorithms and computational hardware have greatly enhanced the efficiency and accuracy of alignments, a significant percentage of reads often remain unmapped. RESULTS: We generated de novo assemblies of unmapped reads from the DNA and RNA sequencing of the Bos taurus reference individual and identified the closest matching sequence to each contig by alignment to the NCBI non-redundant nucleotide database using BLAST. As expected, many of these contigs represent vertebrate sequence that is absent, incomplete, or misassembled in the UMD3.1 reference assembly. However, numerous additional contigs represent invertebrate species. Most prominent were several species of Spirurid nematodes and a blood-borne parasite, Babesia bigemina. These species are either not present in the US or are not known to infect taurine cattle and the reference animal appears to have been host to unsequenced sister species. CONCLUSIONS: We demonstrate the importance of exploring unmapped reads to ascertain sequences that are either absent or misassembled in the reference assembly and for detecting sequences indicative of parasitic or commensal organisms.

Global liver gene expression differences in Nelore steers with divergent residual feed intake phenotypes.

BACKGROUND: Efficiency of feed utilization is important for animal production because it can reduce greenhouse gas emissions and improve industry profitability. However, the genetic basis of feed utilization in livestock remains poorly understood. Recent developments in molecular genetics, such as platforms for genome-wide genotyping and sequencing, provide an opportunity to identify genes and pathways that influence production traits. It is known that transcriptional networks influence feed efficiency-related traits such as growth and energy balance. This study sought to identify differentially expressed genes in animals genetically divergent for Residual Feed Intake (RFI), using RNA sequencing methodology (RNA-seq) to obtain information from genome-wide expression profiles in the liver tissues of Nelore cattle. RESULTS: Differential gene expression analysis between high Residual Feed Intake (HRFI, inefficient) and low Residual Feed Intake (LRFI, efficient) groups was performed to provide insights into the molecular mechanisms that underlie feed efficiency-related traits in beef cattle. A total of 112 annotated genes were identified as being differentially expressed between animals with divergent RFI phenotypes. These genes are involved in ion transport and metal ion binding; act as membrane or transmembrane proteins; and belong to gene clusters that are likely related to the transport and catalysis of molecules through the cell membrane and essential mechanisms of nutrient absorption. Genes with functions in cellular signaling, growth and proliferation, cell death and survival were also differentially expressed. Among the over-represented pathways were drug or xenobiotic metabolism, complement and coagulation cascades, NRF2-mediated oxidative stress, melatonin degradation and glutathione metabolism. CONCLUSIONS: Our data provide new insights and perspectives on the genetic basis of feed efficiency in cattle. Some previously identified mechanisms were supported and new pathways controlling feed efficiency in Nelore cattle were discovered. We potentially identified genes and pathways that play key roles in hepatic metabolic adaptations to oxidative stress such as those involved in antioxidant mechanisms. These results improve our understanding of the metabolic mechanisms underlying feed efficiency in beef cattle and will help develop strategies for selection towards the desired phenotype.

Detection of quantitative trait loci for mineral content of Nelore longissimus dorsi muscle.

BACKGROUND: Beef cattle require dietary minerals for optimal health, production and reproduction. Concentrations of minerals in tissues are at least partly genetically determined. Mapping genomic regions that affect the mineral content of bovine longissimus dorsi muscle can contribute to the identification of genes that control mineral balance, transportation, absorption and excretion and that could be associated to metabolic disorders. METHODS: We applied a genome-wide association strategy and genotyped 373 Nelore steers from 34 half-sib families with the Illumina BovineHD BeadChip. Genome-wide association analysis was performed for mineral content of longissimus dorsi muscle using a Bayesian approach implemented in the GenSel software. RESULTS: Muscle mineral content in Bos indicus cattle was moderately heritable, with estimates ranging from 0.29 to 0.36. Our results suggest that variation in mineral content is influenced by numerous small-effect QTL (quantitative trait loci) but a large-effect QTL that explained 6.5% of the additive genetic variance in iron content was detected at 72 Mb on bovine chromosome 12. Most of the candidate genes present in the QTL regions for mineral content were involved in signal transduction, signaling pathways via integral (also called intrinsic) membrane proteins, transcription regulation or metal ion binding. CONCLUSIONS: This study identified QTL and candidate genes that affect the mineral content of skeletal muscle. Our findings provide the first step towards understanding the molecular basis of mineral balance in bovine muscle and can also serve as a basis for the study of mineral balance in other organisms.

Identification of genomic regions associated with feed efficiency in Nelore cattle.

BACKGROUND: Feed efficiency is jointly determined by productivity and feed requirements, both of which are economically relevant traits in beef cattle production systems. The objective of this study was to identify genes/QTLs associated with components of feed efficiency in Nelore cattle using Illumina BovineHD BeadChip (770 k SNP) genotypes from 593 Nelore steers. The traits analyzed included: average daily gain (ADG), dry matter intake (DMI), feed-conversion ratio (FCR), feed efficiency (FE), residual feed intake (RFI), maintenance efficiency (ME), efficiency of gain (EG), partial efficiency of growth (PEG) and relative growth rate (RGR). The Bayes B analysis was completed with Gensel software parameterized to fit fewer markers than animals. Genomic windows containing all the SNP loci in each 1 Mb that accounted for more than 1.0% of genetic variance were considered as QTL region. Candidate genes within windows that explained more than 1% of genetic variance were selected by putative function based on DAVID and Gene Ontology. RESULTS: Thirty-six QTL (1-Mb SNP window) were identified on chromosomes 1, 2, 3, 5, 6, 7, 8, 9, 10, 12, 14, 15, 16, 18, 19, 20, 21, 22, 24, 25 and 26 (UMD 3.1). The amount of genetic variance explained by individual QTL windows for feed efficiency traits ranged from 0.5% to 9.07%. Some of these QTL minimally overlapped with previously reported feed efficiency QTL for Bos taurus. The QTL regions described in this study harbor genes with biological functions related to metabolic processes, lipid and protein metabolism, generation of energy and growth. Among the positional candidate genes selected for feed efficiency are: HRH4, ALDH7A1, APOA2, LIN7C, CXADR, ADAM12 and MAP7. CONCLUSIONS: Some genomic regions and some positional candidate genes reported in this study have not been previously reported for feed efficiency traits in Bos indicus. Comparison with published results indicates that different QTLs and genes may be involved in the control of feed efficiency traits in this Nelore cattle population, as compared to Bos taurus cattle.

Characterization of the rat developmental liver transcriptome.

Gene regulation and transcriptome studies have been enabled by the development of RNA-Seq applications for high-throughput sequencing platforms. Next generation sequencing is remarkably efficient and avoids many issues inherent in hybridization-based microarray methodologies including the exon-specific dependence of probe design. Biologically relevant transcripts including messenger and regulatory RNAs may now be quantified and annotated regardless of whether they have previously been observed. We used RNA-Seq to investigate global patterns of gene expression in early developing rat liver. Liver samples from timed-pregnant Lewis rats were collected at six fetal and neonatal stages [embryonic day (E)14, E16, E18, E20, postnatal day (P)1, P7], transcripts were sequenced using an Illumina HiSeq 2000, and data analysis was performed with the Tuxedo software suite. Genes and isoforms differing in abundance were queried for enrichment within functionally related gene groups using the Functional Annotation Tool of the DAVID Bioinformatics Database. While hematopoietic gene expression is initiated by E14, hepatocyte maturation is a gradual process involving clusters of genes responsible for response to nutrients and enzymes responsible for glycolysis and fatty acid catabolism. Following birth, a large cluster of differentially abundant genes was enriched for mitochondrial gene expression and cholesterol synthesis indicating that by 1 wk of age, the liver is engaged in lipid sensing and bile production. Clustering results for differentially abundant genes and isoforms were similar with the greatest difference for the E14/E16 comparison. Finally, a bioinformatic approach was used to annotate 1,307 novel liver transcripts assembled from sequences that aligned to intergenic regions of the rat genome.

Identification of KCNJ11 as a functional candidate gene for bovine meat tenderness.

The potassium inwardly rectifying channel, subfamily J, member 11 (KCNJ11) gene was investigated as a candidate for meat tenderness based on the effects reported on muscle for KCNJ11 gene knockout in rat models and its position in a quantitative trait locus (QTL) for meat tenderness in the bovine genome. Sequence variations in the KCNJ11 gene were described by sequencing six amplified fragments, covering almost the entire gene. We identified single nucleotide polymorphisms (SNP) and validated them by different approaches, taking advantage of simultaneous projects that are being developed with the same Nelore population. By sequencing the KCNJ11 in Nelore steers representing extreme phenotypes for Warner-Bratzler shear force (WBSF), it was possible to identify 22 SNPs. We validated two of the identified markers by genotyping the whole population (n = 460). Analysis of association between genotypes and WBSF values revealed a significant additive effect of a SNP at different meat aging times (P ≤ 0.05). In addition, an association between the expression levels of KCNJ11 and WBSF was found, with lower expression levels of KCNJ11 associated with more tender meat (P ≤ 0.05). The results showed that the KCNJ11 gene is a candidate mapped to a QTL for meat tenderness previously identified on BTA15 and may be useful to identify animals with genetic potential to produce tender meat. The effect of KCNJ11 observed on muscle is potentially due to changes in activity of KATP channels, which in turn influence the flow of potassium in the intracellular space, allowing establishment of the membrane potential necessary for muscle contraction.

Metagenomics analysis of the morphological aspects and bacterial composition of broiler feces.

In this descriptive study, we used metagenomics to analyze the relationship between the morphological aspects of chicken feces and its respective bacterial compositions. The microbiota composition was determined by sequencing the V4 region of the 16S rRNA genes collected from fresh broiler feces at 19 d old. In total, 48 samples were collected and divided into 8 groups of 6 samples each. The morphological changes studied were feed passage (FP) and reddish mucus (RM). Each was classified into 3 levels of intensity: 1 (slight), 2 (moderate), or 3 (intense). Thus, the 8 groups studied were feed passage (FP-1; FP-2; FP-3), reddish mucus (RM-1; RM-2; RM-3), normal ileal feces (NIF), and cecal discharge (CD). The alpha diversity (Shannon's index) revealed that the CD group showed greater diversity, and was significantly different from FP-2, FP-3, and RM-1. The beta diversity showed that the CD group samples were more homogeneous than the ileal feces groups. The relative abundance analysis revealed that Firmicutes and Proteobacteria were the most abundant phyla in the ileal feces groups. In CD, Firmicutes and Bacteroidetes were the most abundant. The relative abundance at the genus level revealed 136 different bacterial genera. In the ileal feces groups, the two most abundant genera were Lactobacillus and Escherichia/Shigella, except in the FP-1 and RM-2 groups, which had the opposite order. Unlike the others, the CD group had a higher abundance of Bacteroides and Faecalibacterium. When comparing the NIF group with the others, significant changes were found in the fecal microbiota, with nine genera for the FP groups, 19 for the RM groups, and 61 when compared to CD. The results of the present study suggest that evaluation of fecal morphology is a fundamental task that makes it possible to act quickly and assertively, as the morphological aspects of the feces may be related to the composition and structure of fecal microbiota.

The structure of microbial populations in Nelore GIT reveals inter-dependency of methanogens in feces and rumen.

BACKGROUND: The success of different species of ruminants in the colonization of a diverse range of environments is due to their ability to digest and absorb nutrients from cellulose, a complex polysaccharide found in leaves and grass. Ruminants rely on a complex and diverse microbial community, or microbiota, in a unique compartment known as the rumen to break down this polysaccharide. Changes in microbial populations of the rumen can affect the host's development, health, and productivity. However, accessing the rumen is stressful for the animal. Therefore, the development and use of alternative sampling methods are needed if this technique is to be routinely used in cattle breeding. To this end, we tested if the fecal microbiome could be used as a proxy for the rumen microbiome due to its accessibility. We investigated the taxonomic composition, diversity and inter-relations of two different GIT compartments, rumen and feces, of 26 Nelore (Bos indicus) bulls, using Next Generation Sequencing (NGS) metabarcoding of bacteria, archaea and ciliate protozoa. RESULTS: We identified 4265 Amplicon Sequence Variants (ASVs) from bacteria, 571 from archaea, and 107 from protozoa, of which 143 (96 bacteria and 47 archaea) were found common between both microbiomes. The most prominent bacterial phyla identified were Bacteroidetes (41.48%) and Firmicutes (56.86%) in the ruminal and fecal microbiomes, respectively, with Prevotella and Ruminococcaceae UCG-005 the most relatively abundant genera identified in each microbiome. The most abundant archaeal phylum identified was Euryarchaeota, of which Methanobrevibacter gottschalkii, a methanogen, was the prevalent archaeal species identified in both microbiomes. Protozoa were found exclusively identified in the rumen with Bozasella/Triplumaria being the most frequent genus identified. Co-occurrence among ruminal and fecal ASVs reinforces the relationship of microorganisms within a biological niche. Furthermore, the co-occurrence of shared archaeal ASVs between microbiomes indicates a dependency of the predominant fecal methanogen population on the rumen population. CONCLUSIONS: Co-occurring microorganisms were identified within the rumen and fecal microbiomes, which revealed a strong association and inter-dependency between bacterial, archaeal and protozoan populations of the same microbiome. The archaeal ASVs identified as co-occurring between GIT compartments corresponded to the methanogenic genera Methanobrevibacter and Methanosphaera and represented 26.34% of the overall archaeal sequencesdiversity in the rumen and 42.73% in feces. Considering that these archaeal ASVs corresponded to a significant part of the overall diversity of both microbiomes, which is much higher if one includes the interactions of these co-occurring with other rumen archaea ASVs, we suggest that fecal methanogens could be used as a proxy of ruminal methanogens.

Potential Biomarkers for Feed Efficiency-Related Traits in Nelore Cattle Identified by Co-expression Network and Integrative Genomics Analyses.

Feed efficiency helps to reduce environmental impacts from livestock production, improving beef cattle profitability. We identified potential biomarkers (hub genes) for feed efficiency, by applying co-expression analysis in Longissimus thoracis RNA-Seq data from 180 Nelore steers. Six co-expression modules were associated with six feed efficiency-related traits (p-value ≤ 0.05). Within these modules, 391 hub genes were enriched for pathways as protein synthesis, muscle growth, and immune response. Trait-associated transcription factors (TFs) ELF1, ELK3, ETS1, FLI1, and TCF4, were identified with binding sites in at least one hub gene. Gene expression of CCDC80, FBLN5, SERPINF1, and OGN was associated with multiple feed efficiency-related traits (FDR ≤ 0.05) and were previously related to glucose homeostasis, oxidative stress, fat mass, and osteoblastogenesis, respectively. Potential regulatory elements were identified, integrating the hub genes with previous studies from our research group, such as the putative cis-regulatory elements (eQTLs) inferred as affecting the PCDH18 and SPARCL1 hub genes related to immune system and adipogenesis, respectively. Therefore, our analyses contribute to a better understanding of the biological mechanisms underlying feed efficiency in bovine and the hub genes disclosed can be used as biomarkers for feed efficiency-related traits in Nelore cattle.

Detection of Co-expressed Pathway Modules Associated With Mineral Concentration and Meat Quality in Nelore Cattle.

Meat quality is a complex trait that is influenced by genetic and environmental factors, which includes mineral concentration. However, the association between mineral concentration and meat quality, and the specific molecular pathways underlying this association, are not well explored. We therefore analyzed gene expression as measured with RNA-seq in Longissimus thoracis muscle of 194 Nelore steers for association with three meat quality traits (intramuscular fat, meat pH, and tenderness) and the concentration of 13 minerals (Ca, Cr, Co, Cu, Fe, K, Mg, Mn, Na, P, S, Se, and Zn). We identified seven sets of co-expressed genes (modules) associated with at least two traits, which indicates that common pathways influence these traits. From pathway analysis of module hub genes, we further found an over-representation for energy and protein metabolism (AMPK and mTOR signaling pathways) in addition to muscle growth, and protein turnover pathways. Among the identified hub genes FASN, ELOV5, and PDE3B are involved with lipid metabolism and were affected by previously identified eQTLs associated to fat deposition. The reported hub genes and over-represented pathways provide evidence of interplay among gene expression, mineral concentration, and meat quality traits. Future studies investigating the effect of different levels of mineral supplementation in the gene expression and meat quality traits could help us to elucidate the regulatory mechanism by which the genes/pathways are affected.

A comprehensive manually-curated compendium of bovine transcription factors.

Transcription factors (TFs) are pivotal regulatory proteins that control gene expression in a context-dependent and tissue-specific manner. In contrast to human, where comprehensive curated TF collections exist, bovine TFs are only rudimentary recorded and characterized. In this article, we present a manually-curated compendium of 865 sequence-specific DNA-binding bovines TFs, which we analyzed for domain family distribution, evolutionary conservation, and tissue-specific expression. In addition, we provide a list of putative transcription cofactors derived from known interactions with the identified TFs. Since there is a general lack of knowledge concerning the regulation of gene expression in cattle, the curated list of TF should provide a basis for an improved comprehension of regulatory mechanisms that are specific to the species.

An integrative transcriptome analysis indicates regulatory mRNA-miRNA networks for residual feed intake in Nelore cattle.

Residual Feed Intake (RFI) is an economically relevant trait in beef cattle. Among the molecular regulatory mechanisms, microRNAs (miRNAs) are an important dimension in post-transcriptional regulation and have been associated with different biological pathways. Here, we performed differential miRNAs expression and weighted gene co-expression network analyses (WGCNA) to better understand the complex interactions between miRNAs and mRNAs expressed in bovine skeletal muscle and liver. MiRNA and mRNA expression data were obtained from Nelore steers that were genetically divergent for RFI (N = 10 [low RFI or feed efficient]; N = 10 [high RFI or feed inefficient]). Differentially expressed and hub miRNAs such as bta-miR-486, bta-miR-7, bta-miR15a, bta-miR-21, bta-miR 29, bta- miR-30b, bta-miR-106b, bta-miR-199a-3p, bta-miR-204, and bta-miR 296 may have a potential role in variation of RFI. Functional enrichment analysis of differentially expressed (DE) miRNA's target genes and miRNA-mRNA correlated modules revealed that insulin, lipid, immune system, oxidative stress and muscle development signaling pathways might potentially be involved in RFI in this population. Our study identified DE miRNAs, miRNA - mRNA regulatory networks and hub miRNAs related to RFI. These findings suggest a possible role of miRNAs in regulation of RFI, providing new insights into the potential molecular mechanisms that control feed efficiency in Nelore cattle.

Tissue Tropism in Host Transcriptional Response to Members of the Bovine Respiratory Disease Complex.

Bovine respiratory disease (BRD) is the most common infectious disease of beef and dairy cattle and is characterized by a complex infectious etiology that includes a variety of viral and bacterial pathogens. We examined the global changes in mRNA abundance in healthy lung and lung lesions and in the lymphoid tissues bronchial lymph node, retropharyngeal lymph node, nasopharyngeal lymph node and pharyngeal tonsil collected at the peak of clinical disease from beef cattle experimentally challenged with either bovine respiratory syncytial virus, infectious bovine rhinotracheitis, bovine viral diarrhea virus, Mannheimia haemolytica or Mycoplasma bovis. We identified signatures of tissue-specific transcriptional responses indicative of tropism in the coordination of host's immune tissue responses to infection by viral or bacterial infections. Furthermore, our study shows that this tissue tropism in host transcriptional response to BRD pathogens results in the activation of different networks of response genes. The differential crosstalk among genes expressed in lymphoid tissues was predicted to be orchestrated by specific immune genes that act as 'key players' within expression networks. The results of this study serve as a basis for the development of innovative therapeutic strategies and for the selection of cattle with enhanced resistance to BRD.

Genome-Enabled Prediction of Breeding Values for Feedlot Average Daily Weight Gain in Nelore Cattle.

Nelore is the most economically important cattle breed in Brazil, and the use of genetically improved animals has contributed to increased beef production efficiency. The Brazilian beef feedlot industry has grown considerably in the last decade, so the selection of animals with higher growth rates on feedlot has become quite important. Genomic selection (GS) could be used to reduce generation intervals and improve the rate of genetic gains. The aim of this study was to evaluate the prediction of genomic-estimated breeding values (GEBV) for average daily weight gain (ADG) in 718 feedlot-finished Nelore steers. Analyses of three Bayesian model specifications [Bayesian GBLUP (BGBLUP), BayesA, and BayesCπ] were performed with four genotype panels [Illumina BovineHD BeadChip, TagSNPs, and GeneSeek High- and Low-density indicus (HDi and LDi, respectively)]. Estimates of Pearson correlations, regression coefficients, and mean squared errors were used to assess accuracy and bias of predictions. Overall, the BayesCπ model resulted in less biased predictions. Accuracies ranged from 0.18 to 0.27, which are reasonable values given the heritability estimates (from 0.40 to 0.44) and sample size (568 animals in the training population). Furthermore, results from Bos taurus indicus panels were as informative as those from Illumina BovineHD, indicating that they could be used to implement GS at lower costs.