Exploring the genetic control of sweat gland characteristics in beef cattle for enhanced heat tolerance.

BACKGROUND: Thermal stress in subtropical regions is a major limiting factor in beef cattle production systems with around $369 million being lost annually due to reduced performance. Heat stress causes numerous physiological and behavioral disturbances including reduced feed intake and decreased production levels. Cattle utilize various physiological mechanisms such as sweating to regulate internal heat. Variation in these traits can help identify genetic variants that control sweat gland properties and subsequently allow for genetic selection of cattle with greater thermotolerance. METHODS: This study used 2,401 Brangus cattle from two commercial ranches in Florida. Precise phenotypes that contribute to an animal's ability to manage heat stress were calculated from skin biopsies and included sweat gland area, sweat gland depth, and sweat gland length. All animals were genotyped with the Bovine GGP F250K, and BLUPF90 software was used to estimate genetic parameters and for Genome Wide Association Study. RESULTS: Sweat gland phenotypes heritability ranged from 0.17 to 0.42 indicating a moderate amount of the phenotypic variation is due to genetics, allowing producers the ability to select for favorable sweat gland properties. A weighted single-step GWAS using sliding 10 kb windows identified multiple quantitative trait loci (QTLs) explaining a significant amount of genetic variation. QTLs located on BTA7 and BTA12 explained over 1.0% of genetic variance and overlap the ADGRV1 and CCDC168 genes, respectively. The variants identified in this study are implicated in processes related to immune function and cellular proliferation which could be relevant to heat management. Breed of Origin Alleles (BOA) were predicted using local ancestry in admixed populations (LAMP-LD), allowing for identification of markers' origin from either Brahman or Angus ancestry. A BOA GWAS was performed to identify regions inherited from particular ancestral breeds that might have a significant impact on sweat gland phenotypes. CONCLUSIONS: The results of the BOA GWAS indicate that both Brahman and Angus alleles contribute positively to sweat gland traits, as evidenced by favorable marker effects observed from both genetic backgrounds. Understanding and utilizing genetic traits that confer better heat tolerance is a proactive approach to managing the impacts of climate change on livestock farming.

Exploring genomic inbreeding and selection signatures in a commercial Brangus herd through functional annotation.

Composite breeds, including Brangus, are widely utilized in subtropical and tropical regions to harness the advantages of both Bos t. taurus and Bos t. indicus breeds. The formation and subsequent selection of composite breeds may result in discernible signatures of selection and shifts in genomic population structure. The objectives of this study were to 1) assess genomic inbreeding, 2) identify signatures of selection, 3) assign functional roles to these signatures in a commercial Brangus herd, and 4) contrast signatures of selection between selected and non-selected cattle from the same year. A total of 4035 commercial Brangus cattle were genotyped using the GGP-F250K array. Runs of Homozygosity (ROH) were used to identify signatures of selection and calculate genomic inbreeding. Quantitative trait loci (QTL) enrichment analysis and literature search identified phenotypic traits linked to ROH islands. Genomic inbreeding averaged 5%, primarily stemming from ancestors five or more generations back. A total of nine ROH islands were identified, QTL enrichment analysis revealed traits related to growth, milk composition, carcass, reproductive, and meat quality traits. Notably, the ROH island on BTA14 encompasses the pleiomorphic adenoma (PLAG1) gene, which has been linked to growth, carcass, and reproductive traits. Moreover, ROH islands associated with milk yield and composition were more pronounced in selected replacement heifers of the population, underscoring the importance of milk traits in cow-calf production. In summary, our research sheds light on the changing genetic landscape of the Brangus breed due to selection pressures and reveals key genomic regions impacting production 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.

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.

Identification of quantitative trait loci and associated candidate genes for pregnancy success in Angus-Brahman crossbred heifers.

