FABP4 is a leading candidate gene associated with residual feed intake in growing Holstein calves

Milk transcriptome biomarker identification to enhance feed efficiency and reduce nutritional costs in dairy ewes

In recent years, rising prices for high-quality protein-based feeds have significantly increased nutrition costs. Consequently, investigating strategies to reduce these expenses and improve feed efficiency (FE) have become increasingly important for the dairy sheep industry. This research investigates the impact of nutritional protein restriction (NPR) during prepuberty and FE on the milk transcriptome of dairy Assaf ewes (sampled during the first lactation). To this end, we first compared transcriptomic differences between NPR and control ewes. Subsequently, we evaluated gene expression differences between ewes with divergent FE, using feed conversion ratio (FCR), residual feed intake (RFI), and consensus classifications of high- and low-FE animals for both indices. Lastly, we assess milk gene expression as a predictor of FE phenotype using random forest. No effect was found for the prepubertal NPR on milk performance or FE. Moreover, at the milk transcriptome level, only one gene, HBB, was differentially expressed between the NPR (n = 14) and the control group (n = 14). Further, the transcriptomic analysis between divergent FE sheep revealed 114 differentially expressed genes (DEGs) for RFI index (high-FERFI = 10 vs low-FERFI = 10), 244 for FCR (high-FEFCR = 10 vs low-FEFCR = 10), and 1 016 DEGs between divergent consensus ewes for both indices (high-FEconsensus = 8 vs low-FEconsensus = 8). These results underscore the critical role of selected FE indices for RNA-Seq analyses, revealing that consensus divergent animals for both indices maximise differences in transcriptomic responses. Genes overexpressed in high-FEconsensus ewes were associated with milk production and mammary gland development, while low-FEconsensus genes were linked to higher metabolic expenditure for tissue organisation and repair. The best prediction accuracy for FE phenotype using random forest was obtained for a set of 44 genes consistently differentially expressed across lactations, with Spearman correlations of 0.37 and 0.22 for FCR and RFI, respectively. These findings provide insights into potential sustainability strategies for dairy sheep, highlighting the utility of transcriptomic markers as FE proxies.

A Novel Missense SNP in the Fatty Acid-Binding Protein 4 (FABP4) Gene is Associated with Growth Traits in Karakul and Awassi Sheep

The fatty acid-binding protein 4 (FABP4) plays a crucial role in the transportation and metabolism of fatty acids. It binds to long-chain fatty acids and facilitates their transport within cells. FABP4 is involved in lipid metabolism, insulin sensitivity, inflammation, and energy homeostasis. This study was conducted to assess the association between the FABP4 gene and growth traits in Karakul and Awassi sheep. A PCR-single strand conformation polymorphism (SSCP) protocol was utilized to assess the polymorphism of FABP4 PCR products with growth traits measured at birth, three, six, nine, and twelve-month intervals. One non-synonymous SNP was identified in the second exon, in which threonine was converted to aspartate in the 61st position in FABP4 (p.61Thr > Asp). This novel SNP showed significant associations with all growth traits measured at all age intervals. The results showed that lambs with the TT genotype exhibited higher growth traits than those with the GT and GG genotypes, respectively. The conducted prediction showed a clearly deleterious effect of p.61Thr > Asp on FABP4 structure, which was accompanied by reduced fatty acid binding efficiency. Owing to the predicted damaging effects caused by p.61Thr > Asp on FABP, lower levels of lipid transport and its consequent increased weight gain and other growth trait indices are expected. Therefore, a putative mechanism through which lambs with these genotypes exhibit higher growth traits is proposed. The FABP4 gene is suggested as a promising marker to improve growth traits in Karakul and Awassi sheep. However, more research is required to validate this mechanism.

Unveiling the Genetic Landscape of Feed Efficiency in Holstein Dairy Cows: Insights into Heritability, Genetic Markers, and Pathways via Meta-Analysis

