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.

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.

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.

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.

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.

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.

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.

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.

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.

Gene expression differences in Longissimus muscle of Nelore steers genetically divergent for residual feed intake.

Residual feed intake (RFI), a measure of feed efficiency (FE), is defined as the difference between the observed and the predictable feed intake considering size and growth of the animal. It is extremely important to beef production systems due to its impact on the allocation of land areas to alternative agricultural production, animal methane emissions, food demand and cost of production. Global differential gene expression analysis between high and low RFI groups (HRFI and LRFI: less and more efficient, respectively) revealed 73 differentially expressed (DE) annotated genes in Longissimus thoracis (LT) muscle of Nelore steers. These genes are involved in the overrepresented pathways Metabolism of Xenobiotics by Cytochrome P450 and Butanoate and Tryptophan Metabolism. Among the DE transcripts were several proteins related to mitochondrial function and the metabolism of lipids. Our findings indicate that observed gene expression differences are primarily related to metabolic processes underlying oxidative stress. Genes involved in the metabolism of xenobiotics and antioxidant mechanisms were primarily down-regulated, while genes responsible for lipid oxidation and ketogenesis were up-regulated in HRFI group. By using LT muscle, this study reinforces our previous findings using liver tissue and reveals new genes and likely tissue-specific regulators playing key-roles in these processes.

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.

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.

Expression of the Bovine NK-Lysin Gene Family and Activity against Respiratory Pathogens.

Unlike the genomes of many mammals that have a single NK-lysin gene, the cattle genome contains a family of four genes, one of which is expressed preferentially in the lung. In this study, we compared the expression of the four bovine NK-lysin genes in healthy animals to animals challenged with pathogens known to be associated with bovine respiratory disease (BRD) using transcriptome sequencing (RNA-seq). The expression of several NK-lysins, especially NK2C, was elevated in challenged relative to control animals. The effects of synthetic peptides corresponding to functional region helices 2 and 3 of each gene product were tested on both model membranes and bio-membranes. Circular dichroism spectroscopy indicated that these peptides adopted a more helical secondary structure upon binding to an anionic model membrane and liposome leakage assays suggested that these peptides disrupt membranes. Bacterial killing assays further confirmed the antimicrobial effects of these peptides on BRD-associated bacteria, including both Pasteurella multocida and Mannhemia haemolytica and an ultrastructural examination of NK-lysin-treated P. multocida cells by transmission electron microscopy revealed the lysis of target membranes. These studies demonstrate that the expanded bovine NK-lysin gene family is potentially important in host defense against pathogens involved in bovine respiratory disease.

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.

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.

Immunological Response to Single Pathogen Challenge with Agents of the Bovine Respiratory Disease Complex: An RNA-Sequence Analysis of the Bronchial Lymph Node Transcriptome.

Susceptibility to bovine respiratory disease (BRD) is multi-factorial and is influenced by stress in conjunction with infection by both bacterial and viral pathogens. While vaccination is broadly used in an effort to prevent BRD, it is far from being fully protective and cases diagnosed from a combination of observed clinical signs without any attempt at identifying the causal pathogens are usually treated with antibiotics. Dairy and beef cattle losses from BRD are profound worldwide and genetic studies have now been initiated to elucidate host loci which underlie susceptibility with the objective of enabling molecular breeding to reduce disease prevalence. In this study, we employed RNA sequencing to examine the bronchial lymph node transcriptomes of controls and beef cattle which had individually been experimentally challenged with bovine respiratory syncytial virus, infectious bovine rhinotracheitis, bovine viral diarrhea virus, Pasteurella multocida, Mannheimia haemolytica or Mycoplasma bovis to identify the genes that are involved in the bovine immune response to infection. We found that 142 differentially expressed genes were located in previously described quantitative trait locus regions associated with risk of BRD. Mutations affecting the expression or amino acid composition of these genes may affect disease susceptibility and could be incorporated into molecular breeding programs. Genes involved in innate immunity were generally found to be differentially expressed between the control and pathogen-challenged animals suggesting that variation in these genes may lead to a heritability of susceptibility that is pathogen independent. However, we also found pathogen-specific expression profiles which suggest that host genetic variation for BRD susceptibility is pathogen dependent.

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.