Skin transcriptomic analysis reveals candidate genes and pathways associated with thermotolerance in hair sheep.

The skin plays an important role in thermoregulation. Identification of genes on the skin that contribute to increased heat tolerance can be used to select animals with the best performance in warm environments. Our objective was to identify candidate genes associated with the heat stress response in the skin of Santa Ines sheep. A group of 80 sheep assessed for thermotolerance was kept in a climatic chamber for 8 days at a stress level temperature of 36 °C (10 am to 04 pm) and a maintenance temperature of 28 °C (04 pm to 10 am). Two divergent groups, with seven animals each, were formed after ranking them by thermotolerance using rectal temperature. From skin biopsy samples, total RNA was extracted, quantified, and used for RNA-seq analysis. 15,989 genes were expressed in sheep skin samples, of which 4 genes were differentially expressed (DE; FDR < 0.05) and 11 DE (FDR 0.05-0.177) between the two divergent groups. These genes are involved in cellular protection against stress (HSPA1A and HSPA6), ribosome assembly (28S, 18S, and 5S ribosomal RNA), and immune response (IGHG4, GNLY, CXCL1, CAPN14, and SAA-4). The candidate genes and main pathways related to heat tolerance in Santa Ines sheep require further investigation to understand their response to heat stress in different climatic conditions and under solar radiation. It is essential to verify whether these genes and pathways are present in different breeds and to understand the relationship between heat stress and other genes identified in this study.

Brain fatty acid and transcriptome profiles of pig fed diets with different levels of soybean oil.

BACKGROUND: The high similarity in anatomical and neurophysiological processes between pigs and humans make pigs an excellent model for metabolic diseases and neurological disorders. Lipids are essential for brain structure and function, and the polyunsaturated fatty acids (PUFA) have anti-inflammatory and positive effects against cognitive dysfunction in neurodegenerative diseases. Nutrigenomics studies involving pigs and fatty acids (FA) may help us in better understanding important biological processes. In this study, the main goal was to evaluate the effect of different levels of dietary soybean oil on the lipid profile and transcriptome in pigs' brain tissue. RESULTS: Thirty-six male Large White pigs were used in a 98-day study using two experimental diets corn-soybean meal diet containing 1.5% soybean oil (SOY1.5) and corn-soybean meal diet containing 3.0% soybean oil (SOY3.0). No differences were found for the brain total lipid content and FA profile between the different levels of soybean oil. For differential expression analysis, using the DESeq2 statistical package, a total of 34 differentially expressed genes (DEG, FDR-corrected p-value < 0.05) were identified. Of these 34 DEG, 25 are known-genes, of which 11 were up-regulated (log2 fold change ranging from + 0.25 to + 2.93) and 14 were down-regulated (log2 fold change ranging from - 3.43 to -0.36) for the SOY1.5 group compared to SOY3.0. For the functional enrichment analysis performed using MetaCore with the 34 DEG, four pathway maps were identified (p-value < 0.05), related to the ALOX15B (log2 fold change - 1.489), CALB1 (log2 fold change - 3.431) and CAST (log2 fold change + 0.421) genes. A "calcium transport" network (p-value = 2.303e-2), related to the CAST and CALB1 genes, was also identified. CONCLUSION: The results found in this study contribute to understanding the pathways and networks associated with processes involved in intracellular calcium, lipid metabolism, and oxidative processes in the brain tissue. Moreover, these results may help a better comprehension of the modulating effects of soybean oil and its FA composition on processes and diseases affecting the brain tissue.

Transcriptome profile of skeletal muscle using different sources of dietary fatty acids in male pigs.

