Genome-wide association study and ancestral origins of the slick-hair coat in tropically adapted cattle

Genomic patterns of selection in morphometric traits across diverse Indian cattle breeds

This study seeks a comprehensive exploration of genome-wide selective processes impacting morphometric traits across diverse cattle breeds, utilizing an array of statistical methods. Morphometric traits, encompassing both qualitative and quantitative variables, play a pivotal role in characterizing and selecting livestock breeds based on their external appearance, size, and physical attributes. While qualitative traits, such as color, horn structure, and coat type, contribute to adaptive features and breed identification, quantitative traits like body weight and conformation measurements bear a closer correlation with production characteristics. This study employs advanced genotyping technologies, including the Illumina BovineSNP50 Bead Chip and next-generation sequencing methods like Reduced Representation sequencing, to identify genomic signatures associated with these traits. We applied four intra-population methods to find evidence of selection, such as Tajima's D, CLR, iHS, and ROH. We found a total of 40 genes under the selection signature, that were associated with morphometric traits in five cattle breeds (Kankrej, Tharparkar, Nelore, Sahiwal, and Gir). Crucial genes such as ADIPDQ, DPP6, INSIG1, SLC35D2 in Kankrej, LPL, ATP6V1B2, CDC14B in Tharparkar, HPSE2, PLAG1 in Nelore, PCSK1, PRKD1 in Sahiwal, and GNAQ, HPCAL1 in Gir were identified in our study. This approach provides valuable insights into the genetic basis of variations in body weight and conformation traits, facilitating informed selection processes and offering a deeper understanding of the evolutionary and domestication processes in diverse cattle breeds.

Genome-Wide Analysis of Genetic Diversity and Selection Signatures in Zaobei Beef Cattle

This investigation provides a comprehensive analysis of genomic diversity and selection signatures in Zaobei beef cattle, an indigenous breed known for its adaptation to hot and humid climates and superior meat quality. Whole-genome resequencing was conducted on 23 Zaobei cattle, compared with 46 Simmental cattle to highlight genetic distinctions. Population structure analysis confirmed the genetic uniqueness of Zaobei cattle. Using methods such as DASDC v1.01, XPEHH, and θπ ratio, we identified 230, 232, and 221 genes through DASDC, including hard sweeps, soft sweeps, and linkage sweeps, respectively. Coincidentally, 109 genes were identified when using XPEHH and θπ ratio methods. Together, these analyses revealed eight positive selection genes (ARHGAP15, ZNF618, USH2A, PDZRN4, SPATA6, ROR2, KCNIP3, and VWA3B), which are linked to critical traits such as heat stress adaptation, fertility, and meat quality. Moreover, functional enrichment analyses showed pathways related to autophagy, immune response, energy metabolism, and muscle development. The comprehensive genomic insights gained from this study provide valuable knowledge for breeding programs aimed at enhancing the beneficial traits in Zaobei cattle.

The prolactin receptor: A cross-species comparison of gene structure, transcriptional regulation, tissue-specificity, and genetic variation

The conserved and multifaceted functions of prolactin (PRL) are coordinated through varied distribution and expression of its cell-surface receptor (PRLR) across a range of tissues and physiological states. The resultant heterogeneous expression of PRLR mRNA and protein across different organs and cell types supports a wide range of PRL-regulated processes including reproduction, lactation, development, and homeostasis. Genetic variation within the PRLR gene also accounts for several phenotypes impacting agricultural production and human pathology. The goal of this review is to highlight the many elements that control differential expression of the PRLR across tissues, and the various phenotypes that exist across species due to variation in the PRLR gene.

Selective signatures in composite MONTANA TROPICAL beef cattle reveal potential genomic regions for tropical adaptation

Genomic regions related to tropical adaptability are of paramount importance for animal breeding nowadays, especially in the context of global climate change. Moreover, understanding the genomic architecture of these regions may be very relevant for aiding breeding programs in choosing the best selection scheme for tropical adaptation and/or implementing a crossbreeding scheme. The composite MONTANA TROPICAL® population was developed by crossing cattle of four different biological types to improve production in harsh environments. Pedigree and genotype data (51962 SNPs) from 3215 MONTANA TROPICAL® cattle were used to i) characterize the population structure; ii) identify signatures of selection with complementary approaches, i.e. Integrated Haplotype Score (iHS) and Runs of Homozygosity (ROH); and iii) understand genes and traits related to each selected region. The population structure based on principal components had a weak relationship with the genetic contribution of the different biological types. Clustering analyses (ADMIXTURE) showed different clusters according to the number of generations within the composite population. Considering results of both selection signatures approaches, we identified only one consensus region on chromosome 20 (35399405-40329703 bp). Genes in this region are related to immune function, regulation of epithelial cell differentiation, and cell response to ionizing radiation. This region harbors the slick locus which is related to slick hair and epidermis anatomy, both of which are related to heat stress adaptation. Also, QTLs in this region were related to feed intake, milk yield, mastitis, reproduction, and slick hair coat. The signatures of selection detected here arose in a few generations after crossbreeding between contrasting breeds. Therefore, it shows how important this genomic region may be for these animals to thrive in tropical conditions. Further investigations on sequencing this region can identify candidate genes for animal breeding and/or gene editing to tackle the challenges of climate change.

Heat stress effects on milk yield traits and metabolites and mitigation strategies for dairy cattle breeds reared in tropical and sub-tropical countries