Improving the feeding efficiency of dairy cows is a key component to improve the utilization of land resources and meet the demand for high-quality protein. Advances in genomic methods and omics techniques have made it possible to breed more efficient dairy cows through genomic selection. The aim of this review is to obtain a comprehensive understanding of the biological background of feed efficiency (FE) complex traits in purebred Holstein dairy cows including heritability estimate, and genetic markers, genes, and pathways participating in FE regulation mechanism. Through a literature search, we systematically reviewed the heritability estimation, molecular genetic markers, genes, biomarkers, and pathways of traits related to feeding efficiency in Holstein dairy cows. A meta-analysis based on a random-effects model was performed to combine reported heritability estimates of FE complex. The heritability of residual feed intake, dry matter intake, and energy balance was 0.20, 0.34, and 0.22, respectively, which proved that it was reasonable to include the related traits in the selection breeding program. For molecular genetic markers, a total of 13 single-nucleotide polymorphisms and copy number variance loci, associated genes, and functions were reported to be significant across populations. A total of 169 reported candidate genes were summarized on a large scale, using a higher threshold (adjusted P value < 0.05). Then, the subsequent pathway enrichment of these genes was performed. The important genes reported in the articles were included in a gene list and the gene list was enriched by gene ontology (GO):biological process (BP), and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways analysis. Three GO:BP terms and four KEGG terms were statistically significant, which mainly focused on adenosine triphosphate (ATP) synthesis, electron transport chain, and OXPHOS pathway. Among these pathways, involved genes such as ATP5MC2, NDUFA, COX7A2, UQCR, and MMP are particularly important as they were previously reported. Twenty-nine reported biological mechanisms along with involved genes were explained mainly by four biological pathways (insulin-like growth factor axis, lipid metabolism, oxidative phosphorylation pathways, tryptophan metabolism). The information from this study will be useful for future studies of genomic selection breeding and genetic structures influencing animal FE. A better understanding of the underlying biological mechanisms would be beneficial, particularly as it might address genetic antagonism.

Re-Evaluation of Genotyping Methodologies in Cattle: The Proficiency of Imputation

In dairy cattle, identifying polymorphisms that contribute to complex economical traits such as residual feed intake (RFI) is challenging and demands accurate genotyping. In this study, we compared imputed genotypes (n = 192 cows) to those obtained using the TaqMan and high-resolution melting (HRM) methods (n = 114 cows), for mutations in the FABP4 gene that had been suggested to have a large effect on RFI. Combining the whole genome sequence (n = 19 bulls) and the cows’ BovineHD BeadChip allowed imputing genotypes for these mutations that were verified by Sanger sequencing, whereas, an error rate of 11.6% and 10.7% were encountered for HRM and TaqMan, respectively. We show that this error rate seriously affected the linkage-disequilibrium analysis that supported this gene candidacy over other BTA14 gene candidates. Thus, imputation produced superior genotypes and should also be regarded as a method of choice to validate the reliability of the genotypes obtained by other methodologies that are prone to genotyping errors due to technical conditions. These results support the view that RFI is a complex trait and that searching for the causative sequence variation underlying cattle RFI should await the development of statistical methods suitable to handle additive and epistatic interactions.

Genome Wide Association Study of Beef Traits in Local Alpine Breed Reveals the Diversity of the Pathways Involved and the Role of Time Stratification