Pork is of great importance in world trade and represents the largest source of fatty acids in the human diet. Lipid sources such as soybean oil (SOY), canola (CO), and fish oil (FO) are used in pig diets and influence blood parameters and the ratio of deposited fatty acids. In this study, the main objective was to evaluate changes in gene expression in porcine skeletal muscle tissue resulting from the dietary oil sources and to identify metabolic pathways and biological process networks through RNA-Seq. The addition of FO in the diet of pigs led to intramuscular lipid with a higher FA profile composition of C20:5 n-3, C22:6 n-3, and SFA (C16:0 and C18:0). Blood parameters for the FO group showed lower cholesterol and HDL content compared with CO and SOY groups. Skeletal muscle transcriptome analyses revealed 65 differentially expressed genes (DEG, FDR 10%) between CO vs SOY, and 32 DEG for CO vs FO, and 531 DEG for SOY vs FO comparison. Several genes, including AZGP1, PDE3B, APOE, PLIN1, and LIPS, were found to be down-regulated in the diet of the SOY group compared to the FO group. The enrichment analysis revealed DEG involved in lipid metabolism, metabolic diseases, and inflammation between the oil groups, with specific gene functions in each group and altered blood parameters. The results provide mechanisms to help us understand the behavior of genes according to fatty acids.

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.

Heat stress on breeding value prediction for milk yield and composition of a Brazilian Holstein cattle population.

Due to the high milk production of Holstein cows, many countries have chosen to import semen to improve local dairy herds. This strategy would be more effective if this semen was used in the same environment conditions in which the bulls were selected. If the effect of genotype by environment (G × E) interaction is not considered, the estimated breeding values (EBVs) may vary, potentially reducing the selection response. We evaluate the impact of heat stress on selection for milk yield and composition of Holstein cows using random regression models. To verify the interference of heat stress in milk yield (MY) and composition traits (fat, protein, total saturated, and total unsaturated fatty acids content in milk), temperature-humidity index (THI) on test-day milk records was used. The threshold value to divide the environments using test-day information from Brazilian Holstein cows was 72 units of THI, i.e., < 72 represented no heat stress and > 72 represented heat stress. Legendre polynomials of second-order (Leg 2) model and two lactation points (33 and 122 DIM) were used to estimate heritabilities and EBVs for five important dairy traits. The heritabilities of milk components and fatty acids were low (0.09-0.29), regardless of lactation period and degree of heat stress, with the exception of protein content (0.30-0.35). Fat content was the only milk component that was reduced according to the degree of heat stress and lactation period. The EBVs tended to decrease in heat stress conditions, thus animals with high genetic potential demonstrated evidence of G × E interaction. However, acclimatization of dairy cows to heat stress in the farm production systems may have been responsible for the low differences among genetic parameters and EBVs with and without heat stress found in this study.

Influence of successive heat waves on the thermoregulatory responses of pregnant and non-pregnant ewes.

The frequency of heat waves has increased over the last years, with an impact on animal production and health, including the death of animals. Therefore, the aim of this study was to evaluate the dynamics of thermoregulation and hormonal responses in non-pregnant and pregnant ewes exposed to successive heat waves. Twenty-four non-pregnant and 18 pregnant Santa Ines ewes with black coat color (live weight: 55 ± 9.03 kg; age: 60 months) were used. Weather variables such air temperature, relative humidity, and solar radiation were continuously recorded. The rectal and tympanic temperatures and respiratory rate were measured daily. Serum triiodothyronine (T3) and prolactin were evaluated during the heat wave and thermoneutral periods. The physiological variables were higher under the heat wave conditions and were related to the activation of the thermoregulatory system for maintaining homeothermy (P < 0.05). The core body temperature was higher during successive heat waves (P < 0.05), as was the tympanic temperature, which are both affected by changes in air temperature (P < 0.05). T3 and prolactin levels were not influenced by successive heat waves (P < 0.05) and rectal temperature and respiratory rate were highest in non-pregnant ewes (P < 0.05). Prolactin was not affected by temperature. The results indicate that the Santa Ines breed overcomes the thermal challenge during a heat wave without showing severe signs of thermal stress regardless of being pregnant or not.

Genetic parameters for production, quality, and colors from eggs in Brazilian lineages of chickens.