Heat stress is an important problem for dairy industry in many parts of the world owing to its adverse effects on productivity and profitability. Heat stress in dairy cattle is caused by an increase in core body temperature, which affects the fat production in the mammary gland. It reduces milk yield, dry matter intake, and alters the milk composition, such as fat, protein, lactose, and solids-not-fats percentages among others. Understanding the biological mechanisms of climatic adaptation, identifying and exploring signatures of selection, genomic diversity and identification of candidate genes for heat tolerance within indicine and taurine dairy breeds is an important progression toward breeding better dairy cattle adapted to changing climatic conditions of the tropics. Identifying breeds that are heat tolerant and their use in genetic improvement programs is crucial for improving dairy cattle productivity and profitability in the tropics. Genetic improvement for heat tolerance requires availability of genetic parameters, but these genetic parameters are currently missing in many tropical countries. In this article, we reviewed the HS effects on dairy cattle with regard to (1) physiological parameters; (2) milk yield and composition traits; and (3) milk and blood metabolites for dairy cattle reared in tropical countries. In addition, mitigation strategies such as physical modification of environment, nutritional, and genetic development of heat tolerant dairy cattle to prevent the adverse effects of HS on dairy cattle are discussed. In tropical climates, a more and cost-effective strategy to overcome HS effects is to genetically select more adaptable and heat tolerant breeds, use of crossbred animals for milk production, i.e., crosses between indicine breeds such as Gir, white fulani, N'Dama, Sahiwal or Boran to taurine breeds such as Holstein-Friesian, Jersey or Brown Swiss. The results of this review will contribute to policy formulations with regard to strategies for mitigating the effects of HS on dairy cattle in tropical countries.

Identification of Runs of Homozygosity Islands and Genomic Estimated Inbreeding Values in Caqueteño Creole Cattle (Colombia)

The Caqueteño Creole (CAQ) is a native breed of cattle from the Caquetá department (Colombia), adapted to tropical conditions, which is extremely important to production systems in those regions. However, CAQ is poorly studied. In this sense, population structure studies associated with runs of homozygosity (ROH) analysis would allow for a better understanding of CAQ. Through ROH analysis, it is possible to reveal genetic relationships between individuals, measure genome inbreeding levels, and identify regions associated with traits of economic interest. Samples from a CAQ population (n = 127) were genotyped with the Bovine HD BeadChip (777,000 SNPs) and analyzed with the PLINK 1.9 program to estimate FROH and ROH islands. We highlighted a decrease in inbreeding frequency for FROH 4−8 Mb, 8−16 Mb, and >16 Mb classes, indicating inbreeding control in recent matings. We also found genomic hotspot regions on chromosomes 3, 5, 6, 8, 16, 20, and 22, where chromosome 20 harbored four hotspots. Genes in those regions were associated with fertility and immunity traits, muscle development, and environmental resistance, which may be present in the CAQ breed due to natural selection. This indicates potential for production systems in tropical regions. However, further studies are necessary to elucidate the CAQ production objective.

Genetics of Base Coat Colour Variations and Coat Colour-Patterns of the South African Nguni Cattle Investigated Using High-Density SNP Genotypes

Nguni cattle are a Sanga type breed with mixed B. taurus and B. indicus ancestry and proven resistance to ticks, diseases and other harsh conditions of the African geographical landscape. The multi-coloured Nguni coats have found a niche market in the leather industry leading to breeding objectives towards the promotion of such diversity. However, there is limited studies on the genomic architecture underlying the coat colour and patterns hampering any potential breeding and improvement of such trait. This study investigated the genetics of base coat colour, colour-sidedness and the white forehead stripe in Nguni cattle using coat colour phenotyped Nguni cattle and Illumina Bovine HD (770K) genotypes. Base coat colour phenotypes were categorised into eumelanin (n = 45) and pheomelanin (n = 19). Animals were categorised into either colour-sided (n = 46) or non-colour-sided (n = 94) and similarly into presence (n = 15) or absence (n = 67) of white forehead stripe. Genome-wide association tests were conducted using 622,103 quality controlled SNPs and the Efficient Mixed Model Association eXpedited method (EMMAX) implemented in Golden Helix SNP Variation Suite. The genome-wide association studies for base coat colour (eumelanin vs. pheomelanin) resulted into four indicative SNPs on BTA18 and a well-known gene, MC1R, was observed within 1 MB from the indicative SNPs (p < 0.00001) and found to play a role in the melanogenesis (core pathway for melanin production) and the MAPK signalling pathway. GWAS for colour-sidedness resulted in four indicative SNPs, none of which were in close proximity to the KIT candidate gene known for colour-sidedness. GWAS for the white forehead stripe resulted in 17 indicative SNPs on BTA6. Four genes MAPK10, EFNA5, PPP2R3C and PAK1 were found to be associated with the white forehead stripe and were part of the MAPK, adrenergic and Wnt signalling pathways that are synergistically associated with the synthesis of melanin. Overall, our results prove prior knowledge of the role of MC1R in base coat colours in cattle and suggested a different genetic mechanism for forehead stripe phenotypes in Nguni cattle.

Long-term artificial selection of Hanwoo (Korean) cattle left genetic signatures for the breeding traits and has altered the genomic structure

Indigenous Korean breeds such as Hanwoo (Korean) cattle have adapted to their local environment during the past 5000 years. In the 1980s, the National Genetic Improvement Program was established to develop a modern economic breed for beef production in Korea through artificial selection. This process is thought to have altered the genomic structure of breeding traits over time. The detection of genetic variants under selection could help to elucidate the genetic mechanism of artificial selection in modern cattle breeds. Indigenous Hanwoo cattle have adapted in response to local natural and artificial selection during a 40-year breeding program. We analyzed genomic changes in the selection signatures of an unselected population (USP; n = 362) and a selected population (KPN; n = 667) of Hanwoo cattle. Genomic changes due to long-term artificial selection were identified using a genome-wide integrated haplotype score (iHS) and a genome-wide association study (GWAS). Signatures of recent selection were detected as positive (piHS > 6) or negative (piHS < –6) iHS scores spanning more than 46 related genes in KPN cattle, but none in USP cattle. A region adjacent to the PLAG1 gene was found to be under strong selection for carcass weight. The GWAS results also showed a selection signature on BTA14, but none on BTA13. Pathway and quantitative trait locus analysis results identified candidate genes related to energy metabolism, feed efficiency, and reproductive traits in Hanwoo cattle. Strong selection significantly altered Hanwoo cattle genome structural properties such as linkage disequilibrium (LD) and haplotypes through causal mutation for target traits. Haplotype changes of genome structure which are changes of ancestral allele to derived alleles due to selection were clearly identified on BTA13 and BTA14; however, the structure of the LD block was not clearly observed except BTA14. Thus, selection based on EBVs would be working very well in Hanwoo cattle breeding program appears to have been highly successful.