BACKGROUND: In beef cattle, more than 50% of the energy input to produce a unit of beef is consumed by the female that produced the calf. Development of genomic tools to identify females with high genetic merit for reproductive function could increase the profitability and sustainability of beef production. RESULTS: Genome-wide association studies (GWAS) were performed using a single-step genomic best linear unbiased prediction approach on pregnancy outcome traits from a population of Angus-Brahman crossbred heifers. Furthermore, a validation GWAS was performed using data from another farm. Heifers were genotyped with the Bovine GGP F250 array that contains 221,077 SNPs. In the discovery population, heifers were bred in winter breeding seasons involving a single round of timed artificial insemination (AI) followed by natural mating for 3 months. Two phenotypes were analyzed: pregnancy outcome to first-service AI (PAI; n = 1,481) and pregnancy status at the end of the breeding season (PEBS; n = 1,725). The heritability was estimated as 0.149 and 0.122 for PAI and PEBS, respectively. In the PAI model, one quantitative trait locus (QTL), located between 52.3 and 52.5 Mb on BTA7, explained about 3% of the genetic variation, in a region containing a cluster of γ-protocadherin genes and SLC25A2. Other QTLs explaining between 0.5% and 1% of the genetic variation were found on BTA12 and 25. In the PEBS model, a large QTL on BTA7 was synonymous with the QTL for PAI, with minor QTLs located on BTA5, 9, 10, 11, 19, and 20. The validation population for pregnancy status at the end of the breeding season were Angus-Brahman crossbred heifers bred by natural mating. In concordance with the discovery population, the large QTL on BTA7 and QTLs on BTA10 and 12 were identified. CONCLUSIONS: In summary, QTLs and candidate SNPs identified were associated with pregnancy outcomes in beef heifers, including a large QTL associated with a group of protocadherin genes. Confirmation of these associations with larger populations could lead to the development of genomic predictions of reproductive function in beef cattle. Moreover, additional research is warranted to study the function of candidate genes associated with QTLs.

Impact of Brahman genetics on skin histology characteristics with implications for heat tolerance in cattle.

Cattle lose heat predominantly through cutaneous evaporation at the skin-hair coat interface when experiencing heat stress. Sweating ability, sweat gland properties, and hair coat properties are a few of the many variables determining the efficacy of evaporative cooling. Sweating is a significant heat dissipation mechanism responsible for 85% of body heat loss when temperatures rise above 86°F. The purpose of this study was to characterize skin morphological parameters in Angus, Brahman, and their crossbred cattle. Skin samples were taken during the summer of 2017 and 2018 from a total of 319 heifers from six breed groups ranging from 100% Angus to 100% Brahman. Epidermis thickness decreased as the percentage of Brahman genetics increased where the 100% Angus group had a significantly thicker epidermis compared to the 100% Brahman animals. A more extended epidermis layer was identified in Brahman animals due to more pronounced undulations in this skin layer. Breed groups with 75% and 100% Brahman genes were similar and had the largest sweat gland area, indicative of superior resilience to heat stress, compared to breed groups with 50% or lower Brahman genetics. There was a significant linear breed group effect on sweat gland area indicating an increase of 862.0 µm2 for every 25% increase in Brahman genetics. Sweat gland length increased as the Brahman percentage increased, while the sweat gland depth showed an opposite trend, decreasing from 100% Angus to 100% Brahman. The number of sebaceous glands was highest in 100% Brahman animals which had about 1.77 more sebaceous glands (p < 0.05) per 4.6 mm2area. Conversely, the sebaceous gland area was greatest in the 100% Angus group. This study identified significant differences in skin properties related to heat exchange ability between Brahman and Angus cattle. Equally important, these differences are also accompanied by significant levels of variation within each breed, which is indicative that selection for these skin traits would improve the heat exchange ability in beef cattle. Further, selecting beef cattle for these skin traits would lead to increased resilience to heat stress without disrupting production traits.

Carcass and meat quality traits in Brangus steers.

The quality grade system used in the United States to identify carcasses with superior eating satisfaction to consumers is based on the amount of marbling within the ribeye and the maturity of the carcass. However, the most important quality attribute for consumers is tenderness. The objective of this study was to investigate the phenotypic correlations between carcass and meat quality traits of strip loin steaks from Brangus steers, particularly the relation between the United States Department of Agriculture (USDA) quality grade and tenderness. Warner-Bratzler shear force (WBSF) values in this study averaged 5.10 ± 0.96 kg, slightly higher than the national average of 4.55 ± 1.14 kg. Average WBSF across all quality grades ranged from 4.90 to 5.27 kg with standard deviations ranging from 0.78 to 1.40 kg. In the present Brangus steer population, there was a weak negative (-0.13) but favorable correlation (P < 0.05) between marbling score and tenderness measured through WBSF. The USDA quality grade had a significant (P = 0.02) effect on WBSF. The WBSF least square means were significantly higher in the Select compared with the Choice¯, and Choiceº quality grades. The Choice⁺ and Prime quality grades were not significantly different from any quality grades regarding the WBSF. The standard quality grade did not have significantly different WBSF least square means from any other quality grade. The range of WBSF values was large, especially in the lower quality grade categories, indicating that there is considerable variation in tenderness, even within quality grade. The high level of variation in tenderness within USDA quality grades highlights the limitation of the USDA grading system to predict eating quality or tenderness.