Genetic parameters were estimated for egg production, egg quality, and eggshell colors in eight lineages of Brazilian laying hens. Age at first egg (AFE), total egg production up to the 45th week (PROD), egg weight (EW), albumen height (AH), yolk color (YC), the Haugh units (HU), eggshell strength (ESS), eggshell thickness (EST), yolk weight (YW), eggshell weight (ESW), and eggshell color (L*, a*, and b*) were measured in 2030 eggs obtained from 645 laying hens. Variance components were estimated from a mixed animal model, which included the fixed effects of contemporary groups, cage location, and hen line, and the additive genetic, permanent environmental, and residual as random effects. In general, heritabilities were low to moderate (h2 = 0.11 to 0.48). Genetic correlations among eggshell quality traits were moderate to high (0.36 and 0.69). High genetic correlations were obtained between the eggshell color traits [rg = -0.90 (L* and a*); rg = -0.64 (L* and b*); and rg = 0.65 (a* and b*)]. Results suggest that EW is strongly correlated with ESW, but the genetic correlations between EW and ESS and between EW and EST were low. Genetic correlations between L* and eggshell quality traits were low to moderate, suggesting that L* has little or no relation with external egg quality. However, genetic correlations between a* and b* values and eggshell quality traits were high. The genetic correlations between eggshell color and eggshell quality traits were low, suggesting that the eggshell color does not influence external egg quality. Genetic correlations between PROD and egg quality traits were negative and varied between -0.42 and -0.05. This antagonistic relationship emphasizes the importance of adopting breeding schemes that allow the simultaneous genetic progress of these traits by considering their genetic correlation and economic relevance, such as the selection index.

Prune homolog 2 with BCH domain (PRUNE2) gene expression is associated with feed efficiency-related traits in Nelore steers.

Animal feeding is a critical factor in increasing producer profitability. Improving feed efficiency can help reduce feeding costs and reduce the environmental impact of beef production. Candidate genes previously identified for this trait in differential gene expression studies (e.g., case-control studies) have not examined continuous gene-phenotype variation, which is a limitation. The aim of this study was to investigate the association between the expression of five candidate genes in the liver, measured by quantitative real-time PCR and feed-related traits. We adopted a linear mixed model to associate liver gene expression from 52 Nelore steers with the following production traits: average daily gain (ADG), body weight (BW), dry matter intake (DMI), feed conversion ratio (FCR), feed efficiency (FE), Kleiber index (KI), metabolic body weight (MBW), residual feed intake (RFI), and relative growth ratio (RGR). The total expression of the prune homolog 2 (PRUNE2) gene was significantly associated with DMI, FCR, FE, and RFI (P < 0.05). Furthermore, we have identified a new transcript of PRUNE2 (TCONS_00027692, GenBank MZ041267) that was inversely correlated with FCR and FE (P < 0.05), in contrast to the originally identified PRUNE2 transcript. The cytochrome P450 subfamily 2B (CYP2B6), early growth response protein 1 (EGR1), collagen type I alpha 1 chain (COL1A1), and connective tissue growth factor (CTGF) genes were not associated with any feed efficiency-related traits (P > 0.05). The findings reported herein suggest that PRUNE2 expression levels affects feed efficiency-related traits variation in Nelore steers.

Variants in myostatin and MyoD family genes are associated with meat quality traits in Santa Inês sheep.

Myostatin and MyoD family genes play vital roles in myogenesis and this study aimed to identify association of variants in MyoD1, MyoG, MyF5, MyF6, and MSTN genes with meat traits in Santa Inês sheep. A dataset with 44 variants and records of seven meat traits in 192 lambs (pH0, pH24, a*, b*, L*, tenderness assessed by shear force, and water-holding capacity) was used. Single-locus and haplotype association analyses were performed, and the significance threshold was established according to Bonferroni's method. Single-locus analysis revealed two associations at a Bonferroni level, where the variant c.935-185C > G in MyoD1 had an additive effect (-4.31 ± 1.08 N) on tenderness, while the variant c.464 + 185G > A in MyoG had an additive effect (-2.86 ± 0.64) on a*. Additionally, the haplotype replacement GT>AC in MSTN was associated with pH0 (1.26 ± 0.31), pH24 (1.07 ± 0.27), a* (-1.40 ± 0.51), and tenderness (3.83 ± 1.22 N), while the replacement GT > AG in MyoD1 was associated with pH0 (1.43 ± 0.26), pH24 (1.25 ± 0.22), b* (-1.06 ± 0.39), and tenderness (-4.13 ± 1.16 N). Our results have demonstrated that some variants in MyoG, MyF6, MyoD1, and MSTN can be associated with physicochemical meat traits in Santa Inês sheep.