Coat score. A possible explanation for the zebuine selective sweep located on bovine chromosome 5: 47,670,001-48,100,000 bp

Over 65% of the world's cattle population resides in warm areas where heat stress conditions limit the breed of European taurine cattle. Composite breeds were developed to retain the main traits of both parental breeds. The skin plays a central role in animal response to heat stress. Research on the genetic architecture of skin traits has identified genes and regions related to warm resistance skin features. The aim of this study was to determine whether the indicine proportion accounted for coat type or whether there were genes of large effect segregating in Brangus. Bulls (n = 108) were genotyped using microarrays and their coat score and hair length were evaluated. Indicine-taurine genome-wide composition was estimated and GWAS was performed. Although significant correlations between indicine proportion and traits were not observed, four windows of SNPs on BTA4 and BTA5 explained more than 2% of the trait variance. The GWAS for coat score in summer showed the main peak on BTA5:46,941,446-48,030,219 bp, accounting for 4.65% of the variance. Our results suggest that the variation in coat score and undercoat hair length in Argentinian Brangus bulls is associated with the presence of some particular gene variants, rather than with the whole indicine genetic content.

Haplotype-Based Single-Step GWAS for Yearling Temperament in American Angus Cattle

Behavior is a complex trait and, therefore, understanding its genetic architecture is paramount for the development of effective breeding strategies. The objective of this study was to perform traditional and weighted single-step genome-wide association studies (ssGWAS and WssGWAS, respectively) for yearling temperament (YT) in North American Angus cattle using haplotypes. Approximately 266 K YT records and 70 K animals genotyped using a 50 K single nucleotide polymorphisms (SNP) panel were used. Linkage disequilibrium thresholds (LD) of 0.15, 0.50, and 0.80 were used to create the haploblocks, and the inclusion of non-LD-clustered SNPs (NCSNP) with the haplotypes in the genomic models was also evaluated. WssGWAS did not perform better than ssGWAS. Cattle YT was found to be a highly polygenic trait, with genes and quantitative trait loci (QTL) broadly distributed across the whole genome. Association studies using LD-based haplotypes should include NCSNPs and different LD thresholds to increase the likelihood of finding the relevant genomic regions affecting the trait of interest. The main candidate genes identified, i.e., ATXN10, ADAM10, VAX2, ATP6V1B1, CRISPLD1, CAPRIN1, FA2H, SPEF2, PLXNA1, and CACNA2D3, are involved in important biological processes and metabolic pathways related to behavioral traits, social interactions, and aggressiveness in cattle. Future studies should further investigate the role of these candidate genes.

Frequency of genotypic markers for genetic disorders, colour, polledness, and major genes in Blanco Orejinegro cattle

The purpose of this work was to establish DNA marker frequencies for genetic disorders, colour, horned/polled trait, and major genes of importance for productive and reproductive traits in Blanco Orejinegro (BON) cattle. The Blanco Orejinegro breed is a Colombian creole breed characterized by a white hair coat on black skin with black ears, black hair on the middle part of the legs, and absence of horns. We genotyped 420 animals of Colombia, 70 with the GGP-Bovine 150 K chip and 350 with the GGP-Bovine 50 K chip. The markers were associated with 50 genetic diseases, 52 major gene variants related to productive traits, and 12 variants related to coat coloration, presence of horns, and adaptation, selected from the information contained in the chips. Genotype frequencies were estimated using the R statistical program. Genetic disorder annotations were derived using the Online Mendelian Inheritance in Animals tool (OMIA) and the average inbreeding coefficient (F) (n = 7799) using the MTDFNRM program. Carriers were found for 16 of the genetic disorders evaluated but with low frequencies (0.24 to 2.46%); no homozygous animals were found for the disorders. Carriers were associated with disorders such as bovine leukocyte adhesion deficiency (BLAD), deficiency of uridine monophosphate synthase (DUMPS), syndactyly, and epidermolysis bullosa (EB). The F was 4.41%. Concerning the genes associated with colour (TYR, MC1R, and PMEL), alleles related to black pigmentation, the absence of horns (polledness), and slick coat (an adaptive trait) were highly frequent (> 81.90%). Markers associated with milk production and quality, yellow fat, and fertility showed variable frequencies, indicating selection potential. Allele frequency of genetic disorders in BON cattle was low, suggesting few genetic disorder problems, with syndactyly being the most frequent condition. The markers associated with colour and polledness were almost fixed, with a frequency at or near 100%. Production and reproduction markers showed variability for selection.

Genome-wide association study of resistance/susceptibility to infectious bovine keratoconjunctivitis in Brazilian Hereford cattle

Genome-wide association studies were conducted to identify the more informative genomic regions and SNPs, as well as to identify candidate genes associated with infectious bovine keratoconjunctivitis (IBK) resistance/susceptibility in Hereford cattle. A Bayes B statistical approach was initially applied in genome-wide association studies by using deregressed estimated breeding values for IBK resistance/susceptibility. To estimate the combined effect of a genomic region that is potentially associated with QTL, 2504 non-overlapping 1-Mb windows that varied in SNP number were defined, with the most informative 24 windows including 427 SNPs and explaining more than 20% of the estimated genetic variance for IBK resistance/susceptibility. These regions were explored with respect to their biological functions through functional analysis to map potential candidate genes. The significant SNPs were mapped on chromosomes 1, 3, 5, 6, 7, 8, 10, 11, 12, 13, 14, 15, 18, 20, 23, and 28, and candidate genes were detected as related to the IBK. Most informative SNPs in term of genetic variance were located in proximity of genes related to phenotypic expression of lesions and biological processes associated to the IBK. Knowledge about phenotypic and genomic variation generated in the present study can be used to on design selection strategies to improve the resistance to IBK of Hereford cattle herds.