What a 31-yr multibreed herd taught us about the influence of B. indicus genetics on reproductive performance of cows.

Bos taurus × Bos indicus crosses are widespread in tropical and subtropical regions, nonetheless, quantitative information about the influence of B. indicus genetics on the reproductive performance of beef cattle is lacking. Herein, we determined the association between level of B. indicus genetics and reproduction from a 31-yr dataset comprising sequential breeding seasons of the University of Florida multibreed herd (n = 6,503 Angus × Brahman cows). The proportion of B. indicus genetics in this herd is evenly distributed by each 1/32nd or approximately 3-percentage points. From 1989 to 2020, the estrous cycle of cows was synchronized for artificial insemination (AI) based on detected estrus or timed-AI (TAI) using programs based on gonadotropin-releasing hormone and prostaglandin, and progestin/progesterone. All cows were exposed to natural service after AI and approximately 90-d breeding seasons, considering the day of AI as day 0. The proportion of B. indicus genetics of cows was associated negatively with pregnancy per AI, ranging from 51.6% for cows with 0%-19% of B. indicus genetics to 37.4% for cows with 81%-100% of B. indicus genetics. Similar association was found for estrous response at the end of the synchronization protocol, ranging from 66.3% to 38.4%, respectively. This reduced estrous response helped to explain the pregnancy results, once the pregnancy to AI of cows showing estrus was 2.3-fold greater than for those not showing estrus and submitted to TAI. Despite reduced pregnancy per AI, the increase in the proportion of B. indicus genetics of cows was not associated with a reduction in the proportion of pregnant cows at the end of the breeding season. Nevertheless, the interval from entering the breeding season to pregnancy was lengthened as the proportion of B. indicus genetics of cows increased. The median days to pregnancy was extended by 25 when the proportion of B. indicus genetics surpassed 78% compared with less than 20%. Thus, the increase in the proportion of B. indicus genetics of cows was related to a reduction in pregnancy per AI and lengthening the interval to attain pregnancy during the breeding season, but not with the final proportion of pregnant cows. As a result, reproductive management strategies directed specifically to cows with a greater proportion of B. indicus genetics are needed to improve the rate of pregnancy in beef herds.

Cow­calf operations in the tropics and sub-tropics have benefited from the environmental adaptation provided by Bos indicus genetics. However, reproductive performance has been a cause of concern, although poorly quantified. This study characterized how much the B. indicus genetics in crossbred cows influence herd reproduction. We analyzed data from cows with known proportions of Angus and Brahman genetics, from the same crossbred herd, for 31 sequential breeding seasons. The increase in the proportion of B. indicus genetics reduced estrous response and pregnancy per artificial insemination after estrous synchronization, but not the proportion of pregnant cows at the end of the breeding season. Interval from the beginning of the breeding season to pregnancy was extended by 25 d when the proportion of B. indicus genetics surpassed 78%. In conclusion, reproductive management strategies directed specifically to cows with a greater proportion of B. indicus genetics are needed to improve the rate of pregnancy in beef herds.

Genomic Approaches Reveal Pleiotropic Effects in Crossbred Beef Cattle.