Genomic integration to identify molecular biomarkers associated with indicator traits of gastrointestinal nematode resistance in sheep.

This study aimed to integrate GWAS and structural variants to propose possible molecular biomarkers related to gastrointestinal nematode resistance traits in Santa Inês sheep. The phenotypic records FAMACHA, haematocrit, white blood cell count, red blood cell count, haemoglobin, platelets and egg counts per gram of faeces were collected from 700 naturally infected animals, belonging to four Brazilian flocks. A total of 576 animals were genotyped using the Ovine SNP12k BeadChip and were imputed using a reference population with Ovine SNP50 BeadChip. The GWAS approaches were based on SNPs, haplotypes, CNVs and ROH. The overlapping between the significant genomic regions detected from all approaches was investigated, and the results were integrated using a network analysis. Genes related to the immune system were found, such as ABCB1, IL6, WNT5A and IRF5. Genomic regions containing candidate genes and metabolic pathways involved in immune responses, inflammatory processes and immune cells affecting parasite resistance traits were identified. The genomic regions, biological processes and candidate genes uncovered could lead to biomarkers for selecting more resilient sheep and improving herd welfare and productivity. The results obtained are the start point to identify molecular biomarkers related to indicator traits of gastrointestinal nematode resistance in Santa Inês sheep.

CNV detection and their association with growth, efficiency and carcass traits in Santa Inês sheep.

Copy number variations (CNV) are an important source of genetic variation. CNV has been increasingly studied and frequently associated with diseases and productive traits in livestock animals. However, CNV-based genome-wide association studies (GWAS) in Santa Inês sheep, one of the principal sheep breeds in Brazil, have not yet been reported. Thus, the aim of this study was to investigate the association between CNV and growth, efficiency and carcass traits in sheep. The Illumina OvineSNP50 BeadChip array was used to detect CNV in 491 Santa Inês individuals. Then, CNV-based GWAS was performed with a linear mixed model approach considering a genomic relationship matrix, for ten traits: (1) growth: body weight at three (W3) and six (W6) months of age; (2) efficiency: residual feed intake (RFI) and feed efficiency (FE) and (3) carcass: external carcass length (ECL), leg length (LL), carcass yield (CY), commercial cuts weight (CCW), loin eye area (LEA) and subcutaneous fat thickness (SFT). We identified 1,167 autosomal CNV in 438 sheep, with 294 non-redundant CNV, ranging from 21.8 to 861.9 kb, merged into 216 distinct copy number variation regions (CNVRs). One significant CNV segment (pFDR -value<0.05) in OAR3 was associated with CY, while another significant CNV in OAR6 was associated with RFI. Additionally, another 5 CNV segments were considered relevant for investigation in the future studies. The significant segments overlapped 4 QTLs and spanned 8 genes, including the SPAST,TGFA and ADGRL3 genes, involved in cell differentiation and energy metabolism. Therefore, the results of the present study increase knowledge about CNV in sheep, their possible impacts on productive traits, and provide information for future investigations, being especially useful for those interested in structural variations in the sheep genome.

Feeding restriction in the pre and postpartum period of hair ewes raised in the semi-arid region: implications on performance and carcass traits of the progeny.