Knowledge of the genetic architecture of key growth and beef traits in livestock species has greatly improved worldwide thanks to genome-wide association studies (GWAS), which allow to link target phenotypes to Single Nucleotide Polymorphisms (SNPs) across the genome. Local dual-purpose breeds have rarely been the focus of such studies; recently, however, their value as a possible alternative to intensively farmed breeds has become clear, especially for their greater adaptability to environmental change and potential for survival in less productive areas. We performed single-step GWAS and post-GWAS analysis for body weight (BW), average daily gain (ADG), carcass fleshiness (CF) and dressing percentage (DP) in 1,690 individuals of local alpine cattle breed, Rendena. This breed is typical of alpine pastures, with a marked dual-purpose attitude and good genetic diversity. Moreover, we considered two of the target phenotypes (BW and ADG) at different times in the individuals' life, a potentially important aspect in the study of the traits' genetic architecture. We identified 8 significant and 47 suggestively associated SNPs, located in 14 autosomal chromosomes (BTA). Among the strongest signals, 3 significant and 16 suggestive SNPs were associated with ADG and were located on BTA10 (50-60 Mb), while the hotspot associated with CF and DP was on BTA18 (55-62 MB). Among the significant SNPs some were mapped within genes, such as SLC12A1, CGNL1, PRTG (ADG), LOC513941 (CF), NLRP2 (CF and DP), CDC155 (DP). Pathway analysis showed great diversity in the biological pathways linked to the different traits; several were associated with neurogenesis and synaptic transmission, but actin-related and transmembrane transport pathways were also represented. Time-stratification highlighted how the genetic architectures of the same traits were markedly different between different ages. The results from our GWAS of beef traits in Rendena led to the detection of a variety of genes both well-known and novel. We argue that our results show that expanding genomic research to local breeds can reveal hitherto undetected genetic architectures in livestock worldwide. This could greatly help efforts to map genomic complexity of the traits of interest and to make appropriate breeding decisions.

FABP4 gene has a very large effect on feed efficiency in lactating Israeli Holstein cows

Improving feed efficiency (FE) is a major goal for the livestock industry. Previously, we have identified 48 SNP markers distributed over 32 genes significantly associated with residual feed intake (RFI) in Israeli Holstein male calves, the most significant of which are located in the bovine FABP4 gene. In the present study, we tested associations of eight of the FABP4 markers with RFI and feed conversion ratio (FCR), along with milk composition and feeding behavioral traits, in 114 lactating Israeli Holstein cows. Large allele effects were found, along with large contributions of FABP4 markers to the phenotypic variation [mean contribution of all significant markers (P < 0.05), 15.4 and 12.0% for RFI and FCR, respectively] and genotypic variation [means of all significant markers (P < 0.05), 75.7 and 32.4% in RFI and FCR, respectively]. However, the association of all significant FABP4 markers with FE and milk content traits was found in opposite directions, such that improved FE was accompanied by decreased milk content. Hence, before inclusion in breeding programs, the gain in FE must be economically balanced with the loss in milk contents. On the other hand, these findings imply that in any current improvement program concentrated on milk traits alone, without taking into account the effect on FE, the progress in milk composition is probably accompanied by deterioration of FE. These results, if confirmed in other populations and breeds, set FABP4 as a prime candidate in any marker-assisted selection program targeting FE as a whole and RFI in particular.

Identification of the association between FABP4 gene polymorphisms and milk production traits in Sfakia sheep

The aim of this study was to estimate the effect of variation in the fatty acid binding protein 4 gene (FABP4) on milk production traits in Greek Sfakia sheep. Polymerase chain reaction – single-stranded conformational polymorphism (PCR-SSCP) analysis was used to genotype a total of 374 Sfakia ewes for two regions of FABP4 located around exon 2–intron 2 (Region 1) and exon 3–intron 3 (Region 2). Each month, for a period of 6 months, milk samples were collected from the ewes to measure total milk yield, fat content, protein content, lactose content, non-fat solid content, pH, and somatic cell count (SCC). A general linear model was used to test the association between the variation observed in FABP4 and milk production traits. Four gene variants (A1–A4) were found in Region 1 and two variants (C1–C2) were found in Region 2. In the first region, the FABP4 genotype significantly affected (P<0.05) non-fat solid levels, fat content, and SCC. The presence of the A2 variant was significantly associated (P<0.05) with decreased SCC, while the presence of A4 was significantly associated with decreased milk yield (P<0.01), increased non-fat solid content (P<0.05), decreased fat content (P<0.01), increased lactose content (P<0.05), and increased pH (P<0.05). In the second region, FABP4 genotype had an effect (P<0.05) on protein content and the presence of the C2 variant was associated (P<0.05) with increased protein content, decreased SCC, and lower pH. The results suggest an association between variation in ovine FABP4 and milk production traits in Greek Sfakia sheep. Nevertheless, further analyses in independent sheep populations of increased size will strengthen these findings.