Carcass and meat quality are two important attributes for the beef industry because they drive profitability and consumer demand. These traits are of even greater importance in crossbred cattle used in subtropical and tropical regions for their superior adaptability because they tend to underperform compared to their purebred counterparts. Many of these traits are challenging and expensive to measure and unavailable until late in life or after the animal is harvested, hence unrealistic to improve through traditional phenotypic selection, but perfect candidates for genomic selection. Before genomic selection can be implemented in crossbred populations, it is important to explore if pleiotropic effects exist between carcass and meat quality traits. Therefore, the objective of this study was to identify genomic regions with pleiotropic effects on carcass and meat quality traits in a multibreed Angus-Brahman population that included purebred and crossbred animals. Data included phenotypes for 10 carcass and meat quality traits from 2,384 steers, of which 1,038 were genotyped with the GGP Bovine F-250. Single-trait genome-wide association studies were first used to investigate the relevance of direct additive genetic effects on each carcass, sensory and visual meat quality traits. A second analysis for each trait included all other phenotypes as covariates to correct for direct causal effects from identified genomic regions with pure direct effects on the trait under analysis. Five genomic windows on chromosomes BTA5, BTA7, BTA18, and BTA29 explained more than 1% of additive genetic variance of two or more traits. Moreover, three suggestive pleiotropic regions were identified on BTA10 and BTA19. The 317 genes uncovered in pleiotropic regions included anchoring and cytoskeletal proteins, key players in cell growth, muscle development, lipid metabolism and fat deposition, and important factors in muscle proteolysis. A functional analysis of these genes revealed GO terms directly related to carcass quality, meat quality, and tenderness in beef cattle, including calcium-related processes, cell signaling, and modulation of cell-cell adhesion. These results contribute with novel information about the complex genetic architecture and pleiotropic effects of carcass and meat quality traits in crossbred beef cattle.

Genome-wide scan reveals important additive and non-additive genetic effects associated with resistance to Haemonchus contortus in Florida Native sheep.

Florida Native sheep is among the sheep breeds best adapted to humid and hot climatic conditions such as those of Florida, USA, and have shown a superior ability to regulate nematode burdens. This is one of the oldest sheep breeds in North America and is an endangered species. To ensure genetic diversity and long-term survival of the breed, protection of the current genetic stock is critical and conservation efforts are required to promote its breeding and production. The objective of the present study was to investigate the importance of additive and non-additive genetic effects on resistance to natural Haemonchus contortus infections in Florida Native sheep using a whole genome scan. A total of 200 sheep were evaluated in the present study. Phenotypic records included faecal egg count (FEC, eggs/gram), FAMACHA® score, packed cell volume (PCV, %), body condition score and average daily gain (ADG, kg). Sheep were genotyped using the GGP Ovine 50K SNP chip and 45.2 k single nucleotide polymorphism (SNP) markers spanning the entire genome were available for quality control procedures. Mixed models were used to analyse the response variables and included the identity by state matrix to control for population structure. Bonferroni correction was used to control for multiple testing and a second arbitrary threshold (0.1 × 10-3) was used. Fifteen SNPs with additive and non-additive genetic effects and located in Ovis aries chromosome OAR1, 2, 3, 6, 8, 10, 11, 12, 13 and 21 were associated with FEC, FAMACHA® score, PCV and ADG. These SNPs could be potential genetic markers for resistance to natural H. contortus exposure in Florida Native sheep.

RNA-seq analysis identifies cytoskeletal structural genes and pathways for meat quality in beef.

RNA sequencing (RNA-seq) has allowed for transcriptional profiling of biological systems through the identification of differentially expressed (DE) genes and pathways. A total of 80 steers with extreme phenotypes were selected from the University of Florida multibreed Angus-Brahman herd. The average slaughter age was 12.91±8.69 months. Tenderness, juiciness and connective tissue assessed by sensory panel, along with marbling, Warner-Bratzler Shear Force (WBSF) and cooking loss, were measured in longissimus dorsi muscle. Total RNA was extracted from muscle and one RNA-seq library per sample was constructed, multiplexed, and sequenced based on protocols by Illumina HiSeq-3000 platform to generate 2×101 bp paired-end reads. The overall read mapping rate using the Btau_4.6.1 reference genome was 63%. A total of 8,799 genes were analyzed using two different methodologies, an expression association and a DE analysis. A gene and exon expression association analysis was carried out using a meat quality index on all 80 samples as a continuous response variable. The expression of 208 genes and 3,280 exons from 1,565 genes was associated with the meat quality index (p-value ≤ 0.05). A gene and isoform DE evaluation was performed analyzing two groups with extreme WBSF, tenderness and marbling. A total of 676 (adjusted p-value≤0.05), 70 (adjusted p-value≤0.1) and 198 (adjusted p-value≤0.1) genes were DE for WBSF, tenderness and marbling, respectively. A total of 106 isoforms from 98 genes for WBSF, 13 isoforms from 13 genes for tenderness and 43 isoforms from 42 genes for marbling (FDR≤0.1) were DE. Cytoskeletal and transmembrane anchoring genes and pathways were identified in the expression association, DE and the gene enrichment analyses; these proteins can have a direct effect on meat quality. Cytoskeletal proteins and transmembrane anchoring molecules can influence meat quality by allowing cytoskeletal interaction with myocyte and organelle membranes, contributing to cytoskeletal structure and architecture maintenance postmortem.