Feed restriction (FR) occurs commonly in sheep production systems in the Brazilian semi-arid region and can cause physiological changes in the progeny. We assessed the effects of the FR pre and postnatal on the performance and carcass traits of Morada Nova lambs. Twenty-four lambs born from a group of 68 ewes were distributed in three treatments of 8 replicates as follows: ewes fed ad libitum pre and postnatal (AL-AL); postnatal restriction (POSTN-R): included lambs born from ewes fed ad libitum in the last third of pregnancy with FR postnatal; and prenatal feeding restriction (PREN-AL): comprised lambs born from ewes subjected to FR in the last third of pregnancy but ad libitum postnatal. Slaughter body weight (20.96 kg); total weight gain (11.34 kg); average daily weight (0.096 kg); fasting carcass weight (19.45); hot carcass weight (9.33 kg); and cold carcass weight (9.11 kg) were smaller (P < 0.05) in POSTN-R but similar between AL-AL and PREN-AL. Rib (0.47 kg), shoulder (0.85 kg), loin (0.50 kg), and chest-flank (0.97 kg) had lower weight in lambs under POSTN-R (P < 0.05). The treatments did not affect the tissue composition of the 12th rib. The non-carcass components and carcass traits are the variables that best discriminate animals under FR plans. The postpartum FR in native lambs has an impact on the performance and carcass weight of the progeny, unlike lambs submitted to prepartum FR, which demonstrates the adaptation to intrauterine nutritional deficiency in Morada Nova ewes to produce lambs with heavy carcasses in regions semi-arid.

Genotype by environment interaction for fat and protein yields via reaction norms in Holstein cattle of southern Brazil.

Our objective was to evaluate the genetic merit of Holstein cattle population in southern Brazil in response to variations in the regional temperature by analyzing the genotype by environment interaction using reaction norms. Fat yield (FY) and protein yield (PY) data of 67 360 primiparous cows were obtained from the database of the Paraná Holstein Breeders Association, Brazil (APCBRH). The regional average annual temperature was used as the environmental variable. A random regression model was adopted applying mixed models with Restricted Maximum Likelihood (REML) algorithm using WOMBAT software. The genetic merit of the 15 most representative bulls, depending on the temperature gradient, was evaluated. Heritability ranged from 0.21 to 0.27 for FY and from 0.14 to 0.20 for PY. The genetic correlation observed among the environmental gradients proved to be higher than 0.80 for both traits. Slight reranking of bulls for both traits was detected, demonstrating that non-relevant genotype by environment interaction for FY and PY were observed. Consequently, no inclusion of the temperature effect in the model of genetic evaluation in southern Brazilian Holstein breed is required.

Genetic parameters for milk yield and quality traits of Brazilian Holstein cows as a function of temperature and humidity index.

Measurements of milk yield (MY), somatic cell score (SCS), percentage of fat (FP), protein (PP), lactose (LP), casein (CP) and percentage of palmitic (C16:0), stearic (C18:0), oleic (C18:1), total saturated (SFA), unsaturated (UFA), monounsaturated (MUFA) and polyunsaturated (PUFA) fatty acids in milk from 5,224 Holstein cows were evaluated as a function of a temperature and humidity index (THI). Legendre orthogonal polynomials from second to seventh order were tested. The best fit order for MY, PP and C18:0 was the third, whereas the second for all other traits. The heritability estimates decreased for MY (0.31 to 0.14), FP (0.28 to 0.16), LP (0.43 to 0.30), SCS (0.14 to 0.09), SFA (0.33 to 0.22) and C16:0 (0.31 to 0.26), whereas increased for CP (0.32 to 0.42), MUFA (0.08 to 0.13), UFA (0.07 to 0.11) and C18:1 (0.07 to 0.11) as the THI level increased. For PP, heritabilities (0.26 to 0.39) presented larger values in intermediate THI. For PUFA and C18:0, heritabilities were approximately constant (0.13 to 0.14 and 0.15, respectively). However, the greatest variations may have been the result of the limitations of Legendre polynomials at the extreme points of the curve, and the pattern of heritabilities curves was approximately constant for the evaluated traits. Spearman's rank correlations between breeding values in extreme THI levels were greater than 0.80 for all traits considering all animals, only cows and only bulls. When considering the top 1% and the top 50% animals (only cows, only bulls and all), Spearman correlations smaller than 0.70 were found, suggesting reranking of the animals. Although there was little variation in the variance components over THI, it is possible that there is no heat stress in the animals studied, because, on average, there was no great impact of the thermal load on the traits. One possible explanation is the use of herds with little climatic difference among herds, as well as the use of fans and sprinklers into the barns. However, the THI levels may be important factors in the selection process, as reranking of animals was verified.