Evidence for and localization of proposed causative variants in cattle and pig genomes

BACKGROUND

This paper reviews the localization of published potential causative variants in contemporary pig and cattle reference genomes, and the evidence for their causality. In spite of the difficulties inherent to the identification of causative variants from genetic mapping and genome-wide association studies, researchers in animal genetics have proposed putative causative variants for several traits relevant to livestock breeding.

RESULTS

For this review, we read the literature that supports potential causative variants in 13 genes (ABCG2, DGAT1, GHR, IGF2, MC4R, MSTN, NR6A1, PHGK1, PRKAG3, PLRL, RYR1, SYNGR2 and VRTN) in cattle and pigs, and localized them in contemporary reference genomes. We review the evidence for their causality, by aiming to separate the evidence for the locus, the proposed causative gene and the proposed causative variant, and report the bioinformatic searches and tactics needed to localize the sequence variants in the cattle or pig genome.

CONCLUSIONS

Taken together, there is usually good evidence for the association at the locus level, some evidence for a specific causative gene at eight of the loci, and some experimental evidence for a specific causative variant at six of the loci. We recommend that researchers who report new potential causative variants use referenced coordinate systems, show local sequence context, and submit variants to repositories.

Heat stress response in slick vs normal-haired Criollo Limonero heifers in a tropical environment

To assess the effect of hair type on the heat stress response, 20 Criollo Limonero heifers with slick (n = 11) or normal hair (n = 9) were studied. Under a high temperature-humidity index (THI) environment, heat stress response was assessed through physiological variables that included respiration rate (RR), heart rate (HR), ruminal frequency (RMF), rectal temperature (RT), saliva pH (SPH), and lymphocyte count (LC) in the morning (5:00 AM, 27.4 °C, 64% relative humidity, THI = 77) and afternoon (1:00 PM, 34.5 °C, 70% relative humidity, THI = 88). A case-control study using a split plot design was used. Data were analyzed using ANOVA (PROC MIXED SAS 2010) and a statistical model comprising the fixed effects of hair length, sampling hour, interaction of hair length by sampling hour, and the random effect of animal nested within hair type on physiological variables associated with heat stress response. Sampling hour influenced (P < 0.0001) RR, RT, and (P < 0.003) SPH. Hair length influenced RR (P < 0.01) and RT (P < 0.04) and tended to influence LC (P < 0.07). The interaction of sampling hour by hair influenced RR (P < 0.04), RT (P < 0.0002), and both SPH and LC (P < 0.05). During afternoon hours, slick-haired heifers had lower values for RR (81 ± 4.2 vs 102 ± 4.7 bpm; P < 0.01), RT (39.5 ± 0.1 vs 40.3 ± 0.1 C°; P < 0.002), and LC (60 ± 3.2 vs 72.3 ± 3.6; P < 0.09) than normal-haired heifers. In normal-haired heifers, SPH increased during afternoon compared to morning-hours (8.66 ± 0.1 vs 9.11 ± 0.1; P < 0.04). It was concluded that slick-coated heifers exhibited an enhanced capability to cope with heat stress compared to normal-haired heifers likely due to an enhanced capacity for heat dissipation.

Genome-Wide Association Study and Transcriptome Differential Expression Analysis of the Feather Rate in Shouguang Chickens

The feather rate phenotype in chicks, including early-feathering and late-feathering phenotypes, are widely used as a sexing system in the poultry industry. The objective of this study was to obtain candidate genes associated with the feather rate in Shouguang chickens. In the present study, we collected 56 blood samples and 12 hair follicle samples of flight feathers from female Shouguang chickens. Then we identified the chromosome region associated with the feather rate by genome-wide association analysis (GWAS). We also performed RNA sequencing and analyzed differentially expressed genes between the early-feathering and late-feathering phenotypes using HISAT2, StringTie, and DESeq2. We identified a genomic region of 10.0–13.0 Mb of chromosome Z, which is statistically associated with the feather rate of Shouguang chickens at one-day old. After RNA sequencing analysis, 342 differentially expressed known genes between the early-feathering (EF) and late-feathering (LF) phenotypes were screened out, which were involved in epithelial cell differentiation, intermediate filament organization, protein serine kinase activity, peptidyl-serine phosphorylation, retinoic acid binding, and so on. The sperm flagellar 2 gene (SPEF2) and prolactin receptor (PRLR) gene were the only two overlapping genes between the results of GWAS and differential expression analysis, which implies that SPEF2 and PRLR are possible candidate genes for the formation of the chicken feathering phenotype in the present study. Our findings help to elucidate the molecular mechanism of the feather rate in chicks.

Homozygosity Mapping Reveals Population History and Trait Architecture in Self-Incompatible Pear (Pyrus spp.)

Runs of homozygosity (ROH) have been widely used to study population history and trait architecture in humans and livestock species, but their application in self-incompatible plants has not been reported. The distributions of ROH in 199 accessions representing Asian pears (45), European pears (109), and interspecific hybrids (45) were investigated using genotyping-by-sequencing in this study. Fruit phenotypes including fruit weight, firmness, Brix, titratable acidity, and flavor volatiles were measured for genotype-phenotype analyses. The average number of ROH and the average total genomic length of ROH were 6 and 11 Mb, respectively, in Asian accessions, and 13 and 30 Mb, respectively, in European accessions. Significant associations between genomic inbreeding coefficients (F_ROH) and phenotypes were observed for 23 out of 32 traits analyzed. An overlap between ROH islands and significant markers from genome-wide association analyses was observed. Previously published quantitative trait loci for fruit traits and disease resistances also overlapped with some of the ROH islands. A prominent ROH island at the bottom of linkage group 17 overlapped with a recombination-supressed genomic region harboring the self-incompatibility locus. The observed ROH patterns suggested that systematic breeding of European pears would have started earlier than of Asian pears. Our research suggest that F_ROH would serve as a novel tool for managing inbreeding in gene-banks of self-incompatible plant species. ROH mapping provides a complementary strategy to unravel the genetic architecture of complex traits, and to evaluate differential selection in outbred plants. This seminal work would provide foundation for the ROH research in self-incompatible plants.