Whole Genome Sequence Data Provides Novel Insights Into the Genetic Architecture of Meat Quality Traits in Beef.

Tenderness is a major quality attribute for fresh beef steaks in the United States, and meat quality traits in general are suitable candidates for genomic research. The objectives of the present analysis were to (1) perform genome-wide association (GWA) analysis for marbling, Warner-Bratzler shear force (WBSF), tenderness, and connective tissue using whole-genome data in an Angus population, (2) identify enriched pathways in each GWA analysis; (3) construct a protein-protein interaction network using the associated genes and (4) perform a µ-calpain proteolysis assessment for associated structural proteins. An Angus-sired population of 2,285 individuals was assessed. Animals were transported to a commercial packing plant and harvested at an average age of 457 ± 46 days. After 48 h postmortem, marbling was recorded by graders' visual appraisal. Two 2.54-cm steaks were sampled from each muscle for recording of WBSF, and tenderness, and connective tissue by a sensory panel. The relevance of additive effects on marbling, WBSF, tenderness, and connective tissue was evaluated on a genome-wide scale using a two-step mixed model-based approach in single-trait analysis. A tissue-restricted gene enrichment was performed for each GWA where all polymorphisms with an association p-value lower than 1 × 10-3 were included. The genes identified as associated were included in a protein-protein interaction network and a candidate structural protein assessment of proteolysis analyses. A total of 1,867, 3,181, 3,926, and 3,678 polymorphisms were significantly associated with marbling, WBSF, tenderness, and connective tissue, respectively. The associate region on BTA29 (36,432,655-44,313,046 bp) harbors 13 highly significant markers for meat quality traits. Enrichment for the GO term GO:0005634 (Nucleus), which includes transcription factors, was evident. The final protein-protein network included 431 interations between 349 genes. The 42 most important genes based on significance that encode structural proteins were included in a proteolysis analysis, and 81% of these proteins were potential µ-Calpain substrates. Overall, this comprehensive study unraveled genetic variants, genes and mechanisms of action responsible for the variation in meat quality traits. Our findings can provide opportunities for improving meat quality in beef cattle via marker-assisted selection.

The effect of Brahman genes on body temperature plasticity of heifers on pasture under heat stress.

Bos taurus indicus cattle have the superior ability for the regulation of body temperature during heat stress due to a number of physiological and cellular level adaptive traits. The objectives of this study were to quantify the change in body temperature in heifers with various proportions of Brahman genes per unit increase in heat stress as measured by temperature-humidity index (THI) and to assess how different breed groups responded to varying intensity and duration of heat stress. A total of 299 two-yr-old heifers from six breed groups ranging from 100% Angus to 100% Brahman were evaluated under hot and humid conditions during 2017 and 2018 summer days. Two strategies were used to estimate the plasticity in body temperature of breed groups in response to environmental challenges: 1) a random regression mixed model was used to estimate reaction norm parameters for each breed group in response to a specified environmental heat stress and 2) a repeated measures mixed model was used to evaluate the response to different environmental heat loads. The reaction norm model estimated an intercept and slope measuring the change in body temperature per unit increase in THI environmental heat stress for different breed groups of animals and allowed the identification of genotypes which are robust, with low slope values indicative of animals that are able to maintain normal body temperature across a range of environments. The repeated measures mixed model showed that Brahman cattle have an advantage under moderate or high heat stress conditions but both Angus and Brahman breed groups are greatly affected when heat stress is severe. A critical factor appears to be the opportunity to cool down during the night hours more than the number of hours with extreme THI. With heat stress conditions predicted to intensify and expand into currently temperate zones, developing effective strategies to ensure sustainable beef production systems are imperative. Effective strategies will require the identification of the genes conferring the superior thermotolerance in Brahman cattle.

Identification of eQTLs and sQTLs associated with meat quality in beef.