The effect of breed and diet type on the global transcriptome of hepatic tissue in beef cattle divergent for feed efficiency

Background

Feed efficiency is an important economic and environmental trait in beef production, which can be measured in terms of residual feed intake (RFI). Cattle selected for low-RFI (feed efficient) have similar production levels but decreased feed intake, while also emitting less methane. RFI is difficult and expensive to measure and is not widely adopted in beef production systems. However, development of DNA-based biomarkers for RFI may facilitate its adoption in genomic-assisted breeding programmes. Cattle have been shown to re-rank in terms of RFI across diets and age, while also RFI varies by breed. Therefore, we used RNA-Seq technology to investigate the hepatic transcriptome of RFI-divergent Charolais (CH) and Holstein-Friesian (HF) steers across three dietary phases to identify genes and biological pathways associated with RFI regardless of diet or breed.

Results

Residual feed intake was measured during a high-concentrate phase, a zero-grazed grass phase and a final high-concentrate phase. In total, 322 and 33 differentially expressed genes (DEGs) were identified across all diets for CH and HF steers, respectively. Three genes, GADD45G, HP and MID1IP1, were differentially expressed in CH when both the high-concentrate zero-grazed grass diet were offered. Two canonical pathways were enriched across all diets for CH steers. These canonical pathways were related to immune function.

Conclusions

The absence of common differentially expressed genes across all dietary phases and breeds in this study supports previous reports of the re-ranking of animals in terms of RFI when offered differing diets over their lifetime. However, we have identified biological processes such as the immune response and lipid metabolism as potentially associated with RFI divergence emphasising the previously reported roles of these biological processes with respect to RFI.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5906-8) contains supplementary material, which is available to authorized users.

Consistency of feed efficiency ranking and mechanisms associated with inter-animal variation among growing calves

This study investigated the possible mechanisms for explaining interanimal variation in efficiency of feed utilization in intact male Holstein calves. Additionally, we examined whether the feed efficiency (FE) ranking of calves (n = 26) changed due to age and/or diet quality. Calves were evaluated during three periods (P1, P2, and P3) while fed a high-quality diet (calculated mobilizable energy [ME] of 11.8 MJ/kg DM) during P1 and P3, and a low-quality diet (calculated ME of 7.7 MJ/kg DM) during P2. The study periods were 84, 119, and 127 d, respectively. Initial ages of the calves in P1, P2, and P3 were 7, 11, and 15 mo, respectively, and initial body weight (BW) were 245, 367, and 458 kg, respectively. Individual dry matter intake (DMI), average daily gain (ADG), diet digestibility, and heat production (HP) were measured in all periods. The measured FE indexes were: residual feed intake (RFI), the gain-to-feed ratio (G:F), residual gain (RG), residual gain and intake (RIG), the ratio of HP-to-ME intake (HP/MEI), and residual heat production (RHP). For statistical analysis, animals' performance data in each period, were ranked by RFI, and categorized into high-, medium-, and low-RFI groups (H-RFI, M-RFI, and L-RFI). RFI was not correlated with in vivo digestibility, age, BW, BCS, or ADG in all three periods. The L-RFI group had lowest DMI, MEI, HP, retained energy (RE), and RE/ADG. Chemical analysis of the longissimus dorsi muscle shows that the L-RFI group had a higher percentage of protein and a lower percentage of fat compared to the H-RFI group. We suggested that the main mechanism separating L- from H-RFI calves is the protein-to-fat ratio in the deposited tissues. When efficiency was related to kg/day (DMI and ADG) and not to daily retained energy, the selected efficient L-RFI calves deposited more protein and less fat per daily gain than less efficient H-RFI calves. However, when the significant greater heat increment and maintenance energy requirement of protein compared to fat deposition in tissue were considered, we could not exclude the hypothesis that variation in efficiency is partly explained by efficient energy utilization. The ranking classification of calves to groups according to their RFI efficiency was independent of diet quality and age.