Genome-wide detection of signatures of selection in indicine and Brazilian locally adapted taurine cattle breeds using whole-genome re-sequencing data

BACKGROUND

The cattle introduced by European conquerors during the Brazilian colonization period were exposed to a process of natural selection in different types of biomes throughout the country, leading to the development of locally adapted cattle breeds. In this study, whole-genome re-sequencing data from indicine and Brazilian locally adapted taurine cattle breeds were used to detect genomic regions under selective pressure. Within-population and cross-population statistics were combined separately in a single score using the de-correlated composite of multiple signals (DCMS) method. Putative sweep regions were revealed by assessing the top 1% of the empirical distribution generated by the DCMS statistics.

RESULTS

A total of 33,328,447 biallelic SNPs with an average read depth of 12.4X passed the hard filtering process and were used to access putative sweep regions. Admixture has occurred in some locally adapted taurine populations due to the introgression of exotic breeds. The genomic inbreeding coefficient based on runs of homozygosity (ROH) concurred with the populations' historical background. Signatures of selection retrieved from the DCMS statistics provided a comprehensive set of putative candidate genes and revealed QTLs disclosing cattle production traits and adaptation to the challenging environments. Additionally, several candidate regions overlapped with previous regions under selection described in the literature for other cattle breeds.

CONCLUSION

The current study reported putative sweep regions that can provide important insights to better understand the selective forces shaping the genome of the indicine and Brazilian locally adapted taurine cattle breeds. Such regions likely harbor traces of natural selection pressures by which these populations have been exposed and may elucidate footprints for adaptation to the challenging climatic conditions.

Prospects for gene introgression or gene editing as a strategy for reduction of the impact of heat stress on production and reproduction in cattle

In cattle, genetic variation exists in regulation of body temperature and stabilization of cellular function during heat stress. There are opportunities to reduce the impact of heat stress on cattle production by identifying the causative mutations responsible for genetic variation in thermotolerance and transferring specific alleles that confer thermotolerance to breeds not adapted to hot climates. An example of a mutation conferring superior ability to regulate body temperature is the group of frame-sift mutations in the prolactin receptor gene (PRLR) that lead to a truncated receptor and development of cattle with a short, sleek hair coat. Slick mutations in PRLR have been found in several extant breeds derived from criollo cattle. The slick mutation in Senepol cattle has been introgressed into dairy cattle in Puerto Rico, Florida and New Zealand. An example of a mutation that confers cellular protection against elevated body temperature is a deletion mutation in the promoter region of a heat shock protein 70 gene called HSPA1L. Inheritance of the mutation results in amplification of the transcriptional response of HSPA1L to heat shock and increased cell survival. The case of PRLR provides a promising example of the efficacy of the genetic approach outlined in this paper. Identification of other mutations conferring thermotolerance at the whole-animal or cellular level will lead to additional opportunities for using genetic solutions to reduce the impact of heat stress.

Haplotype-Based Genome-Wide Association Study and Identification of Candidate Genes Associated with Carcass Traits in Hanwoo Cattle

Hanwoo, is the most popular native beef cattle in South Korea. Due to its extensive popularity, research is ongoing to enhance its carcass quality and marbling traits. In this study we conducted a haplotype-based genome-wide association study (GWAS) by constructing haplotype blocks by three methods: number of single nucleotide polymorphisms (SNPs) in a haplotype block (nsnp), length of genomic region in kb (Len) and linkage disequilibrium (LD). Significant haplotype blocks and genes associated with them were identified for carcass traits such as BFT (back fat thickness), EMA (eye Muscle area), CWT (carcass weight) and MS (marbling score). Gene-set enrichment analysis and functional annotation of genes in the significantly-associated loci revealed candidate genes, including PLCB1 and PLCB4 present on BTA13, coding for phospholipases, which might be important candidates for increasing fat deposition due to their role in lipid metabolism and adipogenesis. CEL (carboxyl ester lipase), a bile-salt activated lipase, responsible for lipid catabolic process was also identified within the significantly-associated haplotype block on BTA11. The results were validated in a different Hanwoo population. The genes and pathways identified in this study may serve as good candidates for improving carcass traits in Hanwoo cattle.

A Genetic Investigation of Island Jersey Cattle, the Foundation of the Jersey Breed: Comparing Population Structure and Selection to Guernsey, Holstein, and United States Jersey Cattle

For two centuries, Jersey cattle were exported globally, adapting to varying climates and production systems, yet the founding population remained genetically isolated on the Island of Jersey. The Island of Jersey formally allowed the importation of pure Jersey cattle in 2008. This study characterized the genetic variation of 49 popular bulls from the Island of Jersey born from 1964 to 2004 and compared them to 47 non-Island Jersey bulls and cows, primarily from the United States In addition, 21 Guernsey cattle derived from the Island of Guernsey and 71 Holstein cattle served as reference populations for genetic comparison. Cattle were genotyped on the Illumina BovineHD Beadchip producing 777,962 SNPs spanning the genome. Principal component analysis revealed population stratification within breed reflective of individual animal’s continental origin. When compared to Holstein and Guernsey, all Jersey clustered together by breed. The Jersey breed demonstrated increased inbreeding in comparison to Holstein or Guernsey with slightly higher estimates of inbreeding coefficients and identity-by-descent. The Island and United States Jersey have relatively similar, yet statistically different inbreeding estimates despite vastly different population sizes and gene flow. Signatures of selection within Island Jersey were identified using genome-wide homozygosity association and marker-based F_ST that provided population informative single-nucleotide polymorphism (SNPs). Biological significance of the homozygosity association results identified multiple genes on chromosomes 5, 24, and 27, involved in immune function and cellular processes. Overall, genomic variation was identified between the Island and non-Island Jersey cattle producing population informative SNPs and differing runs of homozygosity (ROH) over immune regulation and metabolic genes. Results on inbreeding measures and ROH may reflect varying effective population size or differential selection with grazing systems promoting natural selection for traits such as parasite resistance, whereas confinement systems demonstrate a more intensive artificial selection. More broadly, differences in breed formation, particularly between the two Channel Island breeds, likely contributed to the variation in ROH and inbreeding. This research provides a reference for the Jersey breed based on the genetic foundation of the Island cattle as compared to the intensively selected United States cattle, and identifies regions of the genome for future investigation of immune regulation and metabolic processes.