BACKGROUND: Transcription has a substantial genetic control and genetic dissection of gene expression could help us understand the genetic architecture of complex phenotypes such as meat quality in cattle. The objectives of the present research were: 1) to perform eQTL and sQTL mapping analyses for meat quality traits in longissimus dorsi muscle; 2) to uncover genes whose expression is influenced by local or distant genetic variation; 3) to identify expression and splicing hot spots; and 4) to uncover genomic regions affecting the expression of multiple genes. RESULTS: Eighty steers were selected for phenotyping, genotyping and RNA-seq evaluation. A panel of traits related to meat quality was recorded in longissimus dorsi muscle. Information on 112,042 SNPs and expression data on 8588 autosomal genes and 87,770 exons from 8467 genes were included in an expression and splicing quantitative trait loci (QTL) mapping (eQTL and sQTL, respectively). A gene, exon and isoform differential expression analysis previously carried out in this population identified 1352 genes, referred to as DEG, as explaining part of the variability associated with meat quality traits. The eQTL and sQTL mapping was performed using a linear regression model in the R package Matrix eQTL. Genotype and year of birth were included as fixed effects, and population structure was accounted for by including as a covariate the first PC from a PCA analysis on genotypic data. The identified QTLs were classified as cis or trans using 1 Mb as the maximum distance between the associated SNP and the gene being analyzed. A total of 8377 eQTLs were identified, including 75.6% trans, 10.4% cis, 12.5% DEG trans and 1.5% DEG cis; while 11,929 sQTLs were uncovered: 66.1% trans, 16.9% DEG trans, 14% cis and 3% DEG cis. Twenty-seven expression master regulators and 13 splicing master regulators were identified and were classified as membrane-associated or cytoskeletal proteins, transcription factors or DNA methylases. These genes could control the expression of other genes through cell signaling or by a direct transcriptional activation/repression mechanism. CONCLUSION: In the present analysis, we show that eQTL and sQTL mapping makes possible positional identification of gene and isoform expression regulators.

Characterization of water intake and water efficiency in beef cattle1,2.

In the future, water may not be as readily available due to increases in competition from a growing human population, wildlife, and other agricultural sectors, making selection for water efficiency of beef cattle increasingly important. Substantial selection emphasis has recently been placed on feed efficiency in an effort to reduce production costs, but no emphasis has been placed on making cattle more water efficient due to lack of data. Thus, the objective of this study was to calculate water efficiency metrics for cattle and evaluate their relationship to growth, feed intake (FI), and feed efficiency. Individual daily FI and water intake (WI) records were collected on 578 crossbred steers over a 70-d test period. Animals with low water intake ate less feed, had lower gains, and were more water efficient (as defined by water to gain ratio, W/G, and residual water intake, RWI). However, the amount of water consumed by animals had minimal phenotypic relationship with feed efficiency (residual feed intake [RFI], R2 = 0.1050 and feed to gain ratio (F/G) ratio R2 = 0.0726). Cattle that had low DMI consumed less water, had lower gains, had lower RFI, and had higher F/G. The level of feed consumed had minimal relationship with water efficiency. WI, W/G, RWI, and ADG had moderate heritability estimates of 0.39, 0.39, 0.37, and 0.37, respectively. High heritability estimates were observed for DMI and RFI (0.67 and 0.65, respectively). Feed to gain had a low heritability estimate of 0.16. WI had a strong positive genetic correlation with W/G (0.99) and RWI (0.88), thus selecting for decreased WI should also make cattle more water efficient. The genetic correlation between WI and ADG was 0.05; thus, selecting for low WI cattle should have little effect on growth. There is a low to moderate genetic correlation between WI and DMI (0.34). RWI has a positive genetic correlation with W/G ratio (0.89) and F/G ratio (0.42) and is negatively genetically correlated with RFI (-0.57). Water to gain and F/G had a strong positive genetic correlation (0.68). RFI has a positive genetic correlation with W/G ratio (0.37) and F/G (0.88). Minimal antagonisms seem to be present between WI and ADG, although it should be noted that standard errors were large and often not significantly different from zero due to the small sample size. However, care should be taken to ensure that unintended changes do not occur in DMI or other production traits and incorporation of WI into a selection index would likely prove to be the most effective method for selection.

Signatures of selection for resistance to Haemonchus contortus in sheep and goats.