Sheep meat production in the Brazilian semi-arid region: crossing between indigenous breeds.

Investigations in the Brazilian semi-arid region evaluating the performance and carcass traits of sheep of indigenous breeds and their crossings have been performed; however, these studies use exotic breeds which produce precocious lambs with heavier carcasses, but not adapted to climatic conditions and the extensive production system, jeopardizing the sustainability of the sheep production system. We crossed between three indigenous breeds: Morada Nova (MN - maternal breed); Rabo Largo (RL - paternal breed), and Santa Inês (SI - paternal breed) with the objective of evaluating the effect of genotype and sex on the performance and carcass traits of purebred and crossbred animals. A total of 30 lambs, males and females, reared in a semi-intensive system were evaluated. Birth and weaning weights were 2.26 ± 0.53 and 7.31 ± 1.85, respectively. All lambs were slaughtered at 10 months of age. A completely randomized design in a 3 × 2 factorial scheme (three genotypes and two sexes) was used. Multivariate techniques were also performed to reduce group and discriminate variables. Birth and weaning weight were similar (P > 0.05) among genetic groups and sexes. The weight gain, carcass and morphometric trait characteristics, and the main commercial cuts were higher in crossbred lambs (P < 0.05). All indicators have discriminatory power between genotypes and sexes, but the carcass traits have a higher discriminatory power (P < 0.001). All genotypes, regardless of sex, have particular characteristics, i.e. MN × SI was characterized by greater forelimb and ham perimeters (P < 0.001), and the MN × RL by higher hot carcass weight and finish (P = 0.001). The cluster analysis and the heatmap plot revealed associations between SI and the size of cuts and RL with the cut commercial yield and the reduction in weight loss due to cooling. Our findings confirm the hypothesis that crossing between indigenous breeds represents an adequate alternative in sheep meat production systems in semi-arid regions. Finally, we encourage the use of indigenous breeds for sheep meat production with breed identity in order to favor the conservation of genetic resources and the sustainability of the production system.

Genome-wide association scan for QTL and their positional candidate genes associated with internal organ traits in chickens.

BACKGROUND: Poultry breeding programs have been focused on improvement of growth and carcass traits, however, this has resulted in correlated changes in internal organ weights and increased incidence of metabolic disorders. These disorders can affect feed efficiency or even cause death. We used a high density SNP array (600 K, Affymetrix) to estimate genomic heritability, perform genome-wide association analysis, and identify genomic regions and positional candidate genes (PCGs) associated with internal organ traits in an F2 chicken population. We integrated knowledge of haplotype blocks, selection signature regions and sequencing data to refine the list of PCGs. RESULTS: Estimated genomic heritability for internal organ traits in chickens ranged from low (LUNGWT, 0.06) to high (GIZZWT, 0.45). A total of 20 unique 1 Mb windows identified on GGA1, 2, 4, 7, 12, 15, 18, 19, 21, 27 and 28 were significantly associated with intestine length, and weights or percentages of liver, gizzard or lungs. Within these windows, 14 PCGs were identified based on their biological functions: TNFSF11, GTF2F2, SPERT, KCTD4, HTR2A, RB1, PCDH7, LCORL, LDB2, NR4A2, GPD2, PTPN11, ITGB4 and SLC6A4. From those genes, two were located within haplotype blocks and three overlapped with selection signature regions. A total of 13,748 annotated sequence SNPs were in the 14 PCGs, including 156 SNPs in coding regions (124 synonymous, 26 non-synonymous, and 6 splice variants). Seven deleterious SNPs were identified in TNFSF11, NR4A2 or ITGB4 genes. CONCLUSIONS: The results from this study provide novel insights to understand the genetic architecture of internal organ traits in chickens. The QTL detection performed using a high density SNP array covered the whole genome allowing the discovery of novel QTL associated with organ traits. We identified PCGs within the QTL involved in biological processes that may regulate internal organ growth and development. Potential functional genetic variations were identified generating crucial information that, after validation, might be used in poultry breeding programs to reduce the occurrence of metabolic disorders.