Gene expression profiles in bovine granulocytes reflect the aberration of liver functions

Liver performs several important functions; however, predicting its functions is difficult. Methods of analyzing gene expression profiles, for example, microarray, provide functional information of tissues. Liver and peripheral blood leukocytes (PBLs) were collected from Holstein cows subjected to two different physiological conditions (non-pregnant and pregnant), and PBLs were fractionated by gradient cell separation. RNA from PBLs and liver were applied to oligo-DNA microarray and reverse transcription-quantitative polymerase chain reaction (RT-qPCR). It revealed a group of stable bovine liver genes under constant physiological conditions. When they applied to physiological conditions including non-pregnant and pregnant, the profiles of some genes in liver were consistent with those in PBLs. Microarray data subjected to a principal component analysis (PCA) showed that the hepatic gene expression profiles were more consistent with those of granulocytes than mononuclear cells. The relationship of gene profiles in liver with granulocytes was confirmed by RT-qPCR and hierarchical cluster analysis. Gene profiles of granulocytes were more reliable than those of mononuclear cells, which reflected liver functions. These results suggest that the genes expressed in PBLs, particularly granulocytes, may be convenient bioindicators for the diagnosis of clinical disorder and/or detecting aberration of liver functions in cows subjected to different physiological conditions.

Genetic Selection for Thermotolerance in Ruminants

Variations in climatic variables (temperature, humidity and solar radiation) negatively impact livestock growth, reproduction, and production. Heat stress, for instance, is a source of huge financial loss to livestock production globally. There have been significant advances in physical modifications of animal environment and nutritional interventions as tools of heat stress mitigation. Unfortunately, these are short-term solutions and may be unsustainable, costly, and not applicable to all production systems. Accordingly, there is a need for innovative, practical, and sustainable approaches to overcome the challenges posed by global warming and climate change-induced heat stress. This review highlights attempts to genetically select and breed ruminants for thermotolerance and thereby sustain production in the face of changing climates. One effective way is to incorporate sustainable heat abatement strategies in ruminant production. Improved knowledge of the physiology of ruminant acclimation to harsh environments, the opportunities and tools available for selecting and breeding thermotolerant ruminants, and the matching of animals to appropriate environments should help to minimise the effect of heat stress on sustainable animal genetic resource growth, production, and reproduction to ensure protein food security.

About the inconvenience of handling mixed-breed herds; aspects of social behavior as a potential source of stress and economic losses

The objective of this study was to assess aspects of the social behavior of a mixed-breed herd of beef cows as a potential source for stress and economic losses. Angus (AN; N = 10), Brahman (BR; N = 10), and Senepol (SE; N = 10) cows were assigned to two groups (N = 15 each containing equal breed numbers) on separate pastures. Agonistic interactions (win/loss) during feeding were recorded daily for 45 days. Dominance values were estimated as the proportion of individuals dominated to total herdmates. From this, individuals were placed into social categories based upon linear ranking as follows: dominants (D), intermediate (I), and subordinates (S). Breed influenced (P < 0.01) social category, with SE cows being dominants (P < 0.05) over AN and BR cows. Interactions between AN and BR cows were less (P < 0.0005) than interactions between AN and SE (53 vs 140, respectively). Within breeds, BR (152) and SE (182) cows had more (P < 0.0005) agonistic interactions than AN (107) cows. Although apparently influenced by breed, agonistic interactions occurred more frequently (P < 0.005) between social categories than within social categories (814 vs 310, respectively). Dominant cows were involved in more agonistic interactions with cows from different social categories than were intermediate and subordinate cows (P < 0.0005). However, intermediate (100) and subordinate (157) cows generated more (P < 0.0005) agonistic interactions within their own social category than dominant cows (53). It was concluded that, in mixed-breed herds, breed influences both social organization and agonistic interactions which could be considered as potential sources of stress and economic losses.

Demography and rapid local adaptation shape Creole cattle genome diversity in the tropics

The introduction of Iberian cattle in the Americas after Columbus' arrival imposed high selection pressures on a limited number of animals over a brief period of time. Knowledge of the genomic regions selected during this process may help in enhancing climatic resilience and sustainable animal production. We first determined taurine and indicine contributions to the genomic structure of modern Creole cattle. Second, we inferred their demographic history using approximate Bayesian computation (ABC), linkage disequilibrium (LD) and N e Slope (NeS) analysis. Third, we performed whole genome scans for selection signatures based on cross-population extended haplotype homozygosity (XP-EHH) and population differentiation (F ST) to disentangle the genetic mechanisms involved in adaptation and phenotypic change by a rapid and major environmental transition. To tackle these questions, we combined SNP array data (~54,000 SNPs) in Creole breeds with their modern putative Iberian ancestors. Reconstruction of the population history of Creoles from the end of the 15th century indicated a major demographic expansion until the introduction of zebu and commercial breeds into the Americas ~180 years ago, coinciding with a drastic N e contraction. NeS analysis provided insights into short-term complexity in population change and depicted a decrease/expansion episode at the end of the ABC-inferred expansion, as well as several additional fluctuations in N e with the attainment of the current small N e only towards the end of the 20th century. Selection signatures for tropical adaptation pinpointed the thermoregulatory slick hair coat region, identifying a new candidate gene (GDNF), as well as novel candidate regions involved in immune function, behavioural processes, iron metabolism and adaptation to new feeding conditions. The outcomes from this study will help in future-proofing farm animal genetic resources (FAnGR) by providing molecular tools that allow selection for improved cattle performance, resilience and welfare under climate change.