BACKGROUND: Gastrointestinal nematode infection (GNI) is the most important disease affecting the small ruminant industry in U.S. The environmental conditions in the southern United States are ideal for the survival of the most pathogenic gastrointestinal nematode, Haemonchus contortus. Host genetic variation for resistance to H. contortus allows selective breeding for increased resistance of animals. This selection process increases the prevalence of particular alleles in sheep and goats and creates unique genetic patterns in the genome of these species. The aim of this study was to identify loci with divergent allelic frequencies in a candidate gene panel of 100 genes using two different approaches (frequentist and Bayesian) to estimate Fst outliers in three different breeds of sheep and goats exposed to H. contortus. RESULTS: Our results for sheep populations showed SNPs under selection in C3AR1, CSF3, SOCS2, NOS2, STAT5B, TGFB2 and IL2RA genes using frequentist and Bayesian approaches. For goats, SNPs in CD1D, ITGA9, IL12A, IL13RA1, CD86 and TGFB2 genes were under selection. Common signatures of selection in both species were observed in NOS2, TGFB2 and TLR4 genes. Directional selection was present in all SNPs evaluated in the present study. CONCLUSIONS: A total of 13 SNPs within 7 genes of our candidate gene panel related to H. contortus exposure were identified under selection in sheep populations. For goats, 11 SNPs within 7 genes were identified under selection. Results from this study support the hypothesis that resistance to H. contortus is likely to be controlled by many loci. Shared signatures of selection related to mechanisms of immune protection against H. contortus infection in sheep and goats could be useful targets in breeding programs aimed to produce resistant animals with low FEC.

Consumer preferences for beef with improved nutrient profile1.

Although beef is a nutrient-rich foodstuff excelling in protein, vitamins, and minerals, there is controversy regarding the fat content of beef and its healthfulness in the diet. Although much of the fatty acid (FA) content in beef is considered "healthy fats," many consumers are confused about the different classifications of FA. The objectives of this study were to determine consumers' knowledge about the nutritional value of beef and its importance in purchasing decisions, and to gain a better understanding about preferences for changes in FA composition. Objectives of the study were completed through 2 consumer studies: 1) an online survey and 2) a taste-panel auction. In the online survey, respondents were asked to choose between 2 steaks that varied in polyunsaturated and saturated FA levels, iron content, and price. Respondents were also asked to categorize "Monounsaturated Fat," "Polyunsaturated Fat," "Saturated Fat," and "Trans Fat," as either "healthy" or "unhealthy" both before and after an educational excerpt was provided. The results from the online survey indicated many consumers are unclear about the differences in beef nutritional value, specifically FA content. Initially, only 66.4%, 69.1%, 79.1%, and 79.2% of respondents correctly categorized the monounsaturated, polyunsaturated, saturated, and trans fat, respectively. However, more than 90% of respondents correctly categorized the various FAs after an educational excerpt was provided. After survey respondents better understood the healthfulness of FA in beef, they were also willing to pay a premium for a steak with improved FA composition. However, these premiums diminished when participants had to actually put forth a monetary value for a steak in the taste-panel auction. Research shows that there is variation among cattle for FA composition. This provides opportunity to identify cattle with a favorable composition and market this product to the increasing population of health-conscious consumers. Our results provide insight for beef promotion and marketing opportunities and indicate that relaying information about FA content is extremely important to collect a premium for healthier beef.

Association study reveals Th17, Treg, and Th2 loci related to resistance to Haemonchus contortus in Florida Native sheep1.

The aim of this study was to identify for the first time single nucleotide polymorphisms (SNPs) associated with Haemonchus contortus resistance in Florida Native sheep, using a targeted sequencing approach. One hundred and fifty-three lambs were evaluated in this study. At the start of the trial, phenotypic records for fecal egg count (FEC), FAMACHA score, body condition score (BCS), and weight were recorded and deworming of sheep with levamisole (18 mg/kg of body weight) was performed. Ten days post-deworming (baseline) and 28 d post-baseline, a full hematogram of each sheep was obtained and FEC, FAMACHA score, BCS, and weight were assessed. Average daily gain was calculated at the end of the trial. Out of 153 animals, 100 sheep were selected for genotyping using a targeted sequencing approach. Targeted sequencing panel included 100 candidate genes for immune response against H. contortus. SNPs were discarded if call rate <95% and minor allele frequency ≤0.05. A mixed model was used to analyze the response variables and included the identity by state matrix to control for population structure. A contemporary group (age, group, and sex) was included as fixed effect. Bonferroni correction was used to control for multiple testing. Eighteen SNPs on chromosomes 1, 2, 3, 4, 6, 7, 11, 15, 18, 20, 24, and 26 were significant for different traits. Our results suggest that loci related to Th17, Treg, and Th2 responses play an important role in the expression of resistant phenotypes. Several genes including ITGA4, MUC15, TLR3, PCDH7, CFI, CXCL10, TNF, CCL26, STAT3, GPX2, IL2RB, and STAT6 were identified as potential markers for resistance to natural H. contortus exposure. This is the first study that evaluates potential genetic markers for H. contortus resistance in Florida Native sheep.