Assessing the accuracy of prediction for milk fatty acids by using a small reference population of tropical Holstein cows.

Fatty acids (FA) have been related to effects on human health, sensory quality and shelf life of dairy products, cow's health and methane emission. However, despite their importance, they are not regularly measured in all dairy herds yet, which can affect the accuracy of estimated breeding values (EBV) for these traits. In this case, an alternative is to use genomic selection. Thus, the aim was to assess the use of genomic information in the genetic evaluation for milk traits in a tropical Holstein population. Monthly records (n = 36,457) of milk FA percentage, daily milk yield and quality traits from 4,203 cows as well as the genotypes of 755 of these cows for 57,368 single nucleotide polymorphisms (SNP) were used. Polygenic and genomic-polygenic models were applied for EBV prediction, and both models were compared through the EBV accuracy calculated from the prediction error and Spearman's correlation among EBV rankings. Prediction accuracy was assessed by using cross-validation. In this case, the accuracy was the correlation between the genomic breeding values (GEBV) obtained as the sum of SNP effects and the EBV obtained in the polygenic model in each validation group. For all traits, the use of the genomic-polygenic model did not alter the animals' ranking, with correlations higher than 0.87. Nevertheless, through this model, the accuracy increased from 1.5% to 6.8% compared to the polygenic model. The correlations between GEBV and EBV varied from 0.52 to 0.68. Therefore, the use of a small group of genotyped cows in the genetic evaluation can increase the accuracy of EBV for milk FA and other traditional milk traits.

Exploring candidate genes for heat tolerance in ovine through liver gene expression.

Thermotolerance has become an essential factor in the prevention of the adverse effects of heat stress, but it varies among animals. Identifying genes related to heat adaptability traits is important for improving thermotolerance and for selecting more productive animals in hot environments. The primary objective of this research was to find candidate genes in the liver that play a crucial role in the heat stress response of Santa Ines sheep, which exhibit varying levels of heat tolerance. To achieve this goal, 80 sheep were selected based on their thermotolerance and placed in a climate chamber for 10 days, during which the average temperature was maintained at 36 °C from 10 a.m. to 4 p.m. and 28 °C from 4 p.m. to 10 a.m. A subset of 14 extreme animals, with seven thermotolerant and seven non-thermotolerant animals based on heat loss (rectal temperature), were selected for liver sampling. RNA sequencing and differential gene expression analysis were performed. Thermotolerant sheep showed higher expression of genes GPx3, RGS6, GPAT3, VLDLR, LOC101108817, and EVC. These genes were mainly related to the Hedgehog signaling pathway, glutathione metabolism, glycerolipid metabolism, and thyroid hormone synthesis. These enhanced pathways in thermotolerant animals could potentially mitigate the negative effects of heat stress, conferring greater heat resistance.

Novel putative causal mutations associated with fat traits in Nellore cattle uncovered by eQTLs located in open chromatin regions.

Intramuscular fat (IMF) and backfat thickness (BFT) are critical economic traits impacting meat quality. However, the genetic variants controlling these traits need to be better understood. To advance knowledge in this area, we integrated RNA-seq and single nucleotide polymorphisms (SNPs) identified in genomic and transcriptomic data to generate a linkage disequilibrium filtered panel of 553,581 variants. Expression quantitative trait loci (eQTL) analysis revealed 36,916 cis-eQTLs and 14,408 trans-eQTLs. Association analysis resulted in three eQTLs associated with BFT and 24 with IMF. Functional enrichment analysis of genes regulated by these 27 eQTLs revealed noteworthy pathways that can play a fundamental role in lipid metabolism and fat deposition, such as immune response, cytoskeleton remodeling, iron transport, and phospholipid metabolism. We next used ATAC-Seq assay to identify and overlap eQTL and open chromatin regions. Six eQTLs were in regulatory regions, four in predicted insulators and possible CCCTC-binding factor DNA binding sites, one in an active enhancer region, and the last in a low signal region. Our results provided novel insights into the transcriptional regulation of IMF and BFT, unraveling putative regulatory variants.