On Estimation of Genome Composition in Genetically Admixed Individuals Using Constrained Genomic Regression

Quantifying the population stratification in genotype samples has become a standard procedure for data manipulation before conducting genome wide association studies, as well as for tracing patterns of migration in humans and animals, and for inference about extinct founder populations. The most widely used approach capable of providing biologically interpretable results is a likelihood formulation which allows for estimation of founder genome proportions and founder allele frequency conditional on the observed genotypes. However, if founder allele frequencies are known and samples are dominated by admixed genotypes this approach may lead to biased inference. In addition, processing time increases drastically with the number of genetic markers. This article describes a simplified approach for obtaining biologically meaningful measures of population stratification at the genotype level conditional on known founder allele frequencies. It was tested on cattle and human data sets with 4,022 and 150,000 genetic markers, respectively, and proved to be very accurate in situations where founder poplations were correctly specified, or under-, over-, and miss-specified. Moreover, processing time was only marginally affected by an increase in the number of markers.

Convergent Evolution of Slick Coat in Cattle through Truncation Mutations in the Prolactin Receptor

Evolutionary adaptations are occasionally convergent solutions to the same problem. A mutation contributing to a heat tolerance adaptation in Senepol cattle, a New World breed of mostly European descent, results in the distinct phenotype known as slick, where an animal has shorter hair and lower follicle density across its coat than wild type animals. The causal variant, located in the 11^th exon of prolactin receptor, produces a frameshift that results in a truncated protein. However, this mutation does not explain all cases of slick coats found in criollo breeds. Here, we obtained genome sequences from slick cattle of a geographically distinct criollo breed, namely Limonero, whose ancestors were originally brought to the Americas by the Spanish. These data were used to identify new causal alleles in the 11^th exon of the prolactin receptor, two of which also encode shortened proteins that remove a highly conserved tyrosine residue. These new mutations explained almost 90% of investigated cases of animals that had slick coats, but which also did not carry the Senepol slick allele. These results demonstrate convergent evolution at the molecular level in a trait important to the adaptation of an animal to its environment.

BREEDING AND GENETICS SYMPOSIUM:Breeding heat tolerant dairy cattle: the case for introgression of the "slick" prolactin receptor variant into dairy breeds

Increasing environmental temperatures are a threat to the sustainability of livestock production and, because of the high metabolic demands of lactation, to dairy production in particular. Summer heat waves in temperate climates reduce feed intake, milk production, and cow comfort. In extreme heat events, there is an increase in cow mortality. In tropical climates, dairy cattle are mostly (zebu) type or zebu crossbred with temperate dairy breeds. Crossbreeding is undertaken to combine the heat tolerance and tick resistance of zebu with the productivity of temperate dairy breeds. In the absence of improved heat tolerance, milk production and fertility of temperate cattle is severely impaired. We have recently identified a key role for the prolactin pathway in regulating heat tolerance. A de novo mutation in prolactin that impairs prolactin activity was discovered in hairy and heat intolerant, New Zealand dairy cattle. The phenotypes produced were remarkably similar to those seen in fescue toxicosis, a syndrome seen in grazing cattle in the U.S. where ingestion of ergovaline, a fungal toxin from infected pasture, inhibits prolactin secretion. Recognition of the role of prolactin in hairy cattle led us to identify a deletion in exon 10 of the long-form of the prolactin receptor in Senepol cattle that causes truncation of the protein and determines the slick coat and heat tolerance traits found in this , beef breed. The short form of the prolactin receptor is predicted to be unaffected by the deletion. Knowledge of this dominant mutation has provided the impetus to begin a crossbreeding program to investigate performance and heat tolerance of temperate dairy cattle carrying the slick, prolactin receptor variant. The perceived opportunity is to introgress this variant into temperate dairy cattle to enable performance and welfare improvement in hot climates. Heat tolerance of cattle with slick coats appears to be mostly associated with coat type although sweating ability may also be enhanced. Further investigation is required of performance traits in cows homozygous for the slick variant because the published data are almost exclusively from heterozygous animals. Combination of the slick mutation with other favorable genes for heat tolerance, especially those for coat color, will be particularly enabled by gene editing technologies, offering opportunities for further improvement in bovine thermotolerance.

Runs of homozygosity: current knowledge and applications in livestock

This review presents a broader approach to the implementation and study of runs of homozygosity (ROH) in animal populations, focusing on identifying and characterizing ROH and their practical implications. ROH are continuous homozygous segments that are common in individuals and populations. The ability of these homozygous segments to give insight into a population's genetic events makes them a useful tool that can provide information about the demographic evolution of a population over time. Furthermore, ROH provide useful information about the genetic relatedness among individuals, helping to minimize the inbreeding rate and also helping to expose deleterious variants in the genome. The frequency, size and distribution of ROH in the genome are influenced by factors such as natural and artificial selection, recombination, linkage disequilibrium, population structure, mutation rate and inbreeding level. Calculating the inbreeding coefficient from molecular information from ROH (F_ROH ) is more accurate for estimating autozygosity and for detecting both past and more recent inbreeding effects than are estimates from pedigree data (F_PED ). The better results of F_ROH suggest that F_ROH can be used to infer information about the history and inbreeding levels of a population in the absence of genealogical information. The selection of superior animals has produced large phenotypic changes and has reshaped the ROH patterns in various regions of the genome. Additionally, selection increases homozygosity around the target locus, and deleterious variants are seen to occur more frequently in ROH regions. Studies involving ROH are increasingly common and provide valuable information about how the genome's architecture can disclose a population's genetic background. By revealing the molecular changes in populations over time, genome-wide information is crucial to understanding antecedent genome architecture and, therefore, to maintaining diversity and fitness in endangered livestock breeds.