Genetic parameters for hair characteristics and core body temperature in a multibreed Brahman-Angus herd1.

Thermal stress in hot humid conditions limits cattle production. The objectives for this study were to estimate genetic parameters for hair characteristics and core body temperature under low and high temperature humidity index (THI) conditions. Hair samples were collected and measured for length and diameter. Core body temperature was measured as vaginal temperature every 15 min over a 5-d period using an iButton temperature measuring device implanted in a blank CIDR in 336 heifers from the University of Florida multibreed herd (ranging from 100% Angus to 100% Brahman). Restricted maximum likelihood procedures were used to estimate heritabilities from multiple bivariate animal models using the WOMBAT program. Estimates of heritability for hair diameter, undercoat length, topcoat length, body temperature under low THI conditions, and body temperature under high THI conditions were 0.50, 0.67, 0.42, 0.32, and 0.26, respectively. The genetic parameters estimated in this study indicate a large, exploitable genetic variance which can be selected upon to improve tolerance in cattle. Breed effects for differing compositions of Brahman and Angus were also estimated. As Brahman breed composition increased by 25% undercoat length, topcoat length, body temperature under low THI conditions, and body temperature under high THI conditions decreased by 1.32 mm, 2.94 mm, 0.11 °C, and 0.14 °C, respectively. Under both low and high THI conditions, cattle with 25% Brahman breed composition or greater maintained a significantly lower body temperature than the 100% Angus breed group. The incorporation of Brahman germplasm is recommended for herds that often experience heat stress conditions in order to increase resilience to heat stress.

Genome wide association and gene enrichment analysis reveal membrane anchoring and structural proteins associated with meat quality in beef.

BACKGROUND: Meat quality related phenotypes are difficult and expensive to measure and predict but are ideal candidates for genomic selection if genetic markers that account for a worthwhile proportion of the phenotypic variation can be identified. The objectives of this study were: 1) to perform genome wide association analyses for Warner-Bratzler Shear Force (WBSF), marbling, cooking loss, tenderness, juiciness, connective tissue and flavor; 2) to determine enriched pathways present in each genome wide association analysis; and 3) to identify potential candidate genes with multiple quantitative trait loci (QTL) associated with meat quality. RESULTS: The WBSF, marbling and cooking loss traits were measured in longissimus dorsi muscle from 672 steers. Out of these, 495 animals were used to measure tenderness, juiciness, connective tissue and flavor by a sensory panel. All animals were genotyped for 221,077 markers and included in a genome wide association analysis. A total number of 68 genomic regions covering 52 genes were identified using the whole genome association approach; 48% of these genes encode transmembrane proteins or membrane associated molecules. Two enrichment analysis were performed: a tissue restricted gene enrichment applying a correlation analysis between raw associated single nucleotide polymorphisms (SNPs) by trait, and a functional classification analysis performed using the DAVID Bioinformatic Resources 6.8 server. The tissue restricted gene enrichment approach identified eleven pathways including "Endoplasmic reticulum membrane" that influenced multiple traits simultaneously. The DAVID functional classification analysis uncovered eleven clusters related to transmembrane or structural proteins. A gene network was constructed where the number of raw associated uncorrelated SNPs for each gene across all traits was used as a weight. A multiple SNP association analysis was performed for the top five most connected genes in the gene-trait network. The gene network identified the EVC2, ANXA10 and PKHD1 genes as potentially harboring multiple QTLs. Polymorphisms identified in structural proteins can modulate two different processes with direct effect on meat quality: in vivo myocyte cytoskeletal organization and postmortem proteolysis. CONCLUSION: The main result from the present analysis is the uncovering of several candidate genes associated with meat quality that have structural function in the skeletal muscle.