A simple strategy for managing many recessive disorders in a dairy cattle breeding program

Background

High-density single nucleotide polymorphism genotypes have recently been used to identify a number of novel recessive mutations that adversely affect fertility in dairy cattle, as well as to track other conditions such as red coat color and polled. Most current methods for mate allocation fail to consider this information, and it will become increasingly difficult to manage matings as the number of recessive mutations to be accounted for increases.

Methods

A modified version of a mating strategy that constrains inbreeding based on genomics (the Pryce method) was developed that also accounts for the economic effects of Mendelian disorders on overall economic merit (modified Pryce method) and compared with random mating, truncation selection, and the Pryce scheme. Several scenarios were considered, including scenarios with six hypothetical recessive alleles and 12 recessive alleles that are currently segregating in the US Holstein population.

Results

The Pryce method and the modified Pryce method showed similar ability to reduce frequencies of recessive alleles, particularly for loci with frequencies greater than 0.30. The modified Pryce method outperformed the Pryce method for low-frequency alleles with small economic value. Cumulative genetic gain for the selection objective was slightly greater when using the Pryce method, but rates of inbreeding were similar across methods.

Conclusions

The proposed method reduces allele frequencies faster than other methods, and also can be used to maintain or increase the frequency of desirable recessives. It can be easily implemented in software for mate allocation, and the code used in this study is freely available as a reference implementation.

Challenges and opportunities in genetic improvement of local livestock breeds

Sufficient genetic variation in livestock populations is necessary both for adaptation to future changes in climate and consumer demand, and for continual genetic improvement of economically important traits. Unfortunately, the current trend is for reduced genetic variation, both within and across breeds. The latter occurs primarily through the loss of small, local breeds. Inferior production is a key driver for loss of small breeds, as they are replaced by high-output international transboundary breeds. Selection to improve productivity of small local breeds is therefore critical for their long term survival. The objective of this paper is to review the technology options available for the genetic improvement of small local breeds and discuss their feasibility. Most technologies have been developed for the high-input breeds and consequently are more favorably applied in that context. Nevertheless, their application in local breeds is not precluded and can yield significant benefits, especially when multiple technologies are applied in close collaboration with farmers and breeders. Breeding strategies that require cooperation and centralized decision-making, such as optimal contribution selection, may in fact be more easily implemented in small breeds.

Genomic data as the "hitchhiker's guide" to cattle adaptation: tracking the milestones of past selection in the bovine genome

The bovine species have witnessed and played a major role in the drastic socio-economical changes that shaped our culture over the last 10,000 years. During this journey, cattle "hitchhiked" on human development and colonized the world, facing strong selective pressures such as dramatic environmental changes and disease challenge. Consequently, hundreds of specialized cattle breeds emerged and spread around the globe, making up a rich spectrum of genomic resources. Their DNA still carry the scars left from adapting to this wide range of conditions, and we are now empowered with data and analytical tools to track the milestones of past selection in their genomes. In this review paper, we provide a summary of the reconstructed demographic events that shaped cattle diversity, offer a critical synthesis of popular methodologies applied to the search for signatures of selection (SS) in genomic data, and give examples of recent SS studies in cattle. Then, we outline the potential and challenges of the application of SS analysis in cattle, and discuss the future directions in this field.

Functionally reciprocal mutations of the prolactin signalling pathway define hairy and slick cattle

Lactation, hair development and homeothermy are characteristic evolutionary features that define mammals from other vertebrate species. Here we describe the discovery of two autosomal dominant mutations with antagonistic, pleiotropic effects on all three of these biological processes, mediated through the prolactin signalling pathway. Most conspicuously, mutations in prolactin (PRL) and its receptor (PRLR) have an impact on thermoregulation and hair morphology phenotypes, giving prominence to this pathway outside of its classical roles in lactation.

The hormone prolactin is a known modulator of mammalian lactation and hair growth. Here, the authors describe two dominant mutations in bovine prolactin and its receptor, demonstrating antagonistic effects on these traits and highlighting a role for this pathway in sweat gland function and thermoregulation.

The SLICK hair locus derived from Senepol cattle confers thermotolerance to intensively managed lactating Holstein cows

The SLICK haplotype (http://omia.angis.org.au/OMIA001372/9913/) in cattle confers animals with a short and sleek hair coat. Originally identified in Senepol cattle, the gene has been introduced into Holsteins. The objectives of the current study were to determine (1) whether lactating Holsteins with the slick hair phenotype have superior ability for thermoregulation compared with wild-type cows or relatives not inheriting the SLICK haplotype, and (2) whether seasonal depression in milk yield would be reduced in SLICK cows. In experiment 1, diurnal variation in vaginal temperature in the summer was monitored for cows housed in a freestall barn with fans and sprinklers. Vaginal temperatures were lower in slick-haired cows than in relatives and wild-type cows. In experiment 2, acute responses to heat stress were monitored after cows were moved to a dry lot in which the only heat abatement was shade cloth. The increases in rectal temperature and respiration rate caused by heat stress during the day were lower for slick cows than for relatives or wild-type cows. Moreover, sweating rate was higher for slick cows than for cows of the other 2 types. In experiment 3, effects of season of calving (summer vs. winter) on milk yield and composition were determined. Compared with milk yield of cows calving in winter, milk yield during the first 90 d in milk was lower for cows calving in the summer. However, this reduction was less pronounced for slick cows than for wild-type cows. In conclusion, Holsteins with slick hair have superior thermoregulatory ability compared with non-slick animals and experience a less drastic depression in milk yield during the summer.